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BLOCK - I: 
STATISTICAL THERMODYNAMICS 


Unit-1: Statistical thermodynamics 
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1.0 Introduction 


In thermodynamics, statistical thermodynamics is the study of the 
microscopic behaviors of thermodynamic systems using probability theory. 
Statistical thermodynamics, generally, provides a molecular level 
interpretation of thermodynamic quantities such as work, heat, free energy, 
and entropy. Statistical thermodynamics was born in 1870 with the work of 
Austrain physicist Ludwig Boltzmann, much of which was collectively 
published in Boltzmann's 1896 Lectures on Gas Theory. 


1.4 Objectives 


After going through this unit, you will be able to: 
e Understand about the statistical thermodynamics 
e Understand the Microscopic behaviour of thermodynamic system 
using probability theory. 
e Explain the concept of Negative Kelvin Temperature 
e Explains the microstates and configuration 


1.2 Maxwell - Boltzmann Distribution 


Postulates: 


1. The system considered is an isolated system of independent, non 
interacting identical particles without spin in equilibrium at a 
definite temperature. 

2. The particles are distinguishable 


Unit -1 Statistical 
Thermodynamics 


NOTE 


Self-Instructional Material 


a 4 There is no restriction to the number of particles 
Unit -1 Statistical which may occupy any given quantum state. Let us consider 


Thermodynamics : t . . CM 
any ideal gas containing N identical distinguishable 
independent particles. let at any instant 

NOTE no particle have energy Eo 


ni particles have energy Ei 


ni particles have energy Ei 


Total number of particles in the system is 


No +N n2 + N3+ ---------- ni + --- =N 


Total number of particles in the system is assumed to 
be constant and hence 


$an, 20 > (1) 
= 


This is the first condition of constraint. 


4. Total energy possessed by all the molecules E is 


E=Em + En, * * En, * => Eg, 


i-0 


5. Total energy is constant and hence OE = 0 


Self-Instructional Material 


S Ean, =0 > (2) 
1=0 


This is the second condition of constraint. 


6. he system is considered at statistical equilibrium and hence 
thermodynamic probability will be maximum. 


1.2.1 Maximization of thermodynamic probability 
For Non - degenerate energy levels: 


The probability of distribution of particles among the quantum states is 
proportional to the number of different ways in which the molecules can be 
arranged under the specified conditions. The number of ways in which ‘N’ 
distinguishable particles can be arranged in different energy levels such that no 
particles are in the ground state , n, in the 1", no in the 2™ and so is given by 


N! 
W= — x cons tan f 


m;! 
i=0 


In W = In N!-In( zn !) + cons tan f 
i=0 


Applying stirlings approximation formula 


In Ni NinN-N 


In W=NinN-N- > In( n;!) + cons tan t 
1=0 


=NInN-N-))(n,Inn,—n,)+ constant 
i-Ü 


oo oo 


- NI N-N -Y nn, Yn, * constant 


i=0 i=0 
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Self-Instructional Material 


Unit -1 Statistical 
Thermodynamics 


NOTE 


Self-Instructional Material 


When the system attains statistical equilibrium there will be 
most probable distribution and thermodynamic probability will 
be maximum and hence variation in W 


1. O0 or Olno is zero. 
alno =0 


din 9—O0(NIn N) - 0 n,n n,) &constant 
i=0 


-0- 5 (n +4innyan 20 


i=0 un 
=— X 'In( n;)ón;- 0 >> an, =0 
i=0 i=0 
(or) 


nC n,)0n;2 0 > (3) 
i=0 


This is the 3" condition of constraint. For a system of maximum 
thermodynamic probability and for a stem in statistical equilibrium 
condition 1, 2&3 must be satisfied. 


When there are more than two conditions, lagrangian method of 
undetermined multiplers is used to solve the problem. 


Equation (1) is multiplied by a, equation (2) is multiplied by B and 
combined with equation (3) we get 


ay ðn, + B» Eon, t »» n,)on, - 0 
i-0 i=0 i=0 


oo 


S (æ+ BE,+ In n,)ðn,=0 


i=0 


Since, the variations Oni, On», On3 etc are independent of each other, 
provided the conditions (1) & (2) are satisfied, the co- efficient of On1’s 
must be zero and On;zO 


a+ BE;* In n;- 0 Unit -1 Statistical 


Thermodynamics 
Inn;= —a-— DE, NOTE 
z4 1l 
A= e^ PE (a= ? p= ) 
RT KT 
n, = co FE — et "T o E, KT 


-EKT 


The factor e is called Boltzmann factor and R is the Boltzmann 


constant, the expression 


n, =% e  =exp™ exp is called Maxwell -— Boltzmann 


distribution law for the system having N distinguishable particles 
distributed in non degenerated energy levels. 


The total energy or any individual form of energy is given as 


Ecrotal) = Eelec +Evibt+ Erot + Etrans 


1.2.2 When the energy levels are degenerate 

The number of ways in which ‘N’ distinguishable particles can be 
arranged in different energy levels having the degeneracy go, n, in the first 
state having the degeneracy gi, n» in the 3" having the degeneracy g2 and 
so on i.e. given by 


7 n 
WENTE X cons tan t 
i=0 n.! 


rj 


oo n. 
In W= In( M7 g, —X cons tan t) 
ES nj 


= ni 

=InN!+In 2 -h =n !xcons tan t 

T. WE. 
i-0 i=0 


=In N+% n,n gi -J hn n;!+ constant 


i=0 i=0 


Self-Instructional Material 


| Applying Stirlings approximation Formula for evaluation 
Unit -1 Statistical 
Thermodynamics of factorials of large number we get 


NOTE In Ni NinN-N 


co oo 


lnW-2NhN-N-«Ynn 8, >, (n,n n, —n,)+ constant 


i=0 i=0 


oo oo oo 


=NinN-N+ > nng,- nnn, + >on, + constant 


i=0 i=0 i=0 


=NinN+ > n,n g, > nin n, +constant 


i=0 i=0 


When the system attains statistical equilibrium there will be 
most probable distribution and thermodynamic probability will 
be maximum, i.e. variation in law will be zero. 


oW -0 
(or) dln W 20 


9In W 2 ON In N - 0? n, In g, -90» n, hnn, + O(cons tan t) 2 0 


i=0 i=0 


=0+ J In gàn,- 9 (14 Inn, )on,+ 0=0 


i=0 i=0 


= Sin g,on,— n non;- 0 
i-0 i=0 


= Yn 2; /ni on, = 0 
i-0 


Self-Instructional Material (or) 


— : n.)on.=0 (3 Unit -1 Statistical 
2 8:/ iJon, (3) Thermodynamics 


This is the 3 condition of constraint 


NOTE 
For a system of maximum thermodynamic probability and for a 


system in statistical equilibrium conditions 1, 2 & 3 must be satisfied 


When there are more than two conditions Lagrange’s methods of 
undetermined multipliers is used to solve the problem. 


Equation (1) is multiplied by a, equation (2) by B & combined with equation 
(3) 


a Y on + BFE ðn- Y]n(550n =0 
i i i n. t 


i-0 i=0 i=0 i 
S (a+ BE,- Ww g, jn)0n,- 0 
i=0 
Since the variations Oni, On», On3 etc are independent of each other, provided 


the conditions 1 &3 are satisfied the coefficient of Oni's must be equal to zero. 
Oni Z0. 


at DE, -In( g/ nj) 20 
In( g;/n;) =a+ BE, 


a+ BE; 


gin =e 


= per ere 


Where, gi is the degree of degeneracy or statistical weight factor 
for the i" level. This expression is called Maxwell Boltzmann 
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Self-Instructional Material 


| 
Unit -1 Statistical 
Thermodynamics 


NOTE 


Self-Instructional Material 


distribution law. It gives the number of particles in the i^ level 
with degeneracy gi having energy Ej and 


Bi(Total) = Eelec +Evibt+ Erot + Erans 


1.3 Negative Kelvin Temperature 


According to Maxwell Boltzmann statistics the equilibrium 
population of two levels is given by 


Pi olg Age 


No 


Where nj is the number of particles in the upper i^ level and no is 
the number of particles in the lowest level at temperature T and Eiis 
the energy of the i" level in excess of zero point level. At T = 0 
ni/no=0 that is all the molecules will be in the ground state as T—^0o 
ni/No=1 and the population of the states will be equal. 


In general in any system at statistical equilibrium the higher 
level will be occupied by lesser number of particles than the lower 
ones or will be equally occupied. To reverse the population ratios 
and have higher level more occupied than the lower levels, 
temperatures even higher than infinity are needed. A population 
ratio larger than one would require T<o or negative. This 
temperature at which population inversion occurs is called negative 
absolute temperature. This temperature must be beyond infinity. 


To understand the phenomenon of negative absolute temperature, 
consider a system of ‘N’ particles existing only in two energy 
levels. i.e in zero and 1“ level. The entropy ‘S’ when plotted as a 
function of energy the following curve is obtained. 


Unit -1 Statistical 
Thermodynamics 


* NOTE 


Entropy (S) 


0 m Ne 


Energy (E) 


At zero energy all the *N' particles are in the zero level (lowest) which the 


state of minimum disorder to minimum entropyis( ) <=1 jT 
QE" 


As the temperature is increased or energy is supplied to the system 
population in the upper level increases. When the two energy levels are 
equally populated the internal energy of the system will correspond to 
NE/2. There is maximum disorder and maximum entropy. It is the most 
probable state and OS/OE =1/T =0 as To or more energy is supplied all 
the ‘N ‘particles will be in the upper level. EZNrz. This is a state of 
maximum energy and minimum entropy with minimum disorder. The left 


half of the curve has a positive slope i.e e» = |/T is positive that is T is 


os 
+ve. The right half of the curve as negative slope ie = )T=-ve that 


is T= -ve. This is the region of negative absolute temperature. Negative 
absolute temperatures are defined by the slope of S vs E curves. 


In the figure as we proceed from left to right in the direction of 
increasing energy there is increasing hotness and therefore increase in 
entropy. At the position of maximum entropy where both the energy levels 
are equally populated, the temperature is infinite. Beyond the maximum the 
temperature must be hotter than infinity. Hence negative absolute 
temperatures are hotter than infinity. 


1.4 Check Your Progress 

1. What is meant by statistical therrmodynamics? 

2. Write the Maxwell — Boltzmann distribution law for the system having 
N distinguishable particles distributed in non degenerated energy levels. 
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| 
Unit -1 Statistical 
Thermodynamics 


NOTE 


Self-Instructional Material 


3. Give the expression for Maxwell distribution law of 
molecular energies. 

4. What are the relative populations of the states of a two-level 
system when the temperature is infinite? 

5. What are the relative populations of the states of a two-level 
system as 

the temperature approaches zero? 


1.5 Answers to Check Your Progress Questions 


1. Statistical thermodynamics is the study of the microscopic 
behaviors of thermodynamic systems using probability theory. 
Statistical thermodynamics, generally, provides a molecular 
level interpretation of thermodynamic quantities such 
as Work, heat, free energy, and entropy. 


(PE) is called Maxwell — 


2. n=e*.e™ =exp™ exp 
Boltzmann distribution law for the system having N 
distinguishable particles distributed in non degenerated energy 


levels. 


z4 d 
inc = g gob (a=—~, p=) 
S a i RT KT 


Where, giis the degree of degeneracy or statistical weight factor 
for the i" level. This expression is called Maxwell Boltzmann 
distribution law. 

4. All the molecules will be in the ground state as T—^oo ni/no =1 
and the population of the states will be equal. 


5. At T = 0 ni/no=0. 


1.6 Summary 


° In thermodynamics, statistical thermodynamics is the 
study of the microscopic behaviors of thermodynamic 
systems using probability theory. Statistical thermodynamics, 
generally, provides a molecular level interpretation of 
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thermodynamic quantities such as work, heat, free energy, 
and entropy 

° Maxwell-Boltzmann statistics gives the average 
number of particles found in a given single-particle microstate. 

° The energies of such particles follow what is known 
as Maxwell-Boltzmann statistics, and the 
statistical distribution of speeds is derived by equating particle 
energies with kinetic energy. 

. One very important conclusion that will emerge from 
the following analysis is that the populations of states depend on 
a single parameter, the ‘temperature’. That is, statistical 
thermodynamics provides a molecular justification for the 
concept of temperature and some insight into this crucially 
important quantity. 


1.7 Keywords 


l. 


. Statistical thermodynamics is the study of the microscopic 


behaviors of thermodynamic systems using probability theory. 


2: 


The weight of a configuration is the number of ways that molecules 


can be distributed over the available states. 


3. 


The most probable distribution, that of the greatest weight, is the 


Boltzmann distribution. 


1.8  Self-assessment questions and exercises 


1. Derive Maxwell distribution law of molecular energies. 
2. Explain Negative Kelvin Temperature. 


1.9 Further readings 


1. 


Statistical Thermodynamics, M.C. Gupta, Wiley Eastern, New 


Delhi, 1990. 


2. 


Introduction to Statistical Thermodynamics, R.P.H.Gasser and 


W.G.Richards, World Scientific, Singapore, 1995. 


3. 


Physical chemistry, Peter Atkins, Julio De Paula, ght edition, 2010. 
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NOTE 
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Unit- II Partition Function 

Structure 

2.0 Introduction 

2.1 Objectives 

2.2 Partition Function 
2.2.1 Vibrtional partiion function 
2.2.2 Rotational partition function 
2.2.3 Electronic partition function 

2.3 Thermodynamic functions in terms of partition function 
2.3.1 Internal energy from partition function 
2.3.2 Heat capacity from partition function 
2.3.3 Entropy from partition function 
2.3.4 Pressure from partition function 
2.3.5 Third law and partition funnction 
2.3.6 Helmholtz Free Enery from Partition Function 
2.3.7 Enthalpy from partition function 
2.3.8 Gibbs free energy from partition function 

24 Check Your Progress 


2.5 Answers to check your progress questions 
2.6 Summary 

2.7 Keywords 

2.8 Self-assessment questions and exercises 


2.9 Further readings 


2.0 Introduction 


In physics,a partition function describes the statistical properties of a 
system in thermodynamic equilibrium. Partition functions are functions of 
the thermodynamic state variables, such as the temperature and volume. 
Most of the aggregate thermodynamic variables of the system, such as 
the total energy, free energy, entropy, and pressure, can be expressed in 
terms of the partition function or its derivatives. The partition function is 
dimensionless, it is a pure number. 


Each partition function is constructed to represent a 
particular statistical ensemble (which, in turn, corresponds to a 
particular free energy). The most common statistical ensembles 
have named partition functions. The canonical partition 
function applies to a canonical ensemble, in which the system is 
allowed to exchange heat with the environment at fixed 
temperature, volume, and number of particles. The grand 
canonical partition function applies to a grand canonical 
ensemble, in which the system can exchange both heat and 
particles with the environment, at fixed temperature, volume, 
and chemical potential. A a thermodynamically large system is 
in thermal contact with the environment, with a temperature T, 
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and both the volume of the system and the number of 
constituent particles are fixed. A collection of this kind of 
systems comprises an ensemble called a canonical ensemble. 


2.1 Objectives 


After going through this unit, you will be able to: 
e Understand the concept behind Partition function 
e Learn about how many Microstates are accessible to your system in a 
given ensembles 
e Explains the thermodynamic functions in terms of partition function 
e The physical basis of equilibrium can be understood by using the 
principles of statistical thermodynamics. 


2.2 + Partition Function 


Quantum theory permits a qualitative explanation of energy, heat capacity 
and other related thermodynamic quantities for di and polyatomic molecules. 
A quantitative interpretation is provided by the use of certain general 
function is called "partition function". Mathematically it was written as 


q= Yee P. 
i=0 


Where ‘q’ is the molecular partition function 


Partition function is derived from Maxwell Boltzmann distribution law. 
According to Maxwell Boltzmann distribution law 


n 0 
Hg E! 
n= "ge es 
i i 
S0 


co 


Total number of molecules 2 N — n 
i=0 


N= no +n; +n2 + n3+ ------- ni + --- 
ny 7" n a n n n =a 

N= °ge Ir °ge BH ge TH Lice + “ge ATHY... 
0 I7 T ie n ee 
Eo Eo So So 


13 


Unit -2 
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No - Eg =f - Ej = 
N= [ge | rcge Forge MR cn ge Ki] 
0 1 2 


Eo 
Unit - 2 
Partition Function 
Hs — -E 
N2— ge ee 
NOTE 80:8 


2 -H 
Where >, ge Aris called “partition function”. The summation is 
i=0 
taken over all integral values of ’1’ from zero to infinity. It is 
represented by the letter q or z 


oo — Ei 
q=} ge ir 
i-0 
_ No 
= = By 
Eo No No 


N N when go=1, gt 


i.e partition function is the ratio of the total number of molecules in 
the ground state. Again 


1 1 
4 = 8) —. when gol, q = 


No N nyN 


Partition function is the reciprocal of the mole fraction of the 
molecules in the ground state. It is a mere number and is a 
dimensionless quantity. 


q= Y mE 
i-0 


Partition function will never be equal to zero. It may vary from 
unity to infinity. As T0 q—1(Eo-0). 


Self-Instructional Materiql 
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Unit -2 
Partition Function 


2.2.1 Vibrational Partition Function 
Partition function is defined mathematically as 


q= » ge Ie 


i=0 


NOTE 


Vibrational partition function is given by 


q,— Ye er 
v=0 


Where gy is the degeneracy factor associated with the vibrational level. It is 
unity and Evi» is the energy of the vibrational state in excess of the zero point 
energy. 


The lowest or ground state energy or zero point energy of simple 


v 


harmonic oscillators is given by £ =0 = 2 where ‘v’ is the fundamental 


frequency of oscillation. 


v= Ku Where K is the force constant and -- the 
A 
reduced mass of the system. 


The actual value of energy of vibration is given by the expression 


1 
E. E )hv 


Where, V is the vibrational quantum number. It can take up value from 0 to 

oo and h is the Planck's constant for v0, E 20-7 hv which is the zero 
i 2 

point energy of the vibrator. Hence energy of the vibrational states in excess 

of the zero point energy is 


E=(V+ lnv- l av= Vhv (gv=1) 
i 2 2 
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co  —yhv 


q,= 3 Ze 'kr— 2e KT 
v=0 = 


Unit - 2 


Partition Funcion qa = e” KT + e^" Kr + e?" KT/ E e?" KT 4 


NOTE 


put hv /KT=x= 1+ e" e? x e? 


1 1 g KT 


]-g*. j2-g "*" pi KT) 4 


—2x 


The series 1+ e™ +e" ce? + converges to 


].. 1l a * a € 
]-e* i_! e*—] e-l 
e 
1o 1] 
Qvi — 7 KT. — 


OVQUMKETN 0 0] 
(1 e ) a —gedker) 


v= wave number of oscillation when zero point energy is also 
included. 


= —hv/2KT = hy PKT 
a en =e 


1— (eg F^) pi KT 1 


This expression is used to determine the vibrational partition 
function of a diatomic molecule at all temperatures provided the 
vibrational frequency is known. Vibrational frequency is obtained 
from the study of the spectrum of the molecule. 


If O, = hv/K , the vibrational characteristic temperature, 


Qyiy — 1/1 -e eT 


L Greater the value of 0, (8, =w /K ) i.e. higher the value 


Self-Instructional Material of the vibrational frequency lower will be the vibrational 
partition function and lesser will be the fraction of the 
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molecules in the excited vibrational states. 
2 Greater the value of greater T will be the value of vibrational partition 
function. 


2.2.2 Rotational Partition Function 


Partition function is given by the expression 


q= Yee 
i-0 


For rotational part of the partition function 


oo 


a=, 86 


where gyis the degeneracy factor or statistical weight factor for the rotational 
levels and is equal to (2J +1) where J is the rotational quantum number. Ej is 
the energy of the rotational level in excess of the zero point value. 


Rotational energy of the rigid diatomic molecule is obtained by 
solving the schrodinger wave equation. 


JO -D9g 


E,- ; Where J is the rotational quantum number. J can taken up 


8x71 
values from 0 to o» when J = 0 E;20 


‘h’ is the Planck’s, constant and I the moment of inertia I -yur? where p is the 
reduced mass of the system and r, internuclear distance . Greater the value of 
‘T lesser will be the value of Ej 


Inserting the value of Eyin the expression for partition function, we get 


z — -J e 
a= > Uae DEF 


Provided moment of inertia is moderately large and the temperature not too 
low for all diatomic molecules except H2 and D» the summation is replaced 
by integration and hence 
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i eJ sq 


-5»QJ-lDe———— 
= PACT De KT 
Unit — 2 
Partition Function 
2 
NOTE Assume nE OQ. where @is called rotational characteristic 


temperature 


o xa ZOID Fe. 
t = = r. 
" SZ IKT T 


T 
X.—=J(J+1) 
0 


r 


2 di = QJ +1)dJ 


à. 
Hence, 
-T 
q,= rs E "dx 
T 
=— f e*dx 
q, 6, | 
q= PE 
r 0, 
T 8m IKT 
q, = h = n 
8z^IK 


For heteronuclear diatomics 
8z^IKT 


K 


q, = 
For symmetric molecules (like homonuclear diatomic molecule) 


a factor called symmetry factor represented by the letter o is 
Self-Instructional Material introduced in the denominator. Hence in general for all diatomic 
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molecules. 


8m IKT Unit -2 
Partition Function 


q,— 
on? 


NOTE 


The symmetry factor o is the number of a identical indistinguishable 
configurations a molecule can assume during one complete rotation. 


For homonuclear diatomic molecule like H202, N2, Cl» etc o =2, for 


heteronuclear diatomic molecules like HCl, ICI, CO, NO etc. 


o =1. For linear molecules like CO2 6 = 2 


Interpretation 
_ 82° IKT 
r on 
1. qris independent of the volume of the container 
d: qr depends on the symmetry of the molecule, greater the symmetry 
of the molecule lesser will be the rotational partition function. 
3. qr depends on the moment of inertia and temperature. It is directly 


proportional to I and T. 


2.2.3 Elecronic Partition function 


Partition function may be defined mathematically as 


q= > ge n 
i-0 


Where, gi is the degeneracy factor for the i^ level having energy Ei 
in excess of the zero point level at temperature T. 
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The electronic partition is given by 


oo 


(Ee) 
qe =F (8.)e Yar 


i=0 


(Ee) Ro» 
Erde je Ln T 


—Ee)o 
q= (8. Joe er (g,).e 


q.— (2, Jo e (g, )o (g, Jo 


The energy spacings of the electronic energy levels are very large 
compared to those of vibrational and rotational levels. For most 
atoms except chlorine and molecules except O2, NO and NO» the 
energy of the next higher electronic state is much greater than £e, 


€e,— €&.. 
sothat( ?) is very large at moderate temperatures. 
KT 


Hence second and other terms can be neglected and hence the 
electronic partition function 


(Ee d 


qe= (8h e 


The degeneracy or statistical weight factor for the electronic level 
normal or excited is (2J +1) 


- À Ei 
In general g = e. )e 


—CE,)o 
In the ground state (Ee)ois zero. Hence " xr =1. The 


contribution of this state to the electronic partition function is thus 
(2J +1). 


He, Ne, Na vapour, Hg etc have single electronic ground state 
and for the atoms of these elements the energy difference between 
the lowest and the next electronic level is very high. 


Many monoatomic substances e.g. Cl and a few 
polyatomic molecules e.g.O», NO have multiplet electronic 
ground states and there may be low lying excited electronic states. 
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Hence one or more electronic states above the ground state are 
appreciably occupied even at moderate temperatures and hence Unit -2 
appropriate terms must be included in the partition function. Partition Function 


NOTE 
Example: the lowest state of chlorine atom having the energy (Ee)o = 


zero , the value of j is 3/2 and there is another state energy (Ee); = Ei with J 
value !^ . Hence the electronic partition function qe at ordinary temperature. 


2 -( Eyy 
— KT 
qe 327 
i-o 


(Ee of 


Fe = (8,)9€ T4 (g,)e 


HUE 
KT 


)p 07 *D-Qx +1)=4 


N| w 


P? 
J=3/2 


(g,), "QJ *D2Q x aja 
e71 2 


J= 1/2 
(E. )= 0, (E) E, 


q, = Hada 049 


where, Ei is the energy of the upper level in excess of the zero point level. 
At higher temperatures other terms would have to be included. 
The value of E; = hv is found from the spectrum of chlorine. Thus 


CAS =4+ 2g KT 


ES 4 4 2g KT 


Internal energy with respect to electronic motion 


à lng 
= 2 e 
(Ee), = NKT * ( 5 ) 
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P) _ = 
NET. — ]n(4 42e "* Kr) 
oT 


Unit - 2 uos 1 
E : cV 
Partition Function =R.( ka 
K Qe ry 
NOTE 


= OE, ja 
= CAR ^v 


d = 
= (Rh) 
oT 
9 cv KT |1 


n. e^ kr 


= 2R.(—)* ee, 
KT . Qe" KT y. 1 


By knowing vat a moderate temperature the value of qe, Ee & 
(Cv)e can be evaluated. The possibility of electronic states would 
fail if applied to atomic chlorine at temperatures, greater than 
250K. The higher the temperature the greater is discrepancy. 


2.3 Thermodynamic functions in terms of partition function 
2.3.1 Internal energy from partition function 

Let E be the total energy of N molecules present in the system. 
Average energy <E> is given by 


Y En, 


<E>= 5 aid a) 


(E= >) En, N 2 Yin, from postulate (3) and (4) of Maxwell 
i-0 


1=0 
Boltzmann Distribution Law) 


Applying Maxwell Boltzmann Distribution Law 


Self-Instructional Materiql 
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i.c TEF gae in equation (1) 
8o 
= n " 
yee 
i=0 Eo 


Molar Partition function (Q) is given by 


Q-Y ge ^ 
i-0 


i.e substituting Q value in equation (2) 
< - E;KT 
Ege” 
Hence < E eee (3) 


Q = Tz 
i=0 


Differentiating ‘Q’ with respect to temperature at constant volume 


oQ T ge BYR 


E 
(y= aa 
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Incorporation of this result in equation (3) 


| kr^ 22 
Unit- 2 <E>= oT 1,00, _,dInQ 
Partition Function 0 . QƏT y =( oT y 
NOTE 
< E >= kri 29, 
oT yY 
Internal energy E = KT? ( "ing ) 
oT Y 


q=molecular Partition function qv =Q 


29lnqg” 
E=KT ( 3 
oT 
E= NKT? (234, — (5) 
oT Y 


Hence it is possible to determine the internal energy of the system 
if the partition function is known. 


2.3.2 Heat capacity from partition function 


Self-Instructional Material 
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2.3.3 Entropy from partition function 


Thermodynamic probability is the numbers of ways in which system 


consisting of N identical indistinguishable molecules can be realized. 


W - f 
i-0 n! 
nW =m z 5) 


-In( £z g")—h( Zn!) 
i=0 i i=0 : 


oo oo 


hW = nan 8; -»Mn!3 (1) 


i-0 i-0 
Applying stirlings approximation formula in equation (1) 


In Ni NinN-N 


Iin W =J n,n £i -> (nin n, =H) 


i=0 i=0 


hW = Sin, In gi -Š n,n n, —n, +n 
E = = 
According to Maxwell Boltzmann distribution law 


n= 0g e ex but o N ‘q’-molecular partition function 
So 8o d 


InW= Jnn g —>,n In( B yg eB 37 


i-0 i-0 q i-0 


InW- Sn In g Fann” -F nhe" + 3n=0 
= ] , E, i q ' i i i 


i=0 
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a ii N x - E 
InW= ndn g; -$ nh - nln g, -Š nihe EKT ra n, =0 
i=0 iO d ix = 


co E. oo 
=5 n mE) + +N (Yn =N) 
fo N KT ' 


i-0 


(In and e gets cancelled, by postulate (3) of Maxwell Boltzmann 
distribution law ) 


E 
nW 2 57n In 2) 4 — 4 N (ME = E) by postulate (4) of 


L 


i=0 N KT 
Maxwell Boltzmann distribution law 


E 
InW-NInq-NInnN- _ +N 
KT 


E 
-NInq-(NInN)-N& 
KT 


In Weng dtt. 
KT 


gau 
In W- In Fu 
( ! KT 


Multiplying throughout by K 


N 
Kin W-S- Ki D+ E 


(7S = K In W) (Boltzmann-Plank entropy probability law) 


We know that internal energy 


- Q for identical indistinguishable particles 
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i.e. equation (2) becomes 


one 
S=KIn 
Q T 
S - KlnQ« xr? 9, 
molar oT V 


Molar entropy of an ideal gas can be determined from the knowledge of 
partition function. 


2.3.4 Pressure from partition function 


dA = —PdV — SdT 


--(95 .. %) (22) > (1) 
ov Y 00" ƏV 


Since Helmholtz free energy 
1 

A--KThQ (inQ= 2 

Q 

-KT 
dQ Q 

Substituting (2) in equation (1), we get 


_AT oQ (or) L 1 ,29Q _ 
P=- av ^ ger T 
B kr 29, 

ov T 


2 (2) 


oe) 


2.3.5 Third law and partion function 


Internal energy E of a system is a function of T and V. 
E=f(T,V) 
Applying rules of partial differentiation, 


E E 
dE = (Ë) ar + (2) dV > (1) 
oT Y ov 
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from the combined 1* and 2" law of thermodynamics. 


dE=q-W 


Unit - 2 
Partition Function dE = TdS — PdV 


NOTE 
TdS = dE + PdV > (2) 


Substituting the value of dE from equation (1) in equation (2), we 


get 
dE 9E 

TdS =(~~) dT (^) dV+ PdV > (3) 
oT" ov’ 

(E E 


At constant volume dV=0, i.e. equation (3) becomes 


TdS = CydT 

dT 
dS =C, rJ 
Entropy change for a finite process 

dT 
fas - [6 — 
T] OE ,,0IlnQ 
S,-S,"[2CO, B= KT? (Gp a 
diii J T ƏT” 
Tið dln Q 


S,— seha Cup 


S,- s= forr ero fuav= uv - [vau 
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d 


no, c 
T v pot 


T d 


lyepuolhQ, * 
--KT Car [KT 
ano + 

3T e f KƏ In Q 


cas - kr (2 M2) +[K InQ]" 
0 


T 0 oT V 


= KT ( 


s-s-kr( 2) Kmo -Kno 


T 0 oT V T-T 


-0 


Comparing left and right hand sides and equating the temperature dependent 
and temperature independent terms we get 


s-xrOP?, guo 


T oT V -T 
Sy K In Qro 
When T=0, Q=1 hence KInQ=0 


This is 3 law i.e. All perfect crystalline solids at 0°K will have zero entropy 


2.3.6 Helmholtz Free Energy from Partition Fuction 


Helmholtz free energy is A = -KTInQ 
According to thermodynamic eqn 


A-E-TS 


Substituting the expressions for E and S in terms of partition function 


dln 
E-KT*( cing, (from internal energy) 
oT " 
dln 
S-KlnQo- kr 9j (from entropy) 
oT Y 
A=KTX 99) rücmnos kr(? P9, 
oT Y oT " 
A=KT (92, _ epo kr: (902) 
oT " oT " 
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A ——KT In Q 
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2.3.7 Enthalpy from partion function 


NOTE 
H=E+PV 


;,9InQ 
ELM 


E-KT/( 7) (From internal energy) 


V 


dlInQ 
P= M (From pressure) 


H kr2 (90, , gr (90, y 
oT "V" ov T 
9 In 9 In 

KT{T.( 2) V( 2 
oT "V" ov T 

H = Kr?” 9 , amg, 
ƏT ovW 

TN Kroes dnd, 
olInT OlV 


2.3.8 Gibbs free energy from partition function 


G=H-TS 
=E+PV-TS 
G=A+PV 


A=-KTInQ 


P= gr (29, 
oT T 


G--KTnQ-«kr(9 9) y 


oT T 


ding 
G = KTl( ) -h 
Əln V7 2l 
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2.6 Check Your Progress 


1. Define partion function. 

2. Define canonical ensemples 

3. Give an expersssion for Helmholtz free energy from partition function. 
4. Give the o value (symmetry factor) for H202, HCl, CO». 


2.7 Answers to check your progress questions 


L. 


partition function describes the statistical properties of a system in 
thermodynamic equilibrium. Partition functions are functions of the 
thermodynamic state variables, such as the temperature and volume. 


Most of the aggregate thermodynamic variables of the system, such as 
the total energy, free energy, entropy, and pressure, can be expressed in 
terms of the partition function or its derivatives. The partition function is 
dimensionless, it is a pure number. 


A quantitative interpretation is provided by the use of certain general 
function is called “partition function”. Mathematically it was written as 


2. 


3. 
4. 


q= Yee a 
i=0 


Where ‘q’ is the molecular partition function 


A a thermodynamically large system is in thermal contact with the 
environment, with a temperature T, and both the volume of the system 
and the number of constituent particles are fixed. A collection of this 
kind of systems comprises an ensemble called a canonical ensemble. 

A — —KT nQ 

c values for H202, HCl, CO»- 2, 1,2. 


2.8 Summary 


The molecular partition function indicates the number of thermally 
accessiblestates of a collection of molecules at a temperature T. 
The translational partition function is calculated by noting that 
translational states form a near continuum. When the energy is a sum of 
contributions from independent modes of motion, the partition function 
is a product of partition functions for each mode of motion. 
The internal energy is proportional to the derivative of the partition 
function with respect to temperature. 
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° The parameter p =1/kT. 

° The statistical entropy is defined by the Boltzmann 
formula but may be expressed in terms of the molecular partition 
function. 

. A canonical ensemble is an imaginary collection of 
replications of the actual system with a common temperature. It is 
used to extend statistical thermodynamics to include interacting 
molecules. 

° The thermodynamic limit is reached when the number of 
replications becomes infinite. 

e Most members of the ensemble have an energy very close 
to the mean value. 

° The internal energy of a system composed of interacting 
molecules is proportional to the derivative of the canonical 
partition function with respect to temperature. 

. The entropy of an interacting system can be calculated 
from the canonical partition function. 


2.9 Keywords 


L. The molecular partition function indicates the number of 
thermally accessiblestates of a collection of molecules at a 
temperature T. 
2. A canonical ensemble is an imaginary collection of 
replications of the actual system with a common temperature. It is 
used to extend statistical thermodynamics to include interacting 
molecules. 
3. The thermodynamic information in partition function 
a) The internal energy of a system composed of interacting 
molecules is proportional to the derivative of the canonical 
partition function with respect to temperature. 
b) The entropy of an interacting system can be calculated 
from the canonical partition function. 
The following functions are written in terms of the canonical 
partition function: 
(a) Helmoltz energy, 
(b) Pressure, 
(c) Enthalpy, 

(d) Gibbs energy. 
4. The molecular partition function factorizes into a product of: 
a) translational, 


b) rotational, 
c) vibrational, and 
d) electronic contributions. 


The contributions to the overall partition function are summarized 
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2.10 Self-assessment questions and exercises 


1. Define Molecular partition function and ensembles. 
2. Define partition function and thermodynamic functions 
interms of partition function. 


2.11 Further readings 


1. Statistical Thermodynamics, M.C. Gupta, Wiley 
Eastern, New Delhi,1990. 


2. Introduction to Statistical Thermodynamics, 


R.P.H.Gasser and W.G.Richards, World Scientific, Singapore, 
1995. 


3. Physical chemistry, Peter Atkins, Julio De Paula, 9™ 
edition, 2010 


33 


Unit -2 
Partition Function 


NOTE 


Self-Instructional Material 


Unit - 3 
Statistical Interpretation 
of Third Law 


NOTE 


UNIT- III statistical interpretation of third law 


Structure 
3.0 Introduction 
3.1 Objectives 
3.2 Bose Einstein Distribution Law 
3.21 Photon Gas — Application of Bose Einstein Law 
3.2.1.1Derivative of Planck's Black Body radiation law 
33 Fermi — Dirac Distribution Law 
3.3.1] Electron Gas — Application of Fermi Dirac Distribution Law 
3.3.1.1 Determination of Average Energy of Free Electrons in Metals 
3.4 Check Your Progress 
3.5 Answers to check your progress questions 
3.6 Summary 
3.7 Keywords 
3.8 Self-assessment questions and exercises 
3.9 Further readings 


3.0 Introduction 


The Third law of thermodynamics is sometimes stated as follows, 
regarding the properties of closed systems in thermodynamic 
equilibrium: 


The entropy of a system approaches a constant value as its temperature 
approaches absolute zero. 


This constant value cannot depend on any other parameters 
characterizing the closed system, such as pressure or applied magnetic 
field. At absolute zero (zero kelvin) the system must be in a state with 
the minimum possible energy. Entropy is related to the number of 
accessible microstates, and there is typically one unique state (called 
the ground state) with minimum energy. ' In such a case, the entropy at 
absolute zero will be exactly zero. If the system does not have a well- 
defined order (if its order is glassy, for example), then there may remain 
some finite entropy as the system is brought to very low temperatures, 
either because the system becomes locked into a configuration with non- 
minimal energy or because the minimum energy state is non- unique. 
The constant value is called the residual entropy of the system."! The 
entropy is essentially a state-function meaning the inherent value of 
different atoms, molecules, and other configurations of particles 
including subatomic or atomic material is defined by entropy, which can 
be discovered near 0 K. The Nernst-Simon statement of the third law of 
thermodynamics concerns thermodynamic processes at a fixed, low 
temperature: 


The entropy change associated with any condensed system 
undergoing a reversible isothermal process approaches zero as 
the temperature at which it is performed approaches 0 K. 


Here a condensed system refers to liquids and solids. A classical 
formulation by Nernst (actually a consequence of the Third Law) is: 

It is impossible for any process, no matter how idealized, to reduce the 
entropy of a system to its absolute-zero value in a finite number of 
operations. 


There also exists a formulation of the Third Law which approaches the 
subject by postulating a specific energy behavior: 


If the composite of two thermodynamic systems constitutes an isolated 
system, then any energy exchange in any form between those two systems 
is bounded 


3.1 Objectives 


After going through this unit, you will be able to: 

Explain about statistical interpretation of third law 

Understand the concept of Bose Einstein distribution law and it's 
applications 

Understand the concept of Ferrmi Dirac distribution law and it's 
application. 


3.2 Bose Einstein distribution law 


The statistics applicable to particles whose total wave function is 
completely symmetric is known as Bose — Einstein distribution law. Bose 
Einstein statistics is applicable to Bosons, which are particles with integral 
spins. i.e 0,1,2,3,4--- ex:photons, ;*He, ?H, HN, 160, CO» etc. In the 
statistics thermodynamic probability is the number of ways in which the 
(ni) indistinguishable particles can be placed in distinguishable boxes or 
energy levels (gi) without limiting the number of particles placed in each 
box or energy level. 


“ (g; I] d M 
w- 4S n) 
HF Us -Dingj 


Where W is the total probability, niis the number of particles in the i^ level 
having the degeneracy gi. 


Consider a system of N particles of an ideal gas, which are 
identical, independent and indistinguishable let no particles have energy Eo. 


ni particles are in the energy state E1. 
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n2 particles are in the energy state E2 


ns particles are in the energy state Ei. 


Such that ` nj=N and 


i=0 


> En =E 
i=0 
Total number of particles in the system is constant and hence there is no 


variation in the total number of particles. 


ON = Oie 


Š Eðn,=0 > (1) 
i-0 


Total energy of the system is constant 


QE -0i.e 


Y En 0 (2) 
i-0 


Equations 1 and 2 are the two conditions of constraint for a system of gas 
in statistical equilibrium. 


Maximization of thermodynamic probability 


According to Bose — Einstein statistics thermodynamic probability is the 
number of ways in which the system of (ni) identical indistinguishable 
particles can be placed in different energy levels (gi) without any 
limitation to the occupation number is 


_ z (gi—1+n,)! 
i0 (g;—1)!n;! 


Where niis the number of particles in the i^ level having 
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the degeneracy gi and gi>>ni 


Since giis very large (gi >>1) one can be neglected and Unit —3 
hence Statistical Interpretation 
of Third Law 
W = i i 
i=0 gi tn ! 
i-0 $; In. ! 


In W= y» g,* n,)!-» ]n( g,)!->_In( nl 
i=0 i=0 i-o 


On using stirlings approximation for factorials of large 
numbers. 


nW = Y en) htg,en)- (8: * n)]— Dime, -8)- Lohn cn) 


-0 


W- Ys noti +n)- va -Zn-Zehs,+ Ys DX 


i=0 i=0 iH 


- (g; +n) g; n) - 9g, n 8,2, nhn, 


i=0 i=0 i=0 


For a system which has attained statistical equilibrium 
there will be most probable distribution of molecules and 
the thermodynamic probability of the system will be 
maximum and hence, 


OW =0 (or) dln W=0 


olhW= Derm j 2 em m 0- Xo tha yn =0 
n i 


i=0 i i i=0 i 
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Unit — 3 i-0 i-0 
Statistical gn 
Interpretation of aln W= > "inc i ‘dn =0 
Third Law as n; i 
| NOTE 9InW- Sint +8")an =0 

i=0 n; ' 


t 


9InW =-FIn(1+ on =0 > (3) 


i=0 n; 


Li 


Equation (3) is the third condition of constraint. 


Since there is more than one constraint in the system Larangian 
method of undetermined multiplier is applied to solve the problem. 


Equation (1) multiplied by a, (2) by p and combined with equation (3), 
we get 


a Yon - BY, En- F nat 9n 20 
io i0 =O n, | 
Where a &f are the lagrangian undetermined multipliers. 


Ya BE - In(1- iðn) =0 


i=0 n; 


L 


dn, s will vary independently and hence 0n,# 0 and coefficient of Om, 
must be equal to zero. 


a+ BE - (1 5)-20 


n. 


L 


nd+)= a+ pE 


$i z (e%PE)—] 


= Si = 
Self-Instructional Materiql nm gf) = a=— |RT B=1kT 


n= e WRT ET 


As T—=œ ~=e**F* >>] and hence ‘1’ in the denominator can be neglected to 
get Maxwell Boltzmann distribution law. At low temperatures and for 
problems concerned with radiation ‘l? in the denominator cannot be 
neglected. 


Applications 


1. Bose Einstein distribution law is useful to determine the number of 
identical and indistinguishable particle in the i^ level having the 
degeneracy gi with energy Ei excess of the zero of the zero point level. 


2 It is used to explain the behavior of Helium at low temperature 
through Bose — Einstein condensation. 
3. It is used to explain the radiation by considering them as photons. 


It is used to derive Planck’s Black body radiation law and all the classical 
laws of black body radiation. 


3.2.1 Photon Gas - Application of Bose - Einstein Distribution Law 


3.2.1.1Derivation of Planck’s Black Body radiation law 

According to Bose-Einstein Distribution law the number of identical 
indistinguishable particles ni having the specified energy Ei in excess of 
Zero point energy is 


8i 
n;— —Bádu——— 2 U 
és B = ( ) 


Where, ‘g; is the degeneracy factor of the translational levels. 


Electromagnetic radiations consist of discrete energy particles 
called photons contained in a container of volume v with a definite 
energy. Photons do not interact with each other and hence a very small 
black body is assumed to the present in the container to absorb and emit 
photons and thus thermal equilibrium is possible. Gibbs free energy of 
radiation is zero at equilibrium and hence p = 0 and @=— ufRT —-0. 
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For this system Bose Einstein distribution law i.e. equation 
(1) becomes. 


"E 8i 


i BE; 
e -l 


> (2) 


(a= - p/RT=0) 


dgi, the degeneracy of the translational states is equal to the 
number of lattice points contained in 1/8" of the spherical 
shell of radius n and thickness dn . the number of lattice 
points corresponding to the specified energy range us given 
by 1/8" the volume of the spherical shell of radius n and 
thickness dn. 


| Amdn  mdn 
8 2 


dg 2 (3) 


The length ‘Il’ of the box is related to the wavelength of the de-Broglie 
wave in particle of one dimensional box as 


I2 n 422 y? 
p ux siey 
we a v2 
c 
3n'dn = SY ym nn = STAY 


Substituting the value of n’dn in equation (3) 

_mvy'dy >() 

z 88 4 
2 c 


dg 
Radiations are polarized in two ways (i.e. equation (4) x2) 


Hence, dg — pc 8vydy 
2 e 
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87v 
dg- .y’dy> (5) 
E 


Replacing gi, by dgi and ni by dni and substituting dg value (5) 
in equation (2) we get 


ne dud 
| 68 (e)-1 


‘dN’ is the number of photons in the energy range from v and v+dv. To 
obtain the total energy in the range v and v+dv, dN is multiplied by hv, 
which is the energy of one photon. 


dE — hy.dn 


8w dy 


— -E 
(ema 


dE =hy. 
"S 


Energy density=energy per unit volume. 


dn P _8my dy 
y c^ (ele -]1) 


dE | 8zhyd 
Lm mY S, L2EQy 


Bd KT 
v cle” — 


This is Planck’s radiation law in terms of frequency conversion in terms of 


À. 

4 c 
y=c/A ye 

A 
8zhc'dÀ b ee 
du- CA ETTI 4ydy- c .—4A dA 
8mc ^ dA c 

du- 4. (7-1) Lgs " 
du=E,dr 
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This expression is Planck's radiation law expressed in terms 


of wavelength. 


Unit- 3 
Statistical. When energy density is plotted as a function of A, the 
e following graph is obtained 
Third Law & grap l 
| NOTE - 
Emax | / i 
(2000K) | fs \ 
/ | 
du ^ 
d | / 4 A 
; Emax | ! AX 
(ene " 
puis d (1500K) y i Fi V - 
| Ww N EC 
Emax | / A oat N Ae 
(1000K) C a DR | 
sd | "t T - 
"<< RR OR Ew 
Amax Amax Amax 
(2000K) (1500K) (1000K) 
A a’. 


Experimental observation 


1. Energy is not distributed uniformly throughout the spectrum. 


2. Energy density is minimum both as A-0 & X—»oo. It increases as À 
increases, attains a maximum then decreases. The wavelength at which 
energy density is maximum is called A max and the maximum energy 


density corresponding to Amax is called Emax. 


Emax increases with increase in temperature and Amax decreases with 
increase in temperature. 


3. Total energy density increases with increase in temperature. 


All the experimental observation is explained by Planck's radiation 


law. 


33 Fermi - Dirac Distribution Law 
Need for Fermi - Dirac distribution law 
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l. Maxwell Boltzmann distribution law fails to explain the 
low temperature behavior of helium and the spectrum of black — 


body radiation. 

2. It does not explain the properties of "electron gas" in 
metals. 

3: It fails to explain the behavior of ideal gas containing 


identical indistinguishable molecules. 


Derivation of Fermi-dirac distribution law 
Consider a system of N identical independent indistinguishable particles 
of an ideal gas. Let there be no particles in the energy state Eo, ni particles in 
the energy state E;, n» particles in the energy state E2 and so on. 
1 Total number of particles > n;=N 
i=0 


Since the total number of particles in the system is constant. 


> an, 20 5 (1) 
1=0 


This is the first condition of constraint. 


2 Total energy of the system is constant 


Y Eðn,=0 (2) 
i-0 


This is the second condition of constraint. 


Maximization of thermodynamic probability 


According to Fermi- dirac statistics the thermodynamic probability or 
the number of ways in which ni identical indistinguishable particles can be 
distributed among gi energy levels of energy Ei. Such that not more than one 
particle can occupy a given level is given by 


The total number of ways for the whole system consisting of N particles 
would be 
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i ii gi! 
TW,=W=2 
= i (g,— n,)!n;! 


nW=Inz gi! — (a) 


i-0(g,— n,)!n;! 
Applying stirlings approximation formula in equation (a) 


In W-In zg;!- In z(g;j- n;)!- In zn! 
i-0 i-0 i-0 


InW= Yn g,! -J In (g,— n)!-Y In n,! 
i=0 i=0 i=0 
In N!=NInN-N 


In W 2 * 8,Ing,->)(g,-n,)In(g,-n,)— >on, In n; 
i-o i-0 i-0 


For a system in statistical equilibrium thermodynamic probability 
must be maximum. Hence maximization of thermodynamic 
probability must be maximum. Hence maximization of 
thermodynamic probability leads to ð In W = 0 


dInW=0- Sgn. an In(g. -n)ón)- Yn dn &In nón)-0 


i 
i-0 St; i=0 n, 


-—M (1+ In( g,- n,)on;- (1+ In n,)ðn,= 0 
i-0 i-0 


= -Yn 5r Uu )on =0 


i20 n; 


t 


--Y In(*'-1an =0 > (3) 
n. i 


i-0 i 


This is the third condition of constraint. In order to solve the 
problem for the system consisting of two or more conditions of 
constraint Lagrangian method of undetermined multipliers used. 
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Equation (1) is multiplied by a and equation (2) by p and 
combined with the condition (3), a and f are the lagrangian 


undetermined multipliers. Unit - 3 
Statistical Interpretation 
of Third Law 
aS on * BY, En- Sin *'—1)dn 20 NOTE 
i i i i= n; ! 


i=0 i=0 
Ya BE -n( 8'-1)dn = 0 


i=0 n; 


t 


Since àn;'s can vary independently ðn, #0 and hence 


a+ BE - In( 5-1) =0 


n. 


1 


In( 5-1) = @+ BE 


un 
Bi | = Qr BED 
n; 
Si | e+ BE) 44 
n; 

Si 
n = 
en 
e(%* BE) 41 

n = 


gi 
i eC RT+ BIET) 4 
Equation (4) is the mathematical expression of Fermi-Dirac distribution law. 


Application 


1. It is applied to explain the behaviour of "electron gas" in metals and 
thermionic emission and semiconductors. 

2. It helps to determine the number of particles in the i" state with 
degeneracy gi, having energy Eiin excess of zero point energy. 


3.3.1 Electron Gas — Application of Fermi-Dirac Distribution Law 


3.3.1.1 Determination of Average Energy of Free Electron Metals 
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A metal, each of whose atoms contribute one or more free 
conducting electrons is regarded as "electron gas". 


The number of electrons having energy in the range between 
E and E«dE is given by Fermi-Dirac distribution law. 


The number of energy levels in the range between E and E+dE 
for electron in a box of volume is twice that for particles in a 
box of volume is twice that for particles in a box of volume V 
because electrons have spin of +1/2 or -1/2. So for electrons, 


dg~ 2” “4, dn 1 

nn? 2 23 
E (aov a =v ) 
2 8nv PE 


1217 V3 p12 
n m) VPE 


h 
3 (8m)}’ VE]? 
ie 
32 
3n2dn = Sr .V.3 2 E! dE 


Substituting the value of n’dn in equation (1), we get 
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8m)! ?.z.V .E!?dE 
an - Gm? m. 
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dN _ (8m) zm V.E" of Third Law 
LEER NOTE 
2h (e * +1) 


This equation is called Fermi — Dirac formula of free electrons. This 
represents energy distribution for free electrons. 


At absolute zero (i.e. T=O). The particles occupy the lowest level upto 
Er, Fermi energy. Fermi energy indicates the maximum energy of fermions in 
the system. 


If there are N electrons in the metal at absolute zero and if the 
maximum energy of an e^ is EFo is for O<Ej «EF; 


EP, 3/2 
8m)” .zXV 
N=| a= e ap 


0 on(e +1) 


At absolute zero when E; <EF, 


"^ (8m) 2. V 


Jav-v-|— p 
3 


(8m) m.V . EF 3? 
N= 
2h 3/2 
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N/V is the number of free electrons per unit volume. It is the 
free electron density (p). Fermi energy for free electrons in 
metals can be calculated by knowing the free electron density. 


3.4 Check Your Progress 

1. What is the third law of thermodynamics? 

2. What is the need of Fermi Dirac distribution law? 

3. Give the applications of Bose Einstein law. 

4. Give the applications of Fermi Dirac distribution law 


3.5 Answers to Check Your Progress Questions 
1. Third law of Thermodynamics 
. The third law of thermodynamics is sometimes stated as 


follows, regarding the properties of closed systems 
in thermodynamic equilibrium. 


° Here a condensed system refers to liquids and solids. A 
classical formulation by Nernst (actually a consequence of the 
Third Law) is: 


e It is impossible for any process, no matter how 
idealized, to reduce the entropy of a system to its absolute-zero 
value in a finite number of operations. 


° There also exists a formulation of the Third Law which 
approaches the subject by postulating a specific energy 
behavior: 


e If the composite of two thermodynamic systems 
constitutes an isolated system, then any energy exchange in any 
form between those two systems is bounded. 


2. Need for Fermi Dirac distribution law 


° Maxwell Boltzmann distribution law fails to explain the 
low temperature behavior of helium and the spectrum of black — 
body radiation. 


e It does not explain the properties of “electron gas” in 
metals. 
. It fails to explain the behavior of ideal gas containing 


identical indistinguishable molecules. 
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3. Applications of Bose Einstein distribution law 


. Bose Einstein distribution law is useful to determine the number of 

identical and indistinguishable particle in the i" level having the degeneracy Unit - 3 

gi with energy Ei excess of the zero of the zero point level. Statistical Interpretation 
. It is used to explain the behavior of Helium at low temperature of Third Law 

through Bose — Einstein condensation. NOTE 

° It is used to explain the radiation by considering them as 


photons. It is used to derive Planck's Black body radiation law and all 
the classical laws of black body radiation. 


4. Fermi Dirac distribution law 

° It is applied to explain the behaviour of “electron gas” in metals 
and thermionic emission and semiconductors. 

. It helps to determine the number of particles in the i" state with 


degeneracy gi, having energy Eiin excess of zero point energy. 


3.6 Summary 


e Bose-Einstein statistics describe one of two possible ways in 
which a collection of non-interacting, indistinguishable particles may 
occupy a set of available discrete energy states at thermodynamic 
equilibrium. 

° The Bose-Einstein statistics apply only to those particles not 
limited to single occupancy of the same state—that 1s, particles that do 
not obey the Pauli exclusion principle restrictions. Such particles have 
integer values of spin and are named bosons, after the statistics that 
correctly describe their behaviour. There must also be no significant 
interaction between the particles. 


. Fermi-Dirac and Bose-Einstein statistics apply when quantum 
effects are important and the particles are "indistinguishable". Quantum 
effects appear if the concentration of particles satisfies 


° Fermi-Dirac statistics apply to fermions (particles that obey 
the Pauli exclusion principle), and Bose-Einstein statistics apply 

to bosons.. Both Fermi-Dirac and Bose-Einstein become Maxwell- 
Boltzmann statistics at high temperature or at low concentration. 


3.7 Keywords 


The Bose Einstein distribution describes the statistical behaviour 
integer spin particles (bosons). At low temperature, bosons can behave 
very differently than fermions because an unlimited number of them 
can collect into the same energy state, a phenomenon called 
"condensation". 

The Fermi Dirac distribution function, also called Fermi 
function, provides the probability of occupancy of energy 
levels by fermions. Fermions are half — integer spin 
particles, which obey the Pauli exclusion principle. The 
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Pauli exclusion principle postulatestht only one fermion 
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3.8 Self-Assessment Questions and Exercises 
| NOTE 1. Derive Bose Einstein distribution law and give its 
applications. 


2. Derive Plank's black body radiation law using Bose Einstein 
distribution law. 

3. Derive Fermi Dirac distribution law and give its applications. 

4. Determind the average energy of free electron metals using 
Ferrmi Dirac distribution law. 


3.9 Further Readings 
1. Statistical Thermodynamics, M.C. Gupta, Wiley 
Eastern, New Delhi,1990. 


2. Introduction to Statistical Thermodynamics, 
R.P.H.Gasser and W.G.Richards, World Scientific, Singapore, 
1995. 

3. Physical chemistry, Peter Atkins, Julio De Paula, 9" 


edition, 2010 
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Unit -4: Heat Capacity of solids 


Structure 

4.0 Introduction 

4.1 Objectives 

4.2 Heat Capacity 

4.2.1] Classical observation of heat capacity of solids 
4.2.2 Einstein theory of heat capacity of solids 

4.2.3 Debye theory of heat capacity of solids 

4.3 Non equilibrium thermodynamics 

4.3.1  Postulates of local equilibrium Onsager formulatio 
4.4 Check your progress questions 

4.5 Answers to check your progress questions 

4.6 Summary 

4.7 Keywords 

4.8 Self-assessment questions and exercises 

4.9 Further readings 


4.0 Introduction 


Heat capacity or thermal capacity is a physical property of matter, 
defined as the amount of heat to be supplied to a given mass of material to 
produce a unit change in its temperature. The SI unit of heat capacity is joule 
per Kelvin (J/K). Heat capacity is an extensive property. The corresponding 
intensive property is the specific heat capacity. Deviding the capacity by the 
amount of substance in moles yields its molar heat capacity. The volumetric 
heat capacity measures the heat capacity per volume. 


4.1 Objectives 
e To understands the heat capacity of the solids and classical 
observation of heat capacity of solids 
e To understand the Einstein theory of heat capacity 
e To undertand the Debye theory of heat capacity 
e To understand the Non equilibrium thermodynamics and 
Onsager reciprocal relation. 


4.2 Heat capacity 


Heat capacity is the quantity of heat required to raise the temperature of the 
whole of the substance through 1K. 


Even at absolute zero the atoms in a solid will be vibrating 
about their equilibrium positions. When the temperature is 
increased the amplitude of vibration increases and vibrational 
energy of the atoms of the solid is increased, it is the rate of 
change of vibrational energy of the solid with temperature that 
determine the heat capacity of solids is 
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4.2.1 Classical observation of heat capacity of solids 


An ideal solid is one which consists of space lattice of 
independent atoms  vibrating about their equilibrium 
positions. They do not interact with each other. Each atom is 
considered to be a simple harmonic oscillator and can vibrate 
in three mutually perpendicular direction and can have three 
vibrational degree of freedom. 


According to equipartition principle of energy each 
vibrational degree of freedom contributes one KT to the total 
energy of crystal per atom. A solid consisting of 3N atoms 
will have 3N vibrational degree of freedom and hence will 
contribute 3NKT or 3RT/gram atom. 


E,,7 3RT 
E. 
TN T _ (BRT aR 
" 9T " oT 


3x1.1987cal/degree/gm.atom 
Cy= 5.941 cal/degree/gm.atom 
Cy = 3x8.314 cal/degree/gm.atom 
Cy = 24.942 cal/degree/gm.atom 


According to the classical principle the heat capacity at 
constant volume must be constant and independent of 
temperature. Heat capacity measurements are usually made at 
constant pressure and hence C, can be converted into Cp 
using thermodynamic relationship. 


Where oy is the coefficient of thermal expansion and f 
compressibility factor , V= atomic volume. 
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4.2.2 


Einstein’s theory of heat capacity of solids 


Classical theories failed to explain its behaviour of variation of atomic 
heat capacities of solid elements with temperature due to the assumption 
that solid elements absorb heat continuously. Einstein used Quantum 
theory to explain all the experimentally observed regarding the variation 
of heat capacity with temperature. 


Assumption of Einstein’s theory 


l. 


2. 


o2 


d 


Heat is absorbed or radiated in the form of discrete packet called 

photons. Each photon has energy equal to hy. 

Each solid consists of atoms arranged in space lattice and are 

independent. 

Each atom is considered as a simple harmonic oscillator. 

The atoms vibrate about their equilibrium position with a 

uniform frequency which is independent of the presence of 

neighboring atoms. This frequency is represented as y and is 

called Einstein's characteristic frequency. 

The frequency of vibration of atom is characteristic of the 

particular solid concerned. 

Each atom has got 3 independent vibrational degrees of freedom. 

The average ‘E’ per degree of freedom is not equal to *Kt" given 
v 

by equipartition principle, but is equal to ZA as calculated 

b planck's using quantum theory. 


Vibrational energy associated with one atom 


hv 


C RT] 


Vibrational energy associated with 3N atom 


3Nhv 


e^ KT E 


Heat capacity at constant volume 


9E, 
kac ar" 
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C= er 
y v/KT — 

Unit -4 oT e 
Heat Capacity of Owe KT -1 
soilds = 3NhvE . (e =I) 

oT 

NOTE 


Multiplied and divided by K 
hvEJKT hvE 


3NhvE (e? Kr ly KT 2 


hvE > give KT 
/ 


e E KT hvE 23 ( ) 


“3R. T, 
3 (PET er? 
g 
K 


Na= Avagadro number 
n = number of gram atoms 
N=Na, Nak=R 


Put 6,= hvE/K, where 6, is called Einstein’s characteristic 
temperature. 


gor Ia () 
Qu RU asl »- 2 


(e™T—1)} 
Equation (1) and (2) are called Einstein’s heat capacity 
equation. 


Factor explained by Einstein’s theory 


According to Einstein’s theory heat capacity equation 


Self-Instructional Material 
l. Cv varies with temperature. It is a function of and it is 
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related to temperature in an exponential manner. 

2. It explains the shape of the curve . since Ó, is 

related to veas 0—- ^, Cv is a function of VE /T, VE 
"OK 

is the Einsteins characteristic of the element under 

consideration. The solid having more or less comparable 

VE values will have the same type of curve. When Cv is 


plotted as a function of T/O.. for many elements the 
following graph is obtained . 


3. As T—0 Cv attains zero value 
Einsteins heat capacity equation is 


eT 0 
C=3R— — uy 


As T0, 8, /T—o , e? |T »»1. 


And hence ‘T in the denominator of Einsteins heat capacity of 
equation can be neglected and hence the equation becomes 


0 
et" Q, p iT ; 
Cy= SR ATEM) = SRS € 
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As T—0 
eT 46? =) 
Hence Cy =0 


4. As T—^o Einsteins equation reduces to Dulong and Petits law 
is Cy= 3R Einsteins Equation is 


er 


C= 3R are 


y 


( Scy 


As Too, 6,/T will tend to zero, /T is expanded as power 
series of 0,/T. 


3 2 
ves (8 Um. (8, /D " 
2! 3! 


(oed eet 


(F ? 244 (F/T) (Gc/T Y. 
2! 3! 


As T—»oo, 6,/T will be very small and except to I terms all other 
terms can be neglected. Hence, 


(e% F zy os 0. a 
T 
As T—, substituting this condition in the Einstein heat capacity 
equation 


: 0 
C,-3R €"" xE =3R@ =3R 
(QIT T 


As Too, C, Z 3R . This is Dulong petits law. 
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Limitations 


1. It predicted satisfactorily the value of specific heat of solids at high 
temperature. But at low temperature the predicted value 
are found to be lesser than the experimental values. 


2. It fails to describe the behaviour of crystals near absolute zero. 
Failure of this theory is due to the assumption that all the 
atoms in the crystals are independent and oscillate with 
uniform frequency. He neglected the mutual forces 
exerted by the atom upon each other Debye modified 
Einstein's theory. 


4.2.3 Debye theory of heat capacity of solids 


Einsteins theory fails to describe adequately the behaviour of crystals 
near absolute zero. Failure of Einstein's theory is due to the 
assumption that all the atoms in a solid oscillate with uniform 
frequency. He neglected the mutual forces between the atoms. Due to 
Proximity of atoms in a solid there will be interaction and hence the 
atoms cannot osciallate with uniform frequency. A solid containing N 
atoms will behave as N coupled oscillators and their frequency will 
vary from zero to a maximum value ( v max). In order to determine the 
distribution of frequencies Debye disregarded the structure of the solid 
and he treated it as a homogeneous continuous elastic medium the 
vibrations of atoms are considered as equivalent to elastic waves 
propagated in the elastic medium with different frequencies. 


When a continous solid is subjected to elastic vibrations, two 
kinds of vibrations are produced. 


1. Transverse vibration and 
2. longitudinal vibrations 
According to vibrational theory 


Number of modes of longitudinal vibrations 


Aav?^dv 3 
l 


Per unit volume with frequencies in the range v+dv is= 


Where --- is the velocity of the longitudinal vibrations. 
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umber of modes of transverse vibrations 


Unit - 4 Per unit volume with frequencies in the range v and v+d v is = 
Heat Capacity of Aav^dv 
soilds =2x 3 
C, 
NOTE Where Cris the velocity of transverse vibrations ( since 
transverse vibrations have two independent directions of 
4mdv 
motion E is multiplied by 2). 
t 


Total number of independent vibrations in the frequency range 


vto v+dv = Anus d 
l t 


If V is the volume of one gram mole of the solid . then, 


Vmax 


. : 1 2.2 
Total number of vibrations _ — + poi dv 
=| 4m Cl Ci 


where Vmax = vpis called Debye characteristic frequency. 


If there are N atom in volume V the possible vibrational 
degrees of freedom = 3N 


Hence 
1 3 
3N= 4av(— + — yas 
Cl Ct 3 
3 
1 2 
v Cl Ct 


The degeneracy associated with a range of frequencies from v 
and v+dv was calculated by Debye using vibration theory. It is 
found to be proportional to vdv. 


dg œ vdv 


1 1 
dg = cv°dv wherecz 4z( + ) | 
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1 2 
dg = a a dv Unit — 4 


d 3 Heat Capacity of 
soilds 
9N 
dg ——— vd 
E 
"y NOTE 


Total number of vibrational degrees of freedom in a crystal consisting 
of N atoms will be 


9N 

dg = —— v!dv 
Vin 
9N 

dg = —— v!dv 
Vp 


According to Debye the total energy of a crystal is obtained by 
adding the energies associated with all the various vibrational 
frequencies property weighed by their degeneracies. 


Average energy of an oscillator is XE,» = hy ES hv 


2 gre L] 
Where hv/2 is the zero point energy. 
Total vibrational energy of a crystal 


Y max 


ke 


The first integrand on integration gives the total ZPE of the crystal, 
which is equal to 


Pmax hv hv 9Ny?dv 
pa oo — 00 
2 (ehv/KT-l) Vinx 


phe. 9Nv'd  ""9Nhvdv — 9Nhv^« ... _ 9Nhv 


J 2 te J oye. 2v 4 9 8 
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apap 4 
put, x= Av/KT v = T "m 
x.KT 3 44 
v= 4v dv= E Ax dx 
h hê 


When v = v max and applying limits from 0 to x, we get 


ONh (KT)* MR kt x? 


E= dx 
E e =] 


Let v o- v max be the Debye characteristic temperature defined 
f hv hv, 
mathematically as 0, = —- = ae 


Where Opis the Debye characteristic frequency. 


3 x-1 
Gp $e 
This is integral can be evaluated in two special cases. 


l. At moderate and high temperature above 30K. and 
2. At very low temperature below 30K. 


Debye theory at moderate and high temperature 
9RT^ 97 x3 


Gn Ve 


For value of T/0 the integral can be evaluated by 
following series expansion. For large value of T/0 the 
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upper limit of the integration becomes small. First e* is 
expanded as a power series and divided it into x^ by long 
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m d x NK EE desse teRR SR E 
gU 2 12 760 PM 
9RT* er gk at us 
E- 2. —t——-——— +)dx 
i | (* —2 12 760 


Upon integration 
1 1 
E - 3RT (1-—@ T Y - —(8 [ry — ) 
i 20 / 560 / 
C 


9E 1 1 
ECT =3R- (0T + (8 mee 
" QT" 20 560 


As T—œo the heat capacity approaches 3R as predicted by Einsteins 
theory 


1. Cv is a function of O0p/T only and a plot of Cv against T/Op or log 
T/0p should yield a curve that is the same for all the solids, the curves 
are called Debye curves. 
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2. From the curve it is evident that heat capacity of an element 
attains it's classical 


Unit — 4 
Heat Capacity of value of 3R when T/0p is unity. 
soilds 
3. From Debye equation it is possible to calculate the 
heat capacity of any solid element at all moderate and high 
NOTE temperatures if Debye characteristic temperature is known. 


Debye equation at low temperature below 30K 


9RT* 9/ dx 


Be—À— ep 
a x 
when T/0—0, 0/T—2oo 
ORT * x? dx 
Hence, E= 
6, J ge 
f xX dx B x 
4 e™! 15 
9RT^ z^ 
gl 
63 "1s 


ðE 
Cv- E upon differentiating of energy with respect to 


temperature 
ð O9RT^ , 


Cvz 
———.T 
oT 0; .15 


) 


ORI 
= T 


Du 
0; 15 
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4 T 
SORT (^y 


Cv 
15 0, 
Unit — 4 
T ; 
Cv = 464.6... y cu UNE 
Op 
This equation is called Debye third power law. NOTE 


1. According to this equation the atomic heat capacity of elements 
is proportional to T? or Cv'?, 

2. It is used to determine the heat capacity at very low temperature 
by extrapolation of the curve. This is used to determine the 
absolute entropies of substances using 3™ law of the 
thermodynamics. 

3. Op can be obtained from the slope of the curve obtained by 
plotting Cv vs. T. 


4. When 0p is known, Debye characteristic frequency vp Or Vmax 
can be calculated using the expression 0p = hyp/ K. 

5. The most significant difference between the two theories i.e. 
Einstein and Debye theories is at low temperature. The debye 
expression for Cv approaches zero much less rapidly as TO 
than does the Einstein's predictions. 
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It is in this respect Debye theory is superior to Einstein theory. 


Limitations of Debye theory 


0p should be independent of the source of measurement but it 
is not. 

0p Vs T should be a horizontal straight line since Op is a 
constant and independent of temperature. But there are 
deviations in the case of Ag and Na I. 


In the case of metals there is electronic contribution to the heat 
capacity due to free electrons, which has not been taken into 
account 

Cv= R(T /T;) , I r-fermi temperature. 


Debye theory considers the solid to be a homogeneous 
continous elastic medium without any definite structure. It 
may not be exactly correct. 


4.3 Non equilibrium thermodynamics 


The branch of science dealing with the study of 
thermodynamic properties of the system which are not in 
equilibrium and involves transport process which are 
irreversible is termed as irreversible or non equilibrium 
thermodynamics. 


Onsager reciprocal relation 


4.3.1 Postulates of local equilibrium Onsager 
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formulation 


Phenomenological laws, onsager reciprocal relations. 


Onsager theory is based on the principle of microscopic 
reversibility which states that under equilibrium conditions 
any molecular process and the reverse of this process will be 
taking place on the average at the same rate. According to 
onsager, the rates or velocities of the various types of 
processes like diffusion, flow o heat etc are linearly related to 
‘thermodynamic forces’. These are driving forces responsible 
for the transport processes. 


For example, temperature gradient is responsible for flow of heat, 
and a gradient of chemical potential is the driving force for diffusion. If J is 
taken as the rate of flow or flux and x as the force, the flux- force 
relationship is of the form 


J=LX 


Linear laws of this kind are called the phenomenological relations. 
At thermodynamic equilibrium for all processes, the forces Xx and hence the 
flows, Jk in the system are zero. 


Consider a system in which a temperature gradient, X; and also a 
concentration gradient X» exist. There will be a flow of heat and flow of 
matter. Let Ji and J2 represent the rate of flow of heat and the rate of flow of 
matter respectively. If these takes place separately. Ji aX; and Jo a X». 
According to onsager, if both he processes occur simultaneously, near 
equilibrium the flows and fluxes are related by the phenomenologicals 
equations. 


Jj LyX + LX, (1) 


J, = L, Xı + LyX, (2) 


In these equations, 


‘Liv’ is the thermal conductivity coefficient which relates Jı and X1. 
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‘L22’ is the diffusion coefficient which relates Jo and X2 


‘Lix and La) the cross coefficients describes the coupling of the two 
irreversible processes | and 2. These are called drag co-efficient. 


‘Li’ represents the heat flow arising from a concentration gradient. 
‘Lai gives the flow of matter in response to a temperature gradient. 
Eqn 1 and 2 can be written as 

J,= L;iX; LX, 

J, = Ly Xj t LyX; 


The coefficients Lik, Lx; Li Lik are called phenomenological 
coefficients Li; Lęk are specialll referred to as the direct co-efficients. 
while Lik and Lxi are referred to as cross coefficients. 


The rate of entropy production is given by equation 
ÓJ,X 0 


Substituting for the Jx from equations 


dis dt = L X°4+(L +L)XX4+L X?»0 
222 


111 12 21 12 


This quadratic equation will be positive if both X1 and X» have the 
same sign and becomes zero when X; = X2=0. 


The phenomenological coefficients must satisfy the following 
conditions 


Li1>0, L22>0 


(+L, «ALL, (b =4ac) 


This can be shown as 
LX^4(L +L)XX+L X^»0 
111 12 21 12 222 


Ly+L L 
L{X?+( ? 7)XX + ?X?)»0 
11 1 1 2 2 

11 Li 
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L 
L{X?4+2aXX ta X*+ ?X?^-gq X?) »0 
11 2 2 


1 1 2 2 
Ly 
Unit - 4 
where, 2a — Aut Heat Capacity of 
Ly soilds 
L L +L 
L{(X +aX Y* ?X?'-( 9? "yx?jo NOTE 
11 1 2 2 2 


L 7779 


11 


yo Apal =r 7 ec) 
Ljt((X, aX,) [—* Ar "—]X,) 0 3 


11 


(X,+ aX y is always positive. From equation (3), it is seen that 


4L L -(L+Ly - 
u 2 2 231  mustalso be positive. 


11 


Otherwise the term within { ) will have a sign which depends on the 
magnitude of X L (ispositive. 


So(L +L)? <4L L 
12 21 


11 22 


Since the term within the brackets { } is positive Li1>0. 
Since (L +L  »0,4L L >0, L »0 
12 21 11.22 22 


Onsager showed that if fluxes and forces are properly chosen, the cross 
coefficients become equal 


i.e Li2= La or Li = Lki —> (4) 


This equation (4) is called onsager’s reciprocal relation. 


4.4 Check Your Progress 


1. Define Heat capacity. 
2. What are all the limitations of debye theory of heat capacity of 
solids. 
67 
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3: Define Non equilibrium thermodynamics. 


4.5 Answers to check your progress questions 


1. Heat capacity is the quantity of heat required to raise 
the temperature of the whole of the substance through 1K. 
Even at absolute zero the atoms in a solid will be vibrating 
about their equilibrium positions. When the temperature is 
increased the amplitude of vibration increases and vibrational 
energy of the atoms of the solid is increased, it is the rate of 
change of vibrational energy of the solid with temperature that 
determine the heat capacity of solids is 


dE, 
( vib’) = C, 
oT ” 
2. Limitations of Debye theory 
. 0p should be independent of the source of 
measurement but it is not. 
° 0p Vs T should be a horizontal straight line since Op is 


a constant and independent of temperature. But there are 
deviations in the case of Ag and Na I. 


e In the case of metals there is electronic contribution to 
the heat capacity due to free electrons, which has not been 
taken into account 

Cv= R(T/T,) , T'r-fermi temperature. 


s Debye theory considers the solid to be a homogeneous 
continous elastic medium without any definite structure. It 
may not be exactly correct. 

3. Non equilibrium thermodynamics 

The branch of science dealing with the study of 
thermodynamic properties of the system which are not in 
equilibrium and involves transport process which are 
irreversible is termed as irreversible or non equilibrium 
thermodynamics. 


68 


4.6 Summary 


Heat capacity is the quantity of heat required to raise the temperature of the 
whole of the substance through 1K. 


Einstein theory of heat capacity of solids: Classical theories failed to 
explain its behaviour of variation of atomic heat capacities of solid elements 
with temperature due to the assumption that solid elements absorb heat 
continuously. Einstein used Quantum theory to explain all the 
experimentally observed regarding the variation of heat capacity with 
temperature. 


Debye theory of heat capacity of solids: Einsteins theory fails to describe 
adequately the behaviour of crystals near absolute zero. 


In order to determine the distribution of frequencies Debye disregarded the 
structure of the solid and he treated it as a homogeneous continuous elastic 
medium the vibrations of atoms are considered as equivalent to elastic 
waves propagated in the elastic medium with different frequencies. 


Non equilibrium thermodynamics: The branch of science dealing with the 
study of thermodynamic properties of the system which are not in 
equilibrium and involves transport process which are irreversible is termed 
as irreversible or non equilibrium thermodynamics. 


Onsager reciprocal relation: In thermodynamics, the Onsager reciprocal 
relations express the equality of certain ratios between flows and forces in 
thermodynamic systems out of equilibrium, but where a nation of local 
equilibrium exist. 


Reciprocal relations occur between different pairs of forces and flows in a 
variety of physical systems. 
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4.7 Keywords 


Heat capacity 

a) The constant-volume heat capacity can be calculated from 
the molecular partition function. 

b) The total heat capacity of a molecular substance is the sum 
of the contributions of each mode. 

Einstein used Quantum theory to explain all the 
experimentally observed regarding the variation of heat 
capacity with temperature. 


In order to determine the distribution of frequencies Debye 
disregarded the structure of the solid and he treated it as a 
homogeneous continuous elastic medium the vibrations of 
atoms are considered as equivalent to elastic waves propagated 
in the elastic medium with different frequencies. 


Non equilibrium thermodynamics: Thermodynamic properties 
of the system which are not in equilibrium and involves 
transport process which are irreversible is termed as 
irreversible or non equilibrium thermodynamics. 


Onsager reciprocal relation: In thermodynamics, It’s 
explains the the equality of certain ratios between flows and 
forces in thermodynamic systems out of equilibrium, but 
where a nation of local equilibrium exist. 


4.8 Self-assessment questions and exercises 


1. Define Heat Capacity and explain the Debye theory of heat 
capacity of solids. 

2. Explain the Einstein theory of heat capacity of solids. 

3. Derive Onsager reciprocal relation. 
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4.9 Further Readings 


1, Statistical Thermodynamics, M.C. Gupta, Wiley 
Eastern, New Delhi, 1990. 


2, Introduction to Statistical Thermodynamics, 
R.P.H.Gasser and W.G.Richards, World Scientific, Singapore, 
1995. 
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Unit-5: Quantum Chemistry 


Structure 

5.0 Introduction 

5.1 Objectives 

5.2 One dimensional harmonic oscillator 

5.3 Rigid rotator 

5.4 Check your progress questions 

5.5 Answers to check your progress questions 
5.6 Summary 

5.7 Keywords 

5.8 Self-assessment questions and exercises 


5.9 Further reading 


5.0 Introduction 


There are several areas of chemistry that require a knowledge of 
quantum 

Mechanics for their explanation and understanding. Therefore 
quantum 

Mechanics at an elementary level is covered in several physics 
and chemistry.We will introduce some of the procedures and 
terminology of quantum mechanics and some of their 
applications.Here complete coverage of quantum mechanics and 
its applications given below. 


5.1 Objectives 


After going through this unit, you will be able to: 


e Understand about the one dimensional harmonic 
Oscillator. 
. Understand the concept of rigid rotator. 


5.2 One Dimensional Harmonic Oscillator 


Vibration of a diatomic molecule can be described by a 
harmonic oscillator. Consider the simple harmonic motion 
of a single particle of mass m attached to a point and 
oscillating to and fro about its equilibrium position. 
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Restoring force (f) acting on the particle is proportional to Unit -5 


the displacement x from the equilibrium position. Quantum 
Chemistry 
fex 
f--Kx NOTE 


This expression is called Hooke's law. A force of this type is called 
"harmonic" and the system obeying Hooke's law is called simple 
harmonic oscillator. 


Potential energy V of such a particle is equal to the integral of this 
force over the distance it acts 


V= E E -[C Kx).dx 


MES 
C24 


Potential energy V= + 
2 


Total energy of Harmonic oscillator H =f +V 


ie S 
“| 87 aS u 


Schrodinger equation representing the system of one dimensional simple 
harmonic oscillator is 


E h? n 


8z^m ax” 
upon dm" 
sam 1 W)= 3() 
2 y+ E Kx ] x O0 1 
ox? k 1 2 
8z^m 82°m K 2 
" a= E um p p= 87 = 
et he ` he 2 , 2h 
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NOTE 


d'y oy sam yO- I FP yc) 


"m m x "m x 0 
Ola- £? *)w()» 20 
Ec x x 0 2 
dx? 


Changing the variable to 


awk 
VB 
vas, 
B 


Equation (2) becomes 


A vO) Ig BS) feb o 


wip \ p) dp 
Dividing through out by P and multiply by VB 
^ (9 -E e= Q9 


dé V 


Eigen function: 


V (£) is the solution of this equation. But it is not easy to 
find the solution of this second order differential equation. 
Hence this equation is compared with the second order 
differential mathematical equation. 


Pu Nep 


X ux 
dà 


It is the mathematical equation. The solution of this 
equation is u(x’) and contains a polynomial in x`. 


x 2 


-e?.H,(x') 


where H;(x') is called Hermite polynomial of degree n and has real 
roots. 


Comparing equation (3) with (4), unnormalised solution of the 
equation representing simple harmonic oscillator system in its n" 
quantum state of vibration can be written as 


y(&)-e T n) 


where H,(&) is called Hermite polynomial of degree n and has real 
roots. 


f= Bx . 
v()-* n Qr) 


Normalization of the wave function for the oscillator is carried out as 
follows 


_ 2" nw 
Wn (EYLE) =N B 


[- +o] 


Hence normalized solution for the harmonic oscillator is 


ae nm 


("nt ja) 


volta y eur e 


"o [zaWa] 


H,(\Bx)=H, (ve) 


5.3 Rigid Rotor 
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Rigid Rotor consists of two spherical particles attached 
together and the particles are separated by a finite fixed 
distance. 


Consider a rigid rotor model consisting of two 
masses m; and mo at a distances rı and r2 respectively from 
the centre of gravity of the system. Let *v' be the distance 
between the two masses and this distance remains 
unchanged during rotation and hence it is a rigid rotor 
model. 


Moment of inertia of a rigid rotor: 


Moment of inertia of the molecule about an axis 
perpendicular to the plane of the molecule is 


psy m" 


To get ‘T in terms of ‘r’. Consider the following quantities. 


r-r tr, 
r, =r —rn 
my mh 


I= mrri mrn 


=nn(m+m) 
mr, -myj-mjr-r) 
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r potih 
2 
m+ m, ) (m+ m; 
m m? * nim 2 
[= 2y 
m Fm 
2 
2 
mm p) 2 
m m+ p ur 
T= atm jy t= r 
2 (m, * m, ) 


Schrodinger wave equation for rigid rotor 


Total energy H = Kinetic energy + Potential energy 


H=T+V Potential energy ‘V’ is taken as zero 
Hence H=T +0 


And H =T 


Let v; and v2 be the velocities of the two particles of mass m; and m2 
respectively then 


Kinetic energy of the rigid rotor 


Linear velocity ‘v’ is related to angular velocity ‘œ’ as w= v/r , where ‘r’ 
is the radius. 


v=@r’; v=0r 
1 2 2 


1 
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where L is total angular momentum 


A, 


2 


Lo es ae ae 
237 for a rigid rotor 


^ 


Schrodinger wave equation for a rigid rotor is 


Hy(0,¢)= Ev(0.9) 
i v(69)- vg) 


2I 


In spherical polar co-ordinate system 
P=-Wf 1 M LA 1 a 
a Ly (op 
4x3 sin 090 36] sin ? ZI 


Hence shcrodinger wave equation becomes, 


“a a alig e | & Yw(6.6) - Ev(0.9) 
Am "eem 98 sin? 000) | 


upon rearranging we get 


Ca asse? y(ag)'+ | y(ooyleserEy(0.0) =0 


NERA 
umm 98 sin" 006 ) m 
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This is Schrodinger wave equation representing rigid rotor. 


This equation contains two angular variables, 0 and ® and is 2" 


order partial differential equation. This equation is separated into two 
equations, which are total differential equations by using the separation 
of variable procedure and solved. 


Separation of variables: 


Total wave function is written as the product of two wave functions, 
which are dependent separately on only one variable and independent of 
the other and inserted in the equation 


y (8.0) = 66d, 
(1 O9fsne948,6,l- L P 8p +871 
Me ag doo sn^899 Ao UE Eo) Ko) =0 
N 
(5 Os. g9 p Y, ko F 8271 
"d ? 38^ + in Aad? fo EG e)%o)=9 
vea 
Multiply by 5 9 
(oto) 
sin il d ) 1 æ% mI ü 
38 sin 8 Ae) + a2 ae Esin 0-0 
ES Àj o 
sin 0 2 272 \ SmI 7 1 9 4 
sin (0) |+ Esn 60--— 2 (9) 
m 98, 00 ) n oF 


Left hand side of this equation contains the function of 0 and right 
hand side contains the terms, which depends only on the variable ®. 
Both sides are equated to a common constant m’. 
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This is called 0 equation. 


1 0 


Gor =m > (2) 


This equation is called ® equation. 
Solution of ® equation: 


Upon rearranging equation (2) 


9? 
et) = m Ko) 


9 
apt + Mho =0 


This equation is called ® equation. 


Øp) is the solution of this equation. 
_ time 
Kp) = Ce 


This is an acceptable solution provided ‘m’ is an integer 


he) =A +27) Since dis an 


angle. 
PELLI = erin(o+ 27) 
gin? — gin) _ etim(2z) 
et" —] 
Note: 
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(e^ — cos Q-F isin 8) 
cos(2mz)+ isin (2mz) =1 
This is true only if ‘m’ can take up 0, +1, +2, +3 etc. 


This function øg) should be normalized. 


(hol ho)! 
[0, 27] 


20 


|  (9)45d9 —1 


) 


(Ao = Ce) 


20 


| Cetin? Ctd ó-1 


2m 


C [a6 -1 
Cp" =0 
C? (22-0) =1 
phar 
2m 
zh 
2m 
1 


Normalised solution of  óequationis $— gue 
V2 
0,1,2,3 


Solution of 0 equation: 
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where m = 


The rearranged 0 equation is that is equation (1) is multiplied by 0(0 and 


divided by sin?0. 


2 


m? ) 
-sin?. 


\ ) 
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0 equation is compared with Legendre equation in 


Hum mathematics. Associated Legendre function is the solution 
Chemistry of the mathematical equation. The normalized solution of a 


equation is 


NOTE axii-my "6 ) 
057 P, cosQ 
2(1-- m)! 


P” (cos@) is a polynomial in cos 0 of degree | and order [n]. 


Hence normalised solution of the schrodiner wave equation 
representing rigid rotor is 


v(80.0)- Go) - Xs) 


ltm itm ltm 


- 21-1 AI — mJ pm (cos 8)—— e+”? 
2(I m) ' J27 


where ‘I’ is a positive quantity. Two quantum conditions 


are 
120, 1,2, 3. erre 
m = 0, 41, +2, 23............ 41 


This solution is called spherical harmonics and is 
represented as 


Y(0.0) = Q6).Go) 


I+m 


A few spherical harmonics are 


Unit -5 


1 0 Y (6.6), = 24 cosa Quantum 
Chemistry 
Energy of Rigid Rotor: | NOTE 


Comparing the two equation 


(1- E —2x 9 u(x) + (C +1)- 0 Nu(x)= 0 


ox? Ox N 1- x? J 


when x=cos@ 


we get B= (1+1) 


(ree 
rot 8 T 27 


where | = 0,1,2 


1. Energy of the rigid rotor depends on the values of | and 
independent of ‘m’. 

2. Energy of rigid rotor is quantized. 

3. Energy increases with increases in ‘I’ value. 


5.4 Check Your Progress 


1. What is meant by Hooks law? 
2. What is wave function? 


5.5 Answers to Check Your Progress Questions 
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1. Hooke's law is named for Robert Hooke,According to the 
Hooke's law ,for a particle executing simple harmonic motion 
in one direction, the restoring force F is directly proportional to 
the displacement x from the equilibrium position 


Fa-x 


F — —kx 
The negative sign indicates that the displacement x and the 


restoring force F are in opposite directions.K is proportionality 
constant called the force constant. 


2. The eigenfunctions of a quantum mechanical operator depend 
on the coordinates upon which the operator acts; these functions 
are called wavefunctions. 


5.6 Summary 


* According to classical mechanics, the state of a system is 
specified by giving the position and velocity of every particle in 
the system. Consider a single particle withoutany internal 
structure so that it cannot rotate or vibrate. 

e [f it can move in three dimensionswe can specify its position 
by the three Cartesian coordinates x, y, and z. 

e These three coordinates are equivalent to a position vector that 
reaches from the origin of coordinates to the location of the 

* particle. 

e The Cartesian coordinates x, y, and z are called the Cartesian 
components of the position vector r.A vector can also be denoted 
by listing its three Cartesian components inside parentheses, as in 
Q, y, z) 

* In order to show how Newton's laws determine the behavior of 
a particle, we apply them to a harmonic oscillator, which is a 
model system designed to represent a mass attached to a 
stationary object by a spring, A model system is designed to 
imitate a real system, but is defined to have simpler properties so 
that it canbe analyzed more easily. 

* Let the horizontal coordinate x of the mass equal zero at its 
equilibrium position and assign it to be positive if the spring is 
stretched and negative if the spring is compressed. 


5.7 Keywords 


e Harmonic oscillator: A System which,when displaced from its 
equilibrium position,experiences a restoring force proportional to 
the displacement according to Hooks Law. 
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° Rigid Rotator: The rigid rotator is a simple model of a rotating Unit -5 
diatomic molecule. Quantum 
Chemistry 


5.8 Self-assessment questions and exercises NOTE 


1. The frequency of vibration of a 1H35CI molecule is 

8.966 x 1013 s-1. 

a. What would the frequency be if the chlorine atom were infinitely 
massive? 

b. What would the frequency be if the hydrogen atom were infinitely 
massive? 


2. What do you mean rigid rotator? Write a note on rigid rotator in a 
plane. 

3. Apply Schródinger wave equation to a system of rigid rotor and 
solve the equation 


5.9.Furthur readings 


1. Quantum Chemistry, LN. Levine, Allyn and Bacon, Boston, 
1983. 

2. Quantum Chemistry, R.K.Prasad, Wiley stern, New 
Delhi,1992. 
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6.0 Introduction. 


In this unit we can elaborately understand Quantum mechanical treatment 
for radial and angular wave function and hydrogen atom like atoms. 


6.1 Objectives 


After going through this unit, you will be able to: 


° Understand the concept of Quantum mechanical treatment. 
e Understand pauli’s exclusion principle and slater determinant. 
° Learn about the Approximation methods. 


6.2 Quantum mechanical treatment 


Hydrogen atom is made up of a proton and an electron revolving the 
nucleus. 


The potential energy of the electron is due to the electrical force 
of attraction between the electron and the nucleus. 


Ze 


r 


V= 


Therefore Schrodinger equation for the hydrogen like atom can be 
written as 
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8z^m 
Vy - 
y+ —_(r-v)ye0 Unit -6 
Quantum Mechanical 
V’ y+ 8m E Ze? Yy- 0 Treatment 
h? ( r rj NOTE 


This can be solved by transforming the Cartesian co-ordinates x, y, z 
into spherical polar co-ordinates r, 0, cb. 


2 z 2 =V? = laf r ð \ 19 ld ( 894 1 d2 
zt Jae Aaa? 
ax” ay’ "E r a xl r? sin oad 36) d x2 


REL a ee (in OW), | Fy 8mm ze lw=0> (1 
r Parl 2 Papua 30) r sin? 899? h2 | :J 


Separation of variables: 


y= R(r).0(8).& 6) (2) 
Substituting (2) in equation (1), we get 


Lal i9? ie à( 9 Odo a EAA yg) 


Aa ze ar ^s gag "^ "n a | r sin? 0/09? 
Multiply by r°, 
o 
(09/773 c), nolis Lfd gto 


8z^mr? e) 


+— l E+ — IR(r).8(8)&(6) =0 
k Ao r) 


Divide by R(r).0 (0). AØ) 
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2 
|o Apa qug) Oeo «e 

Unit -6 RC) ar Pr a AA)sin 030” o A d)sin2 099? 
Quantum 8z?^mr? (E Ze? }=0 
Mechanical Th p + " px 
Treatment 
NOTE 

Multiply by sin?0, 


sin? 00 (59 (9s sin 00 ( g? e(6) ++ 9 Md) 
WO Vra tr "aae ae) agar 
QU mr? sin Ye zezo 


sin? 99 ( ,0 3 singo( d Y 82mr’sin?@(  Ze?\ 1 Fø) 
R(r) |+ rr EY 0(0)|+ — ——| E+— l= 


ET ) €(8)90 j R rc) «0 ag 


Left hand side of the equation contains the variables r and 0 
and right hand side contains the variables ®. Each sides of the 
equation is equated to a common constant say m? such that 


? &0)= 
Hd) 29 


9? 2 E 
ag ht m'(9)-0 > (3) 


The equation (3) is called ® equation. 
sin LIF 2 ð R(r) + sin 0 E (sin 0 ? 9g). 8z^mr? sin "dr Ze jem Ss (4) 


R(r) a D ] PORN 30 h? = r) 


Divide the equation by sin^0 and upon b d 


| 9 [ 2 9 FON l E an 0 j«9s 8z^mr? (E Ze T 
R(r) ar ar ar ] POSTE 00 J h2 7 sin? 0 
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1 af d a RCO)? +82°mr 
ROJU or h? | r 


(E Zeài.m. 1 Ilma? olo) 
+ 
sin? Osin 090 00 J 


Left hand side as well as right hand side of the equation contains only one 
variable and hence both must be equated to a common constant p. 
2 lod 


m — — —7 sin 0 
— A 0(0)- 
sn^9 0(0)sn 090 90 (yep 


-—L— 3 sin 02 
0(0)sin 000 Pgh UE Prin? 


Multiplying by @(@) and rearranging 


il ° sin 6? 0(0)+ I(B- ™ Ye(6)-0— (5) 


snQ 00 00 ( mr 


This equation is called 0 equation 


1 dfa a RC)? m. Ze2\| = 


R(r) or ar\ a 2 
Multiplying by R(r)/r 
“1a 23 RC) -Ê RO) 8 mle ze W(r) =0> (6) 
r? oe oF r h? k r EE 


This is called radial equation or r equation. 


Solution of ® equation: 


® equation is 
d? 2 7 
age miko) =0 


The solution of ® equation is 


timó 


= Ce 


+m 
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(Koao) =1 


[027] 


timó timó 27 timó timó 
Ce /Ce - | Ce Ce =l 


[0,27 | 
2x 


-C' [ao-c'l =1 


Solution of 0 equation: 


The 0 equation ls 5 
sin 0 e(0)-1,p-" > Ye(6) =0 
sn900 00 | sin? 6) 


Associated Legendre function of degree | and order m in the 
variable cos 0 represented as P"cos Gis the solution of 0 


equation. 
Normalised solution of 0 equation 
E — m) 
8(0)- [Emm] pn og 
2(I-m) ! 
l+m 


The product of solutions of 0 and ® equation is spherical 
harmonics that is 
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Ys (8 ó)- A). é 9) 
ltmltm 
Solution of radial equation: 


Radial equation is 


lof 5 ə R(r))- P riya 82° (E Ze jo 0 
r’ SL de or al p 


The solution of this equation is R(r). 


R(r) (E 2c Yn(n- 1-1), w p 


(na, J 2n(n +1) 


where _ 2zr k . and as () is the associated Laguerre 
pP ^ Ay =- n4 P 
Nao Ane Ml 
polynomial. 


Hence the total solution of the wave equation is 


V(r.&.)- R().e(0).90) 


n,lzm  n,l lxm l+m 


v(r.0,9) = [2zY4G- = Ile? 2 p (p) (ENE mp cos |e 
nl tm (na, ) 2n(n +l) n+l (1-4 m)! Jon 


6.3 Hydrogen like atom 


The co-ordinates of electron of Helium atom can be represented as 


Electron 


2 


Nucleus fz 


Hamiltonian operator of Helium atom is 
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For the ground state of Helium atom, the electronic 
configuration is 1s’. The overall wave function will be given 
by a slater determinant. 


Wo 


E 


this on expansion may be written as 


y -—s 15 [a B -Ba | 
42 


According to variation method, the energy, E based on the 
approximate wave function is always higher than true energy 
Eo of the ground state. 


E (y, HIY) 


HS /S 


2 


1 
= i. 1S 


1) 


(2) 


E = (Solse 


0 


1 


Q) 


Q) 


Q) Q) 


Kcu 


qQ) 2) 


H 


Hamiltonian for Helium atom is 


E 


e-siss 


( 


15, 15 
1) 


(2) 


H1 
(1) 


1S 
(1) 


90 


"Du 


Q) (2) 


1S 1S 
Q) 


1 — 


Q) (2) 


Ba la p 
() Q) 


H So ISo) p 


Sy Se ) 


d S 
Hay + Hoy " «Sey 


H 
ol ol 


Q) (2) 


1S 1S 
() 


E 


(1) (2) 


ol) (| 


1S 1S 
Q) 


(2) 
n» 


1 
| 
| 


S 1S 
Q) 


Q) 


| 


E- (so lHo NS MiSolSoh+ 15/818: sq] 18d) 7 


E-Egu(jtEgsQtJ 


J= (bq) As; 


1 
PASOS) 


J is called coulomb integral. 


In a two electron atom, the average potential experienced by an 


electron varies between — _ 


2 E : 
and - that is it depends on the effective 
r. F. 


nuclear charge ‘Z’. The best value of Z will be determined by the 
variation principle. 


The wave function of 1S orbital appropriate to this effective nuclear 
charge Z is 


The result of calculation then leads to 
E--Z? — 2(2- Z') Z'+5/8 Z 
=Z"-27/8 Z' 
The value of Z' which minimizes the energy is found by differentiating E 


with respect to Z’ 


Setting = =0, we get 


az'- 2! <9 
8 

Z=- 160 
16 


S E = —2.8476(a.u) 
(1a.u227.2ev) 
E = —1148eV. 
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6.4 Pauli's Exclusion Principle 


Pauli’s exclusion principle states that “The total wave 
function of electron must be antisymmetric under the 
interchange of any two electrons or no two electrons in one 
and the same orbitals can have all the fact quantum numbers 
same". 


The four possible product combinations of the orbital and 
spin wave function for He atom in Is? state is 


(is, 15, Xo) at) SxS=S 


(18,18, \By/Bo)) SxS=S 


(is, 1S \La B +a B) Sxs-s 


1) (2) 13 (1) Q) (2) o 
1 
(is, 1s )--(« B -a B ) SxA=A 
1) (2) J (1) 2) (2) (1) 


Symmetric (+) x symmetric (+) or 


Antisymmetric (-) x Antisymmetric (-) is symmetric (4-), 
but 


Symmetric (+) x Antisymmetric (-) is Antisymmetric (-) 


Complete description of an electron is given by the spin 
orbital wave function 


P(r, s) = P(r) n(s) 


@(r) is the orbital wave function, which is a function of 
the position variable and m(s) is the spin function involving 
spin variable. 


Spin orbital wave function is obtained by multipling 
each of the orbital function with those of the spin functions. 


Ts) = P(r) n(s) 
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1 (9). 929 (1).9(2))(a a ) 


Eo 1 Q (2) 
SxS =S 
= 1 (9). (2)+0(0). eC). ) 
' m 1 2 1 (1) (2) 
SxS =S 
v= 1 (0090)«0()90). (aB «Ba ) 
J2 * 1 Q) (2) (1) (2) 
SxS =S 
w= 1 (0). 9(Q2)«00. 90). (as -Ba) SxA =A 
ae. 1 Q (2) (1) (2) 
v= 1 (6). o(2)- (1) Da E ) AxS=A 
J2 1 
y= i (1). «)- e) QU 5 ) AxS=A 
4C! 
y- 1 (9 «()9Q)- 940). (ap +Ba ) AS SUA 
42! )0 (Q00 
=i (O. «()- eQ. 42) («B -pa ) Ax A=S 
Eu (1) (2) (1) (2) 


Of the eight possible combinations only four combinations are 
antisymmetric and they are accepted according to Pauli’s exclusion 
principle. 


y (.s)- 1 OMAI (000) (a8. -pa )> Q) 
A J2 1 2 2 1 (Q 09 @@ 
v (r s)= 1 0().9Q)- eoa i r (2) 
42! 2 
y G.s)- 1 (4(0)90)-60)90)(5.8. ) G) 
A J2 1 2 2 1 a) (2) 
v Cr, s)= 1 (9). 92)- IDAD (ab «Ba )>(4) 


X3 © (2) (1) (2) 


To prove Pauli’s exclusion principle consider equation (4) 
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| y (e= 1 (0.62)-0AD@B +B a ) 


Unit -6 1 2 1 (1) (2) (1) (2) 
Quantum 42 
Mechanical This can be written in the Slater determinantal form. 
Treatment 

NOTE 


1 OIO e + p). Po 2 AG a a d ABa ) 


1 (1) (2) 


Sl- 


70s (3-669 Q)- 6(0:96)—96).0 Q)) 
2 1 2 1 


yo.) = 90 E (0 aC 
ZEE ko ez em c] 


When both the electrons are in the same orbitals 


(o, = d». Therefore the determinant becomes 


NOE pn oA 3 A 
42 øl) o,(2 
ó 


when a = p ) 


it t t ad 


When two rows of the determinant or two columns of a 
determinant are equal then the determinant vanishes i.e. it is 
equal to zero. The wave function also vanishes. Hence no 
two electrons in one and the same orbital can have all the 
four quantum numbers same. Pauli's exclusion principle is 
proved. 
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6.5 Slater determinant 


In order to construct antisymmetric wave function for a multi-electron 
atom, the increase in the number of electrons increases the number of 
terms in such large proportions [for ex: 5! i.e 120 terms for an atom with 5 
electrons only]. We must find an abbreviated form to represent a wave 
function. This can be done by writing it in the form a determinant in 
which the spin-orbitals are the elements, each row in the determinant is 
labeled with an electron and each column with a spin orbital. 


The normalized wave function for a three electron atom is written as 


_ 1/|Ag Bo Sw 
y = (1) 

NEG Bur Co 

— 2 2 2 

6143) Ba Co 


For an n electron atom 


The antisymmetry of equation (1) or equation (2) as well as Pauli's 
principle are guaranteed by two properties of determinants. 


(1) If two rows (or) two columns of a determinant are interchanged , 
the resulting 
determinant is just the negative of the original one, i.e 
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Ay Bay Co Ao) Bay Co 
Ao Bey Col=- lAa Ba Cw 
Ag Bey Co Ag) Bay Co 


(2) If any two rows or two columns of a 
determinant are the same, the determinant 

vanishes. 

For example: 


The ground state antisymmetric wave function for He 
atom in the determinantal form and to prove that the two 
electrons cannot have the same spin-orbital. 


(o... = ao) 


v- =[S 01S) 


= [Skobe] Sob Soe 


S0 — 15051; 
185,05 Sobo) 


> (3) 


where 1S, a) is the spin-orbital of electron i (i = 1, 2) 


If both electrons have the same spin-orbitals, the determinant 


215099 15960. 9 
15505. 15:5 


With increasing number of electrons, even the determinant (2) 
will have large size and further abbreviation is necessary. This 
is done by using a bar over the orbital which has f spin, 


96 


absence of a bar will indicate a spin. Thus the determinant 
wave function for Li atom is 


= RET (2) (2) (2) 


[S ISo 256 


6.6 Approximation methods 


Need for approximation methods: 


1. Schrodinger wave equation cannot be solved exactly in many electron 
atoms because the presence of more than one electron introduces electron 
repulsion terms (e?/ri), where ri; is the distance between the two electrons. 
This rij term depends on 6 and 6 coordinates. The separation of variables 
is not possible and is very difficult. 


2. Again the Schrodinger wave equation cannot be solved exactly for one 
electron systems, whose potential field is not spherically symmetrical e.g. 
Hydrogen atom subjected to electric or magnetic field. 


Hence to obtain the solution of the wave equation for the above 
systems approximate methods are used. The two commonly used methods 
are 


(1) Variation method 


(ii) Perturbation method 


6.7 Variation Method 


Variation method is an approximation method. It is used to those systems, 
which differ much from the systems, for which exact solution is known. 
In this system the wave functions can be guessed based on physical and 
chemical combinations and the energy of the system is calculated by 


(Hy) 
E- din d when vis normalized, then 


Wy) 


E = (yH) where H is the complete Hamiltonian of the system. 
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Variation theorem states that “if y be the normalized trial 
wave function satisfying the boundary conditions of a 
system, whose Hamiltonian H has a discrete eigen 
spectrum". 


Then 
( v Hy) 2 Eo 
where Eois the lowest exact Eigen value of H. 


The wave function wis called trial variational function and 
the integral (yH) is called the variational integral. The 


energy calculated by variation method will always be greater 
than the true energy Eoi.e. the lowest eigen value of H' . 


Proof of variation theorem i.e. proof for (uy Hy) 2 Eo 


Let 00,01,02,03.....be the set of normalized and mutually 


orthogonal eigen functions of H^ with discrete eigen values 
Eo,E1, E2, E; ,etc. such that 


Há; = E; ģ; 1-0, L1; 2-3, 5 


Let ybe a normalized trial variational function. Let y be 
expanded in terms of the orthogonal set of functions i.e. y= 


CÓ 
Multiply by y“ and integrate over all space, 
(y- ce», cá) 
= Gl) 


(w|w)- 3, C^ 21 [~y is normalized yy )- 1] 


Considering the variational integral 
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E- (By) 


Unit -6 
E= (È Có Hy, c, ) MADE 
E- oT C,9,+ C0, Hf | C. CØ, + C40, - NOTE 
E=), COROA, CÉ HP i) > (1) 
Moreover, (@H|¢) = Ei — (a) 


(ole) = Ei (Øg) - 0 > ® 
Then (1) becomes 


E 2 C^E -0 [The second summation in (1) vanishes due to (b)] 


If Eois the energy of the lowest state or lowest Eigen value of H then 
E-Eo- 245 a 25678, [5 ae = 1) 
E-Eo= >),C'(E,- &) > (2) 


The quantity £,—£, must be positive or zero for all values of i and Ci? is 


always positive, the R.H.S of the above equation (2) must always be 
positive or zero. 


Hence E-Eo2 0 
E 2 Eo 


(WHY) = Eo 


This is variation theorem. 


6.8 Time Independent Perturbation 


Perturbation method is an approximation 
method used to solve the wave equations representing the 
systems having more than one electron. This method is 
suitable for a system, which are subjected to very weak 
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electric or magnetic field and also for the system whose zero 
order states are non-degenerate. 


Conditions under which perturbation theory is based: 


1. The total Hamiltonian of the 
perturbed system is the sum of the parts (i)unperturbed 
Hamiltonian H and (ii) perturbation 


H=H+Aav 


where AV is the perturbation and / is the expansion 
parameter. 


. The system is only very slightly disturbed and different from 


the system for which the exact solution is known i.e. AV is 
very small compared to Å ® . 
A OAV 


. The eigen values E? and eigen functions y of the 


unperturbed Hamiltonian H® are known. 
H y= EO yO 


. If wis the wave function of the Hamiltonian of the perturbed 


system H with eigen value E then the eigen value equationis 
Hy =Ey 
(H+AV)w=Ew 


The wave function yand energy E are functions of 2 and 
hence it is possible to expand them in terms of / in the form 
of power series as 


w=y® + AyD + PYP RYO | | 2| (1) 


And 
E=E e AE - VEO e EO + (2) 


Where y”, y®...... and E®, E? .... Etc are independent of A. 
y and y? are chosen such that they are orthogonal toy, 


which is normalized. 


5. 'The equation of the perturbed system is 
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(H? - AV (YW e Ay? - y? +...) S (E? AE 4 VEO +...) y HAYO yr? +.) 

The order of accuracy depends on the number of terms chosen in the 
above equation. Since the perturbed system is only slightly different from 
the unperturbed one. y and EO may be chosen as approximate solutions 
of the perturbed system. This is called zeroth order approximation to the 
true solution of the perturbed system. 


If the second term in the series that is y.2 y ? + Aw and 


E= E ? + AE ™ are chosen and the schrodinger wave equation is solved it 
will give rise to first order approximation to the true solution of the 
perturbed system. 


Ay = Ey) -5Zeroth order perturbation equation. 


( 7) — EO) =Vy) + Ey) — First order perturbation equation. 


First order perturbation theory: 
First order perturbation equation is 

(AO — ply = py + Boy 
By solving this equation first order correction to the eigen function y? 
and the first order correction to the eigen value E ? are obtained then y 
and E are calculated. 


yay + Ay” 
E = E? +E” 


y and E® are known for the given system under consideration. 


First order correction to the energy and approximate energy of the 
perturbed system: 
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(4 Oy Ey ) 2 py) +E! T ) is the first order 
perturbation equation for the i" state . y p is expanded in 
terms of the complete set of eigen functions of the 
unperturbed Hamiltonian H © . The unperturbed equation is 


fj oy, = g 0), O 


L L l 


where i = 0,1,2,3 etc. 
yo YW... y y etc belong to orthonormal set of 


functions and E, is non-degenerate. 


Ta Lew” — (a) 


j#i 


(0) 


Multiply by y`” and integrate over all space 


(yi y m = 2c lve Y= 0 


Since the wave function m is orthogonal to v ©. Hence 


y is orthogonal to y? 
i j 
Substituting the value of y O in the first order perturbation 


equation we get 


(47° -po Cy = a vy" +E" y o 
jai 


$6. (H° = Ey J= vy Ey 


j#i 


(0) y,,(0) — FO ,,, (0) 
H yY, =E; Y; 
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bye (E? _ Eo ) y O — md y ® FE Eo yr O 


je Unit -6 
. . Quantum Mechanical 
Multiply by w/? and integrate over all space Treatment 
NOTE 
Ycap go - Eoo )-- (fv y+ E yy 


J*i 


ECE -EO - feme lolo) 


y ? and y etc belong to orthonormal set. Hence 
i j 


yo y — 0 and (Wy 1 


Yu 9+ go 
EO = po Vu 


Hence first order correction to the energy of the i" state is 


0=-y 


EP = oy 


E (? is the first order correction to the true energy. The approximate energy 


of the perturbed system according to the non-degenerate time independent 
first order perturbation theory is 


E = EO + Ag? 


serae 


First order correction to the wave function: 


First order perturbed equation for the i" state is 
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[i - y? vt sey? 


y is expanded in terms of a complete set of eigen 
functions of the unperturbed Hamiltonian 


v" = LCY 


jsi 
yw and y " belong to orthonormal set. Substituting this in 


perturbation equation 


(4° - go Cy) = - Py + pto 
i ij i i i 
j*i 


(0) y 0) — (0) 45, (0) 
HOUR Ey, 


Y CIE E Eo y O - ay y; O + go y O 


j#i 
upon expanding the summation we get, 


putre 1.979 de osito for the i" state. 


(0) (0) (0) (0) (0) (0) 
E" -E + - 
e 1 i Nn C (EO — EY) +C EL -EO 4 c. U^ E Le i 


T d y O 3 go y O 


(0) 
1 


In order to obtain C multiply both sides by y" and 
il 


integrated over all space 
qr - gy ppm. - o veo. Em uy ) 
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Since the zero order wave functions are orthogonal and individually 


normalized 


eer retten da 


z CE 47 i 
= ( Eo _ E) 


il 


- (vp) 
( Eo — go ) 


i2— 


(0) 3, (0) 

" vy) 
jo (0) (0) 
(E; us ) 


Vy") 


Na 
7. (go po 


Substituting the value of C in y ^ 


ij i 


(0) -(0) 


(= T 0) 
yc (EO BOY 
j i 


i*j 


y is the first order correction to the true value function. 


EO + EO 
j i 


Thus according to the non-degenerate time independent first order 
perturbation theory. 


Approximate solution of the first order perturbation equation. 
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(9) (0) 
v ax eu 


izj (E zs E) : 


This relationship is not valid if there is degeneracy in the 
eigen spectrum of H ® that is E — E for i + j . For it to be 


valid ground state should be non-degenerate. 


6.9 SCE Methods 


According to Hartree each electron is assumed to move in 
a spherically symmetrical field due to the nucleus and the 
average potential field is due to all the electrons except the 
one under consideration. This field in which electron move is 
called self consistent field. 


The wave function of an atom containing ‘n’ electrons is 
written as product of ‘n’ number of one electron functions that 
is 


v-ó0)6Q)8G) 6) 


where $,,05.0, ó, are the normalized and mutually 
orthogonal one electron wave functions. 


Ó (1) means the electron (1) is put in the orbital ó................. 
etc. 


The antisymmetric wave function of the atom is expressed 
as a determinant. 


WAD 60) (0) Q0) ..... $, (1) 


_1 PQ) &€ 62 &9 4) 
"I. 3 


n! 


b(n) "me $5 (n) EE "m Ø, (n) 
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For an n-electron atom with closed shell configuration, that 
is in which all atomic orbitals are doubly occupied with paired 
spins. 


1 
y- Hall) di P" d mw à 


where the determinant has been represented through diagonal 
elements. ¢,(1) is the spin orbital with electron 1, ø, (2) with electron 2, in 
the same orbital ø, and so on. 


Fock introduced a new operator F^ known as Fock operator. It is an 
operator for kinetic energy of electron, potential energy of inter electron 
repulsion and potential energy of electron spin exchange between pair of 
electrons. 


The Fock operator for electron is defined as 


Êf =-ty? -445$ (3 -K) 


(1) 2 ! no ja j j 


1 TE 
—  V' represents kinetic energy operator forelectron 
5 1 


— * represents operator for potential energy of attraction between electron 
r 


] and nucleus. 


J, represents operator for potential energy of repulsion between electron 


1 and j^ 

electron. 
K represents operator for potential energy of change of spin between 
electron 1 and j” 

electron. 


The potential energy of repulsion experienced by electron 1 in the field 
of electron j is given by the following operator J; 
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J- 6 l adr 


Jo d 
lij 


j j 
The sum of the operator for the potential energy of 
repulsion experienced by electron 1 in a field of every other 


electron of the atom is expressed as follows 


The potential energy of exchange of spin between 
electron 1 and electron j is given by the following operator 
Ki 


1 
K= fo, $947, 
lj 


There will be exchange of spin between electron 1 and n/2 
number of electron with opposite sign. Therefore there will 
be n/2 number of such operators. The sum of these operator 


Ir- e 


i à j 
j=l j=l n 


Thus the Fock operator for electron 1 Fry is described. 


Similarly the Fock operator for other electrons can also be 
formed. 


A set of orbitals without spin, 1,0», 03 g, with 


2 
variable parameter in each is chosen. Using these functions, a 


set of Fock operators Fy) : fj ; Fo) : Fy F, were 


formed. Substituting the above wave function with respective 
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Fock operators and so having the following “Pseudo- 
Schrodinger wave equation". A set of wave functions 
Ó .À d sss ó, were obtained 


2 


Fo -EÓ 


, Using the new set of wave functions, a new set of Fock operators 
QE 


Ê "E y aesan F are formed and usingó ,@ ,ó ............. ó , 
D) 0 (9 (2 $ 12 3 7 
the pseudo schrodinger equation is solved to get a set of functions 
D sU c aues o . This process is repeated till we get a set of orbitals 
1 2 3 n 


2 
which are not very different from the previous set. This set of orbital is 
called self consistent field (SCF) orbitals g^% 39 9 ee 
1 1 2 


6.10 Check Your Progress 


1. State Pauli’s exclusion principle 
2. What is meant by Slater determinant? 
3. Write notes on self consistent field method . 


6.11 Answers to check your progress questions 


l. Pauli's exclusion principle states that “The total wave function of electron 
must be antisymmetric under the interchange of any two electrons or no two 
electrons in one and the same orbitals can have all the fact quantum numbers 
same". 


2  Inorder to construct antisymmetric wave function for a multi-electron 
atom, the increase in the number of electrons increases the number of terms in 
such large proportions [for ex: 5! i.e 120 terms for an atom with 5 electrons 
only]. We must find an abbreviated form to represent a wave function. This can 
be done by writing it in the form a determinant in which the spin-orbitals are the 
elements, each row in the determinant is labeled with an electron and each 
column with a spin orbital. 
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3. According to Hartree each electron is assumed to move in a 
spherically symmetrical field due to the nucleus and the average 
potential field is due to all the electrons except the one under 
consideration. This field in which electron move is called 
self consistent field. 


6.12 Summary 


Pauli's exclusion principle states that *The total wave function of 
electron must be antisymmetric under the interchange of any two 
electrons or no two electrons in one and the same orbitals can 
have all the fact quantum numbers same". 

For many problems, it is not practical to obtain a wave function by 
the exactsolution of a wave equation that describes the system. It 
is still possible toperform many types of calculations, and one of 
the most useful techniques is that known as the variation method. 
The exact solution of the Schródinger wave equation for complex 
atomsis not possible. However, examination of the form of the 
wave functions obtained for the hydrogen atom suggests that 
approximate wave functionsmight be obtained if we were to take 
into account the mutual electronrepulsion. Such a procedure has 
been devised by J. C. Slater, and theapproximate wave functions 
that result are known as Slater wave functions [or Slater-type 
orbitals (STO)]. 

Perturbation method is an approximation method used to solve 
the wave equations representing the systems having more than one 
electron. This method is suitable for a system, which are subjected 
to very weak electric or magnetic field and also for the system 
which zero order states are non-degenerate. 

According to Hartree each electron is assumed to move in a 
spherically symmetrical field due to the nucleus and the average 
potential field is due to all the electrons except the one under 
consideration. This field in which electron move is called self 
consistent field. 


6.13 Keywords 


Radial wave function: An orbital in a mathematical function 
called a wave function that describes an electron in an 
atom.Radial wave function for a given atom depend only upon the 
distance, r from the nucleus. 

Angular wave function: Angular wave functions depend only 
upon direction and in effect, describes the shape of an orbital. 
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e  Petrubation: The act of disturbance, in physics a secondary 


influence on a system that modifies simple behaviour, such as the Unit -6 
effect of the other electrons on one electron in an atom. Quantum Mechanical 
Treatment 
NOTE 


6.14 Self-assessment questions and exercises 


1. Discuss the principle of variation method. 
2. What is meant by Slater determinant? Obtain its value for helium atom. 
3. “For many electron systems only approximation methods are needed”- 


Explain. 
6.15 Further readings 
1. Quantum Chemistry, I.N. Levine, Allyn and Bacon, Boston, 
1983. 
2. Quantum Chemistry, R.K.Prasad, Wiley Eastern, New 
Delhi, 1992. 
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UNIT- 7 HMO Method 


Structure 

7.0 Introduction 

7.1 Objectives 

7.2 Application of variations methods to hydrogen atom 
duo Application of perturbation methods to helium 


74 HMO method 

7.5 Application to butadiene 

7.6 Check Your Progress 

da Answers to check your progress questions 
7.8 Summary 

7.9 Keywords 

7.10 — Self-assessment questions and exercises 
7.11 Further reading 


7.0 Introduction 


when the phrase “molecular orbital calculations” is first 
encountered, 

the mental image may well be one of hopelessly complicated 
mathematics and piles upon piles of computer output. It is 
interesting to note, however, that sometimes a relatively simple 
calculation may provide useful information that correlates well 
with experimental observations. Such is the case with the method 
known as the Hückel molecular orbital (HMO theory) calculation. 
This method was developed in 1931 by Erich Hückel, a physicist 
in Marburg, Germany, who was trying to understand the concept of 
aromaticity in benzene. The calculational procedures are relatively 
simple and have become known as the “back of an envelope” 
calculations. 


7.1 Objectives 


After going through this unit, you will be able to: 

. Explain about the application of variation methods to 
hyderogen atom. 

e Understand the concept of HMO method and application to 
butadiene. 
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7.2 Application of Perturbation method to hydrogenatom 


Unit 7 
HMO Method 


electron 


Hamiltonian operator of Hydrogen atom in atomic units 


v-la.3, 1 mg, ! > 


ror or r sn690 00 r sin? 0 ag 


Hydrogen atom is spherically symmetrical and hence the wave function 


has no nodes and y is independent of 0 and ®. 


1 
Hence V^ aa. 

r°or or 
And Hamiltonian H= Ee. 9| 1) 

2 

U 2\r° or J J : J 

The wave function for hydrogen atom is chosen by infinitive means. Le 
yr RI e vt 


where *a' is the variable parameter. 
The energy of the system is obtained from the variational integral. 
Mild 
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a 2! 1! I! 
2 Qay" Qa)" Qay 
= 2i Unit 7 
0 ys HMO Method 
a 
dd C+ -)? NOTE 
E-.9a 4a 4a — a 2a 28a 
2 8a? 
8a? 2 
a? — 2a 
E — 
2 
Energy is minimized with respect to a 
OE _ 
Oa 
díu -2au| 1 
| |=—(2a-2)=0 
da\ 2 J 2 ( ) 

a=1 
Hence the wave function is 

y=e" 
Energy 

a’ -2a 
E= 
2 
E= ECT = m au 
2 2 


This represents the exact value of energy and exact eigen function for 
hydrogen atom in its ground state. Any other wave function chosen will 
always lead to an energy greater than the lowest exact energy of the 


system. 


7.3 Application of Perturbation Method to 


Helium 


The coordinates of electrons of Helium atom can be represented as 


Self-Instructional Material 


Unit 7 Jelectron 
HMO Method 


NOTE Mucleus 2 electron 
ra 


Hamiltonian operator of Helium atom is 
1 1 2 2 1 

H=- _ V°- V°- -+ inau 

2! 2? rnr 


where H(1; and Hoare the one electron hydrogen like 
Hamiltonians of the system. 


For the ground state of Helium atom the electronic 
configuration is 1s?. 


Slater determinantal form of overall wave function of 
Helium atom is 


1 md 
V o= Wis (1,2)=Ys (He) = glo 


- l-(s, is \a B -Ba ) 
JZ 2 9 O® 00 


This is the trial approximate wave function. According to 
variation method the energy based on the approximate wave 
function is always higher than the true energy Eo of the 
ground state. 


E - (y |H |v) 
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p ! Gsis aß -Ba ya 1 (isis (a p - zv 
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We yc Y ) Unit 7 


E= iiy HISS O. HMO Method 
1 
=" (15,15 ju (isis. )2 | NOTE 
, 0 oll e o 


= (USSA f5180))) 


Hamiltonian for Helium atom is 


1 
H= Ay + A+ " 


12 


1 
(15,150, ) 


E,s- (suis hr. + Hat — 
12 


(15:85 ) 


E\s= (sis fessis JH a sus Xasus 


Es-((s, [Ho s Js.) (Sian Fy Es KAS 15 y)+7 ) 


E s= Eisat Esot J 


Esa) is the energy of i^ electron in the 1S orbital of Helium and J is 


1 
called Coulomb integral. It is defined as J = isis. (sisa) and it 
r. 
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measures the Coulombic repulsion between electrons. 


Calculation of E, : 


In a two electron atom, the average potential experienced by an electron 
varies between 


Self-Instructional Material 


Unit 7 
HMO Method 


NOTE 


1 2 ea . 
- and-" that is it depends on the effective nuclear 

lj lj 
charge experienced by the i" electron in Helium atom. The 
best value of Z' , the effective nuclear charge will be 


determined by the variation principle. 


The wave function of 1S — orbital appropriate to this 
effective nuclear charge Z' is 


Vx 
This orbital wave function is an eigen function of the 
( 1 > a 
Hamiltonian operator! — mr 
r) 
2 
tpn te) 
with eigen value |- — 
2) 
E o, 
fav (eee (he 


T S2 C3). 


Eis(i), energy of the i^ electron in the 1S orbital of helium 
atom. 


Le Eis(i) = (15G)|H G) 186) 
l 
= AS()p- V?^-24S( 
( GA iu “| Oy 
Lae 2 y 
=O ue ede iso) 
= ASQ V- 
E 
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Eis) = E -(2-z)z 


Eis(1) = Eis(2). 


Since electron (1) and M are equivalent. 
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| 9E 
Minimizing the energy with respect to variable parameter < a E 
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NOTE 


0 
az A 8 J 
. 27 
2z —-__=0 
8 
NR 1.69 
16 


Substitution of this optimum value of z in energy 
expression we get, 


E —- (1.69) it gigo 
8 


= —2.8476au 
lau = 27.2eV 


E- -77.48eV 


This is the approximate ground state energy of Helium 
atom. Experimental value is -2.904au (-79 eV). Thus 
variation principle accounts for 98% of the observed 
energy. 


The effective nuclear charge experienced by an 
electron in Helium atom is 1.69 instead of the true charge 
2. The difference in these two values 0.31 represents the 
screening of the nucleus by other electron and 0.31 is called 
the screening constant or shielding constant. 


7.4 HMO Method 


Huckel adopted an approximation method to estimate 
energy and wave function of the molecular orbital formed 
by the combination of z molecular orbital. The theory 
proposes the following approximation. 


1. In a conjugated molecule the z electron does not 
interact with c electron. The complete wave function yw 
for the molecule is expressed as 
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Y= Vo.Vz 
Where y= complete wave function for o electrons 
V ; —complete wave function for z electrons. 
The energy E of the molecule is expressed as, 
E=E +E, 
2. Each HMO is considered as linear combination of 2P z atomic orbital 


Vi-aipit a 2P2+...... an Pn 


Vi- >) p. 
Where yi -i" HMO of the zelectron system 
pr=2P, atomic orbital of the r"carbon 
n = Total number of carbon atoms in the molecule 


a: = Coefficient of 2P, orbital of r" carbon atom in the i" HMO 


3. The energy of the i" HMO i.e. Eiis calculated by the formula, 


E E B yide 


y'dc 


i 


Application of variation principle with ‘n’ number of carbon 
atoms gives ‘n’ number of secular equation of the following types. 


ai(Hii-ESi1) + a2 (Hio-ES12) +...... an(Hin-ES in) = 0 
ai(H21-ES21) + a? (H22-ES22) +...... an(Han-ES2n) =0 


ai(Hni-ESn1) + a2 (Hn2-ESn2) a re an(Han-ESnn) =0 


3 h s A 
The corresponding n" order secular determinant is 
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NOTE 


Hy, m ES, H; -ES,, on Hi, -ES,, 
H z - ES, H » - ES, aS = H,, - ES, 

- = eee = =0 
Hy 7 ES, Ha 7 ES, Sj + = Ain - ES nn 


4. The columbic integrals Hpi represents the energy of the 
electron in 2P; orbital of i" carbon atom 


i.e. Hi = f pH p,dt 


where piis a wave function of 2P, orbital of i" carbon 
atom 
All such integrals are equal and denoted by symbol @ 
His = Ho...... 
Hm-=@ (say) 
5. The resonance integrals H:s represent the energy of 
interaction between r" carbon atom and s" carbon atom. 


le. Hts =[p, BH pdt 


where p. isa wave function of 2P, orbital of r^ carbon 
atom 
psis a wave function of 2P; orbital of s" carbon atom 


For non-bonded carbon atom Hrs = 0 


For bonded carbon atom all such integrals are equal and 
denoted by symbol 5 


Hı2 =- H23....... H34=....= B 
(say) 
6. Integral of type Sris equal to one. 
i.e. Su- fp, 'dr-1 
Sii = S22....... S33 Sinam sl 
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7. The overlap integrals Srs are assumed to be zero though it 
is not actually zero. This assumption simplifies the 
calculation. 


Ss= f p, p.d 


Si2=S823 CERNI =0 


Hence the secular determinant for a linear conjugated polyene is of the 
following form. 


Dividing each element by fand putting @-E/ B= x 


Then 
x 100-- 0 
1 x 10--0 
0 1 x - — 0 
000 1 -^- x 


The solution of the secular determinant is written as follows 
X = - 2cos (n z/N*1) 
«4 -E/ B= - 2cos (n z/N«1) 
Qa-E=-2 pcos (n z/N+1) 
Energy En = @ +2 pcos (nz /N+1) nz1,23.,...... N 
n = order of energy level (or) order of molecular orbital. 


N = Total number of carbon atoms in the molecule. 


Thus, energy can be written as 
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| En= @+2 f cos (nz 


Unit 7 /N+1) 
HMO Method 
NOTE 7.5 Application to Butadiene 


There are four 2pz atomic orbital to be combined and 
four z electrons in the molecule. The HMO’s are of the 
form, W= aipi+ a2p2+ aspa + a4p4. Where ai, a2, a3, a4 are 
coefficients. There are four coefficients, so 4x4 
determinants. 


The secular determinant of the system is 


x 100 

0 
1 X 1 =0 
0 1 x 1 
0.0 ] x 


x 1 0 1 1 0 
xil x 1-1] x 1=0 
0 1] x O 1 x 


x {[x(x*-1)-1(x)} -1( G0-1-100)) 20 
x (x*-x-x)-1(x”-1) 20 
xx x x +1 =0 
x^3 x*+1 20 
Put yz x y-3y4l-20 
y =-(-3)4 (3)? - 4ayay / 2 [^x = 
bx P? — 4ac /2a] 
= 3+ /9-4/2 
= 3445/2 


= 342.236/ 2 
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. y 20382 (or) y 22.618 


yes Unit 7 
HMO Method 
X sy Xe y 
x=+ 40.382 x = +./2.618 NOTE 
x = 40.618 ; x = -0.618 x= +1.618 ; x = -1.61 


Since @-E/ p= x 
E=a- Bx 
~ E: = @+1.618 2 
E2 = @ +0.618 5 
Es =& -0.618 5 
and = Es=@-1.618 8 
Thus the total energy is, E 7= 2( @+1.618 2) +2(@+0.618 B) -4a@ 
—4q +4.48 p -4a 
z Bond energy = 4.4855 
Delocalization energy = (4a 44.48 B) - 2x2(@+ B) 
—4q44.48 B-4a-4 p 
= 4.48 B-4 p 
= 0.48 5 
Wave functions: 


The wave functions corresponding to the energy levels E),E2,E3 
and E4 may be determined as follows. 


The secular determinant for butadiene is 


x1 0 0 
Ix 10 
01 x 1 
Ü 0. L x 
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| The secular equations are obtained as follows 


Unit? _ 
HMO Method axta 20 — —(LD 
aita2x +a3 = 0 —(2) 
NOTE 


a2*ta3x +a44=0 (3) 
àxtaax 20 —(4) 
For the energy level E1 x 2 -1.618 
Putting x = -1.618 in equation (1) we get, 
a1(-1.618) + a22 0 
7 a2- 1.618 ai 


Substituting the value of a» and x in equation 2 
a1*-1.618 ai(-1.618) +a3 = 0 
a1-2.618ai+a3 = 0 
a3 = -1lai+2.618a1 
a3=1.618a1 
Hence a2= a3 
Substituting the value of a3 and x in equation 4, we get 
1.618 ajt+a4(-1.618) = 0 
<. a1 = a4 
From the condition of normalization, 
altar” +a3*+a4" = 1 
n Wi= apı+1.618 a ıp2+1.618aı 
p3+a4p4 
The wave function yı is normalized as follows 


J far on 


yedt=1=(ap+16l8ap +1.618ap +ap)dt 
2 


[a 
| pdt (1.6182 | 747.6182)? | p/a a] gart=1 
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a2 [1+2.618+2.618+1] = 1 


7.236 4* 21 Unit 7 
HMO Method 
ai =1/ 47.236 
«a1 20.372 NOTE 
a2=1.618a1 
=1.618(0.372) 
a2=0.602 
a3 = 1.618aı =1.618(0.372) 
a3 =0.602 


Thus, 
y 1— 0.372pit 0.602p2+0.602p3+0.372p4 


Similarly we can find for y 2, w3 andy 4 by using the values x = - 
0.618, +0.618 and 1.618 respectively. The four HMO’s of butadiene are 
thus found to be as follows: 


V 1=0.372pi+ 0.602p2+0.602p3+0.372p4 
w 2= 0.602pit 0.372p2 -0.372 p3-0.602p4 

V 3=0.602p1-0.372p2-0.372p3+0.602p.4 

y 4— 0.372pi- 0.602p2+0.602p3-0.372p4 


7.6 Check Your Progress 


1. Explain the HMO method. 
2. Explain the HMO method and apply it to butadiene 


7.7 Answers to Check Your Progress Questions 


1. Huckel adopted an approximation method to estimate energy and 
wave function of the molecular orbital formed by the combination of z 
molecular orbital. 
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NOTE 


| 2 There are four 2p; atomic orbital to be combined and four xz 
electrons in the molecule. The HMO's are of the form, y= 
aipi + a2p2 + a3p3 + a4p4. Where ai, a2, a3, a4 are coefficients. 
There are four coefficients, so 4x4 determinants. Using the 
HMU metnod, we derive, bond energy = 4.48 p 


7.8 Summary 
The variation method to hydrogen atom gives, 


1-2 1 
E= —— au 


2 2 
This represents the exact value of energy and exact eigen function for 
hydrogen atom in its ground state. Any other wave function chosen 
will always lead to an energy greater than the lowest exact energy of 
the system. 


7.9 | Keywords 


e HMO: Huckel Molecular Orbital. 


° Variation method: For many problems, it is not practical to obtain a 
wave function by the exact solution of a wave equation that describes the system. 
It is still possible to perform many types of calculations, and one of the most 
useful techniques is that known as the variation method 


7.10 Self-Assessment Questions and Exercises 


1. Write down the steps involved in applying variation method. 
2. State and explain variation theorem with proof and apply it to 
hydrogen atom. 


7.11 Further readings 


Quantum Chemistry, I.N. Levine, Allyn and Bacon, 
Boston, 1983. 
Quantum Chemistry, R.K.Prasad, Wiley Eastern, New 
Delhi,1992. 
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8.0 Introduction 

8.1 Objectives 

8.2 Symmetry elements and symmetry operations 

8.3 Rules for forming a group 

8.4 Group multiplication table 

8.5 Group classification 

8.6 Point groups and systematic assignment of 
point groups for molecules 

8.7 Check your progress questions 

8.8 Answers to check your progress questions 

8.9 Summary 

8.10 Keywords 

8.11 Self-assessment questions and exercises 

8.12 Further readings 


point groups 


8.0 Introduction 


Group Theory is the mathematical applications of symmetry to an 
object to obtain knowledge of its properties. A point group describes all the 
symmetry operations that can be performed on a molecule that results in a 
conformation indistinguishable from the original. Point groups are used in 
Group Theory,the mathematical analysis of groups, to determine properties 
such as a molecule’s molecular orbital. In this section we are going to focus on 
the Symmetry elements and symmetry operations, rules for forming a group, 
group multiplication table, group classification, Point groups and systematic 
assignment of point groups for molecules. 


8.1 Objectives 


After going through this unit, you will be able to: 
e Understand about the Symmetry elements and symmetry operations. 
e Understand the rules for forming a group, group multiplication table. 
e Explain the concept of point groups and systematic assignment of point 
groups for molecules. 
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8.2 Symmetry elements and symmetry, 
operations 


SYMMETRY ELEMENT 


A symmetry element is a geometrical entity such as a line or 
a plane or a point about which an operation of rotation or 
reflection or inversion is done. 


SYMMETRY OPERATION 


A symmetry operation is a movement of the molecule such that 
the resulting configuration of the molecule is indistinguishable 
from the original. The molecule is taken into an equivalent 
configuration or and identical configuration. 


The symmetry elements and the corresponding 
symmetry operations are listed below: 


Symmetry element 


Symmetry operation 


Proper axis of symmetry 
(Cn) 


Rotation once or several 
times by an angle 0= 
(2z/n) about the axis. 


Plane of symmetry (o) 


One or more reflections 
in the plane. 


Improper axis of 
symmetry (Sn) 


Rotation about the axis 
followed by reflection in 
a plane perpendicular to 
the rotation axis. 


Centre of symmetry (i) 


Inversion of all atoms 
through the centre of 
symmetry. 


Identity element (E) 


This operation leaves the 
molecule unchanged. 


PROPER AXIS OF SYMMETRY (Cn) 


This axis of symmetry can be explained by taking the example 
of triangular planar boron trichloride molecule. In boron 
trichloride molecule an axis of symmetry is located 
perpendicular to the plane containing all the atoms. This is 
known as the C3 axis of symmetry. In general the symbol for 
proper axis of symmetry is Cn, where n is known as the order of 
the axis. 


The order of the axis is given by the number of rotations by 
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0, to get the identical configuration. n is alternative given by 
the formula. 


n = (27/0) 


where 0 is the minimum angle of rotation to obtain the equivalent 
configuration. 


(b) 


PLANE OF SYMMETRY (o) 


The plane containing all the atoms is called as molecular plane. 
[PtCL]" ion contains a molecular plane and four more reflection planes. 
Water has a reflection plane passing through the oxygen atom and another 
one containing all the atoms. 


The reflections plane in water molecule 


The reflection plane is of three types:- 


1) A plane is referred to as horizontal plane (on) if it is 
perpendicular to the principal axis. 
ii) A reflection plane which contains the principal axis is called as 


vertical plane (oy). 
iii) A vertical plane which bisects two perpendicular C» axes is 
called a dihedral plane (oa) e.g. allene. 
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IMPROPER AXIS OF SYMMETRY (Sn) 


Unit-8: 

Assignment It is the line about which a rotation by a specific angle 
of point followed by refletion in a plane perpendicular to the rotation axis 
groups is performed. E.g. Seaxis in staggered form of ethane. 

MP a nai Tall 
NOTE m j u e M 
Ce —— 
- G es : 
br x cc C 
CENTRE OF SYMMETRY (i) 

This is the point such that any line drawn through it meets the 
same atom at equal distances in opposite directions. All 
homonuclear diatomic molecules posses the centre of symmetry. 

x Jf ^ 
j | opaan Pa |. 
S 
"4 T b: 
8.3 Rules for forming a group 
The symmetry elements of a molecule 
must satisfy certain rules in order that they form a group. The 
rules are 
1. The product of any two elements or square of each element 
must be an element of that group. 
2. In each group, there is one element which commutes with every 
other element and leaves it unchanged. In general the identity 
element does so. 
3. The associative law of multiplication hold good. 
4. Each element has an inverse or reciprocal which is also element 
of that group. 
If the element is A and its inverse is A” 
Then, A A” =E 
Self-Instructional For e.g. consider the set of numbers 
Material 
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(1/4, 1/2, 1, 2, 4, 8, 16] 
By rule (1) 
2x4=8 
ded 
4 = 16 
By rule (2) 
A'B-BA 
1x2 =2x1l =2 
1x4 =4x1 =4 
By rule (3) 
A (BC) = (AB) C 
2(4x8) = (2x4)8 
By rule (4) 
A A! 
xe] 
Y4x4=] 


8.4 Group multiplication table 


It consists of rows and columns. Each is represented by a group element. For 
multiplication we can use the following method. i.e. Column element x Row 
element. 


For constructing a group, rearrangement theorem must be known. It 
states, 


“Tn each column or a row a particular elements occur only one". 
For example, 


Consider a group of order 2. This group consists of two elements only. 
Let the elements be E and A. The first step leads to the following table. 


G| E A 
E IE A 
A IA X 
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According to rearrangement theorem ‘X’ should be ‘E’ only and 
cannot be ‘A’. ..the group multiplication table of G2 is 


Example -2 


E Co Ov Ov 
C2 E Ov Ov 
Ov Ov E C2 
Ov Ov C2 E 


The following features of the multiplication table carry 
significance: 


1) E combining with any operation X gives X as the result 
since E does “nothing” i.e., EX=XE =X. 

2) C2 performed on any object, interchanges both left and 
right and front and back. 

3) o. interchanges only left and right while ov, 
interchanges only front and back. 

Thus oy oy’ = oy ovis the same as front-back +left, -right 
interchange 


i.e. C2C20v= Oy Oy o, = OVE = oy and similarly C2 
ov = Ov 


8.5 Group classification 


FINITE GROUP 


A group that contains definite number of symmetry 
elements is called the finite group. For example, water molecule 
has 4 elements of symmetry and Ammonia molecule has 6 
elements of symmetry. Therefore they belong to the finite group. 


INFINITE GROUP 


A group that contains infinite number of symmetry 
elements these are called as infinite group. 


For example: linear molecules H2 (H-H), Br? (B.-B;), HCl 
and CO» can be rotated along the molecular axis to any degree. 
Therefore they have infinite axes of symmetry. Such molecules 
belongs to infinite group. 


ABELIAN GROUP 
Abelian group is a group in which each element commutes 
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with every other elements. 


Unit-8: 
For eg:- In water molecule, each element commutes with every other Assignment of 
element and therefore water molecule belongs to abelian group. Such a point groups 
situation dos not arise in NH3 molecule and hence it is non-abelian. 
CYCLIC GROUP 
A group in which one element generates all the other elements of the NOTES 


group is known as cyclic group. 


For e.g. the elements of symmetry that are present in H20» are E and Co. 
The element C» generates the other element E since C2? = E. 


SUB GROUP 


A sub group is the smaller group of a main group in which the 
elements of symmetry must satisfy the requirements for forming a group. 


The order of a subgroup is given by 'g'. The order of a maingroup is 
the integral multiple of the order of the subgroup. 


h-gxn 
Where ‘n’ is the integer. 
In water molecule, there are four subgroups. They are 
(E) (E, C2) (E, oy) and (E, ov) 


In general, the identity element in any group will form one of the subgroup. 


8.6 Point groups and systematic assignment of point groups 
for molecules 
POINT GROUPS 


Depending upon the symmetry elements molecules are classified in to 
different groups. Molecular groups are called point groups. Because the 
symmetry elements intersect at a point and this point does not move during 
symmetry operations. Molecules are classified into different point groups 
based on certain combinations of symmetry operatins. The system of 
notation used for point group is known as schoenflies system. In this case 
symbols like Cnv , Doh , On and Ta are used. In order to identify the point 
group, the following flow chart will be helpful. 
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Unit-8: 
Assignment 
of point 
groups 


NOTE 


Self-Instructional 
Material 


Is the molecule linear? 


* 
Yes 


* 
Does the molecule 


s 4 
No 


Does the molecule have many C, axes of order 


have 1? greater than i 
4 4 Y + 
Yes No Yes No 
Dor Cay + + 
The molecule Does the molecule have 
belongs to the n C, axes perpendicular 
cubic group to the C, aae 
i i + n Y 
ABs ABs AisBis Yes No 
(tetra- (octa- (icosahedral) 4 
hedral) hedral) | The group is 
4 4 x dihedral 
T, On I. — EX 
+ + + 
+o, «no, no planes 
+ 4 4 
D. D, 2 D. 
$ - 
Does the molecule have the 
C, axis? 
* + 
Yes No 
+ + 
+ + + Does the molecule have the 
+o, +noy no planes reflection plane alone? 
Y + + 
Car Cav Ca + Y 
Yes No 
+ 4 
€, Does the molccule have 


only the centre of symmetry? 


* + 


Yes No 
+ + 
Ci Does it have only the S. axis 
of even order greater than 2? 
Y 
Yes No 
4 
S. Does the molecule have 
only the identity 
element? 
+ 
Yes 
+ 
Cı 


8.7 Check Your Progress 


1. Distinguish symmetry elements and symmetry operations? 


2. Explain Centre of symmetry? 
3. Define Cyclic group? 
4. Define Point group? 
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8.8 Answers To Check Your Progress Questions Unit-8: 


l. Assignment of 
Symmetry element Symmetry operation point groups 
Proper axis of symmetry (Cn) Rotation once or several times by an 
angle 0= (2z/n) about the axis. 
Plane of symmetry (o) One or more reflections in the plane. 
Improper axis of symmetry (Sn) Rotation about the axis followed by NOTED 
reflection in a plane perpendicular to 
the rotation axis. 
Centre of symmetry (i) Inversion of all atoms through the 
centre of symmetry. 
Identity element (E) This operation leaves the molecule 
unchanged. 


2. This is the point such that any line drawn through it meets the same atom at 
equal distances in opposite directions. All homonuclear diatomic molecules 
posses the centre of symmetry. 


3. A group in which one element generates all the other elements of the group 
is known as cyclic group. e.g. the elements of symmetry that are present in 
H20: are E and C». The element C» generates the other element E since C2” = 
E. 


4. Depending upon the symmetry elements molecules are classified in to 
different groups. Molecular groups are called point groups. 


89 Summary 


e A point group describes all the symmetry operations that can be 
performed on a molecule that results in a indistinguishable from the 
original. 

The different types of groups and their definition are symmerized. 
Symmety elements and their operation such as Proper axis of 
symmetry (Cy), Plane of symmetry (o), Improper axis of symmetry 
(Sn) and Centre of symmetry (1) are clearly explained. 


8.10 Keywords 


Point group: point group describes all the symmetry operations that can be 
performed on a molecule that results in a conformation indistinguishable from 
the original. 

Group classification: A subset of a group that is closed under the group 
operation and in which every element has an inverse in the subset. 

Group multiplication table: Every row contains each element exactly 

once and every column contains each element exactly one. 
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8.11 Self-assessment questions and exercises 


Unit-8: 

A w 1. Explain in detail about the rules for forming a group 
ssignment : : : f 
insta 2. Discuss in detail about symmetry group and operations 

r 3. Dicuss the point group 
groups A dU 
4. Explain Group classification 
NOTE 


8.12 Further readings 


3. Chemical Application of Group Theory, F.A. 
Cotton, John Wiley and Sons Inc. New York,1971. 
4. Group theory and its applications to Chemistry, 


K.V. Raman, Tata McGraw-Hill Publishing Company,1990. 


Self-Instructional 
Material 
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Unit- 9 Matrix representation theory Unit -9 


Structure Matrix 
9.0 Introduction Representation 
9.1 Objectives 

9.2 Matrix representation theory 

9.3 Matrix multiplication 

9.4 Inverse of a matrix 

9.5 Matrix diagonalization 

9.6 matrix representation for symmetry operations 

9.7 Check Your Progress 

9.8 Answers to check your progress questions 

9.9 Summary 

9.10 Keywords 

9.11  Self-assessment questions and exercises 

9.12 Further readings 


9.0 Introduction 

The symmetry operations in a group may be represented by a set of 
transformation matrices, one for each symmetry element g. Each 
individual matrix is called a representative of the corresponding 
symmetry operation, and the complete set of matrices is called a matrix 
representation of the group. In this context, we discuss in detail about 
matrix representation theory,matrix multiplication, inverse of a matrix, 
matrixdiagonalization and matrix representation for symmetry operations 


9.1 Objectives 
After going through this unit, you will be able to: 


e Understand the concept Matrix representation theory 
e Explain the Matrix multiplication and Inverse of matrix. 
° Understand the Matrix diagonalization and Matrix representation 


for symmetry operations. 


9.2 Matrix representation theory 

During a symmetry operation, a set of 
coordinates of an atom in a molecule is transformed into a new set of 
coordinates. These two sets are related to one another in the form of a set 
of equations. This set of equation can be formulated in the form of a 
matrix. Each symmetry operation is represented by a matrix. Matrix 
representation of the symmetry operation is useful in the study of 
structural problems. 


Self-Instructional Material 


Unit -9 9.3 Matrix representation for symmetry operations 


Matrix The matrix representation for each kind of symmetry 
Representation operation can be obtained as follows. 
Cn OPERATION 
NOTE 


Consider the figure shown, 


The position of an atom in a molecule is represented by 
the vector OP whose magnitude is ‘r’. Point ‘P’ corresponds 
to the coordinates x,y,z. If the molecule is 


Rotated in the clockways direction through an angle 0 the 
atom changes its position from the point ‘P’ to the point ‘R’. 
point ‘R’ has a new set of coordinates namely x2y222. 


Consider the APOQ, it follows from the knowledge of 
trigonometry that, 


x12 OQ =r Cos Q 
yı =PQ= rsin c 
Consider the AROS, 
Again, 
x2= OS =r cos B =r cos (a - 0) 
y2=RS =r sin B =r sin (a - 0) 


Self-Instructional z2 = Z1 (since the rotation is carried out about the z-axis). 
Material 
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-.X22rcosQ cos ĝ + r sin & cos 0 = xı cos 8+ yı sin 0 +0 


y2=r sin æQ cos 0 - r cos Q sin 0 = yı cos O + xı sin 0 +0 
=040+4+ zi. 
Z2 = Z1 
cos(a tp) = cos a cosp + sin & sin B 


sin (a+) = sin & cosp + cos ac sin D 


This set of equations can be formulated in the form of a matrix. 


|| cosé sinô O] x | d 
m cos Q 0|» =| y| 


| 0 0 alaj ial 


Therefore the matrix representation for Cn axis is given as, 


| cos sind 0| 
C —|-sinOÓ cos8 " 
| 0 0 1 


If the rotation is carried out in the anticlock direction, the matrix will be, 


[cosO —sinQ 0| 
C. 2 and cos@ 0| 
| 0 0 1 


For e.g. the matrix representation for a C» axis is given by 


|| cos180 sin 180 0| j-1 0 0| 
C= -sin 180 cosl80 0-2 0 


-1 0 
2 | 
| 0 o 1| lo o ıl 
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Unit -9 


Matrix 
Representation [ 120 120 qT- 1 Be 0l 
c - Sin i20 mum $ i 0 
NOTE "| i e EI 


For an identity operation, the new set of coordinates will be 
the same as that of the original. 


cos (904-0) = sin 0 
cos (27040) = sin 0 
sin (904-0) = cos 0 
sin (27040) = cos 0 


Therefore, 
xX2=xi+0+0 
y2=0+y +0 


z2=0+0+z2 


1 0 O| x X2 
a 2 TT La 


he eles 


o o ilal liz l] 


-. The matrix representation for the identity operation is 


given as, 
| Oo O0 
E-|O 1 0 
lo o 1| 


INVERSION OPERATION 


Self-Instructional 
Material 
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X2=-Xit0+0 Unit -9 


Matrix 
y2=0-y2+0 Representation 
22=0+4+0-2z2 
|-1 0 olx] | 2 NOTES 
0 -1 
| 0|»[7]X| 


Ilo o -ilal lal 


-. The matrix representation for the identity operation 


[-1 0 d 
UTE 


REFLECTION OPERATION 


During a reflection operation, the coordinates corresponding to the 
reflection plane do not change the sign while that of the other coordinates 
change their sign. Therefore the matrices corresponding to the reflection 
operation is given below, 


[o 0 -1| 
[-1 0 0 

o,=| 0 1 i 
T 0 1| 
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[1.0 0| 


Unit -9 
Matrix 0,7/0 -1 0 
Representation 0 0 i 
NOTE 
Sn OPERATION 
Sn operation involves rotation of the molecule about an 
axis through an angle followed by reflection in a plane 
perpendicular to the axis. .. the matrix for a Sn operation is 
obtained by multiplying the matrices corresponding to Cn and 
Oxy provided z-axis is taken as the rotational axis. 
Scc s" 
|| cosé sin Ó ofi 0 0 | | cosÜ  sinÓ 0| 
S= -sin cos@ n 1 0 poe cos 0 
n | 
| o o ılo o -1| [| o o -1| 
S5 2 Ca.” = Qc 
9.4 Inverse of matrix 
A square matrix B is called the inverse of A if 
AB =BA=1 
Where Iis the unit matrix. Inverse does not exist if the matrix 
is singular. The inverse of the matrix A can be obtained by 
using the following rules: 
|; Compute the determinant of the matrix A. 
A Interchange the elements ai: and a22. 
3. Change the signs of ai? and a21. 
4. Divide each element of the matrix thus 
formed by the determinant of A. 
Self-Instructional A- B i d 
Material 21 22 
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Unit -9 


The resulting matrix is the inverse of A and is Bs s 
represented as A’. The inverse of the matrix B 3 
3 — 
pul? 72 
E 4 
NOTES 


is obtained as follows. The determinant of the matrix B is equal to 14. 
Interchange of the elements 3 and 4 leads to the matrix. 


When the signs of -2 and 1 are changed, the following matrix results: 


[4 2] 


l-1 3] 


Each element of this matrix is divided by 14 to obtain the inverse of B.B”! 
thus obtained is given below: 


9.5 Diagonalisation of a matrix 


The process of reducing a matrix to the diagonal matrix is referred to 
as diagonalisatoin. Let A be the square matrix of order n. P is the 
similarity transformation matrix which reduces A to the diagonal matrix D 
according to the equation. 


p'AP-D 
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Let us consider the matrix A 
Unit -9 
Matrix 
Representation 


NOTE 


And its inverse P`! 


can be used to reduce A to the diagonal matrix according to 
the equation 


e alpes) elm on 


I-3 2 | |-is TIME 2| 


P'AP = 


Self-Instructional 
Material 
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Unit -9 
Matrix 
Representation 


.6 Matrix representation for symmetry operations 
a PETER : y y op NOTES 


Using carthesian coordinates (x,y,z) or some position vector, we 
are able to define an initial position of a point or an atom. 


The initial vector is submitted to a symmetry operation and thereby 
transformed into some resulting vector defined by the coordinates 
x', y and Z'. In an algebraic context, this transformation is expressed 
a matrix which processes the initial position vector. We write 


final vector = Matrix * initial vector. 


The most primitive symmetry operation is the identity and yields a final 
vector identical to the initial vector. It is the unity matrix or identity 
matrix which leaves all coordiates unaffected. 


0 0 | 


|] 
E = |o 1 0 
lo 0 ıl 


If we want to perform a reflection on the xy-plane (analogous to a 
horizontal plane on), coordinate z changes the sign. 


It sO — 5: 
-= | | Ü 
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The matrices which are applied for performing a reflection 
on the yz-plane and xz-plane are the matrices ox and 
oy respectively. 


I-1 0 0| i1 QO O0 
g, = | 0 1 d g,- fo -1 0 
lO 0 0 0 0 ol 


The inversion i relates the coordinates (x,y,z) with (-x,-y,-z) 
and is connected with the following matrix: 


Obviously, a twofold application of the inversion matrix 
yields the coordinates of the initial point (x,y,z) which is 
reflected by E = i*i. 


Li o ola o ol) Ht ool 
TE | 0 -1 0 | 0 -1 Ol= |o l o| =f 
lo o =o o -1| loo 1 


The matrix for a rotation about axis z by an arbitrary angle © 
is derived easily if we imagine two two-dimensional 
coordinate planes with identical origin but an angular 
difference of O between the axes. In our context of 
symmetry, we just need to deal with the discrete values of © 
= 2n/n for the angle of rotation. 


| cos2T/n sin2rr/n 0! 
C, = |-sin2r/n cos2m/n O 
| 0 0 1| 


The matrices for the symmetry operations C2(z), C3(z), Ca(z), 
Cs(z) and Co(z) are obtained easily. The matrices for C," as 
symmetry operation are calculated by an n-fold multiplication 
of matrix Cn. The symmetry operation C2 around axis x 
(x—X,y—-y, z—-2Z) and around axis y are (x-x, y—y, z—- 
z): 
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As we know rotatory-reflection to be a combination of rotation and 
reflection, a matrix representation for this operation is easily to be 
derived.For instance, to obtain the matrix for rotatory reflection Sn(z) we 
multiply the matrices for the fundamental operations &sigmaz and Cn. 


l1 0 Ollcos2m/n sin2m/n 0| |cos2m/n  sin2m/n 0 | 
Sz) = 9C, =|0 1 0 | —sin2 m/n cos2m/n o| = |-sin2m/n cos2m/n 0 
10 0 -1/l 0 0 1} | 0 0 =F] 
9.7 Check Your Progress 


1.What is matrix representation? 
2.Givethe rules of Inverse of matrix? 
3.what is the equation of Sn OPERATION 


4.what are the types of symmetry operations in matrix representation? 


9.8 Answers to check your progress questions 


1. The symmetry operations in a group may be represented by a set of 
transformation matrices, one for each symmetry element g. Each 
individual matrix is called a representative of the corresponding symmetry 
operation, and the complete set of matrices is called a matrix 
representation of the group. 


2s 

Compute the determinant of the matrix A. 
Interchange the elements ai: and a22. 
Change the signs of ai? and a21. 


Divide each element of the matrix thus formed by the determinant 
of A. 


149 


Unit -9 
Matrix 
Representation 


NOTES 


Self-Instructional Material 


Unit -9 
Matrix 
Representation 


NOTE 
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Material 


Sn= C0” 


|| cose sin Q ofi 0 0 | | cos0 sin 0 0 | 
S='-sin Ó cos@ 0° 1 na cos@ P 


| 0 o ılo o -1] | 0 0 -1| 


S3° = Cie” = Go 


Cn operation 
S noperation 
Inverse operation 
Reverse operation 


VVVVF 


9.9 Summary 


° A matrix representation theory describes all the 
symmetry operations that can be performed on a 
mathematical experision 

e The process of reducing a matrix to the 
diagonal matrix is referred to as diagonalisatoin. 

° Matrix representation of the symmetry 
operation is useful in the study of structural problems. 


9.10 Keywords 


Matrix representation-A set matrix of symmetry operations. 


Inverse of matrix -AB = BA = 1 


Operations - Cn S n,Inverse Reverseoperation 


9.11 Self-assessment questions and exercises 


1. Explain in detail about Matrix representation theory and 
matrix multiplication 

2. Give a brief review Inverse of a matrix 

3. Explain matrix representation for symmetry operations 
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9.12 Further readings 


1. Chemical Application of Group Theory, F.A. 
Cotton, John Wiley and Sons Inc. New York, 1971. 
2. Group theory and its applications to 


Chemistry, K.V. Raman, Tata McGraw-Hill Publishing 
Company,1990. 


151 


Unit -9 
Matrix 
Representation 


NOTES 


Self-Instructional Material 


Unit -10 
Character table 


NOTE 


Self-Instructional 
Material 


UNIT- 10 Character table 


Structure 
10.0 Introduction 
10.1 Objectives 


10.2 Reducible and irreducible representations 
10.3 Character table of C2v and C3v point groups. 
10.4 The great orthogonality theorem and character 


10.5 Direct product representation 

10.6 | Check Your Progress 

10.7 Answers to check your progress questions 
10.8 | Summary 

10.9 Keywords 

10.10  Self-assessment questions and exercises 
10.11 Further readings 


10.0 Introduction 


Representation is a set of matrices which represent the operations 
of a point group. It can be classified in to two types, which are 
reducible representations and irreducible representations. A 
character table is a 2 dimensional chart associated with a point 
group that contains the irreducible representations of each point 
group along with their corresponding matrix characters. It also 
contains the Mulliken symbols used to describe the dimensions of 
the irreducible representations, and the functions for symmetry 
symbols for the Cartesian coordinates as well as rotations about the 
Cartesian coordinates. In this context, we discuss in detail about 
reducible and irreducible representations,Character table of C2v 
and C3v point groups,the great orthogonality theorem and character 
and direct product representation. 


10.1 Objectives 


After going through this unit, you will be able to: 


° Explain about the reducible and irreducible representations 
. Understand about the character table of C2v and C3v point 
groups 

° Understand the concept of the great orthogonality theorem 
and character 

. Learn about the direct product representation. 
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10.2 Reducible and irreducible representations 


Unit-10 
Character table 


REDUCIBLE REPRESENTATION 


Reducible representation and its reduction can be undearstood by 
carrying out a similarity transformation. NOTES 


Let A,B,C,D and P be the similarity transformation matrix in this 
group.By similarity transformation the matrices A,B,C,D and P are changed 
into A’,B’,C’,D’ and P' as 

P!AP-A' 
P’ BP =B 
P'CP=C 
P! DP =D’ 
P’ Ppp=P’=P 


If the resulting matrices can be blocked into smaller matrices, then the 
representation T is called a reducible representation. For e.g. A’ can be blocked 


into ay’, a2’, a3, a4 ,......... ae, as sub-matrices as 
la 0 0 0 0 0| 
0 a 000 0 
we 0 9g da 0 0 0 
0 0 0 a 0 0 
0 0 0 a! 0 
o 9 0 0 0 aj 


A reducible representation result when the various symmetry operations are 
performed on all the sigma-bonds of a molecule for e.g. In BF3 molecule (Dsh 
point group) the 12 symmetry operations of this group are divided into 6 
classes as 
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(E) (C3lCO35, (3C25,80),(53,S3) and 


| 
Unit-10 (On). 


Character table 


The three sigma-bonds in BF3 are considered as vectors 

r1, rz and rs. By a symmetry operation R, the vectors ri , r2 and r3 

NOTE ri, 2 and r3 are changed into ri^ ,r?' and r3' according to matrix 
equation. 


f pn 
VIL mem 


In | E 


Where [R] is the matrix for the operation R. The equations 
relating rı’ , r? and r3' and ri, r2 andr3 are 


r! 21r Or +0r 
1 1 2 3 


r! 2 0r-4 lr Or 
2 1 2 3 


r'=0r+0r +1r 
3 1 2 3 


In matrix form the above equation becomes 


"ENT: 
0 1 Mara met 


T = 

2 | 

n lo o tl) [7 

From equations (a) and (b) we get 
100 
0 1 0 
[0 o 1| 


As the matrix for identity operation. 


For C3! operation the resulting vectors ri' , rz and r3’ are related to 


ri, r2 and r3 by the equation. 
Self-Instructional 
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r! 2 0r-4 1r Or 
1 1 2 3 


=e Or TUE 


5j — Or t Or 1r, 
pez 
2 2 


Unit-10 
Character table 


These equations can be represented in the matrix form as 


P 0 1 r, 
ORT 0 Me 
E li 0 ol i| 


From (a) and (c), we get, 


alo 0 r! 
li 0 Ol 


Similarly we get reducible representation for 


S3! 


REPE. 


lo 1 ojo o 1f1 o ojo 


Ov Oh 


NOTES 


on 
0 0 
B 
1 Ol 


The characters of the matrices can be found by using the rule 
"The character of the matrix for a symmetry operation is equal to the 


number of unshifted vectors by the operation". 


The number of vectors unshifted by the identity operation is three. 


k 0 j| 

E=|0 1 0 
i 0 1| 
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The character for this matrix is also 3. 


By using the rule the characters for the different matrices 
of operations of BF3 are obtained. 


Symmetry E | Cs |C2 | on | ov | Ss 
operations 
Character of | 3 0 1 3 1 0 
the matrix 


IRREDUCIBLE REPRESENTATION 


If it is not possible to find a similarity transformation 
matrix which will reduce the matrices of representation I’ then the 
representation is said to be irreducible. All one-dimensional 
representations are always irreducible. 


Example: 


Considering the matrices of transformation for the Z-coordinate 
of a hydrogen atom in hydrogen molecule, which result by the 
symmetry operations of Dh group. The operations of D«h group 
are E, Cæ, ov, C2, Oh, Se and i. From the figure the Z-coordinate is 
unaffected by E, C» and ov operations. 


The Cx axis, vertical plane and horizontal plane in H» molecule 
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The equations and matrices for the transformation of Z-coordinate of 
Character table 


hydrogen atom by these operations are 


E.Zz21Z 


NOTES 
E matrix = [1] 
Co. Z = 1Z Co. 
matrix = [1] 
ov.Z=1Z 


o, matrix = [1] 


All the other operations of this group change the coordinate Z of 
hydrogen atom into —Z, we get 


cC .Z =—1z 
S..-Z=—1z 
O, ox = —1z 
i.z ——1z 


C» matrix =[-1] 
s, matrix = [-1] 
O ,, matrix = [-1] 


i matrix = [+1] 


The matrix representation thus obtained for the z-coordinate of the 
hydrogen atom in the hydrogen molecule is given below. 


E Co Ov C2 So Oh i 
F qug p] pug. BÉ Gt) 
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This representation T is irreducible since it is one- 
dimensional. 


Reducible representations and irreducible representation 
play an important role in obtaining solutions to problems of 
hybridization, molecular vibration, delocalization energies of II 
electron system and so on. In all this applications the first step 
involves point group determination and formation of the 
reducible representations of the group. Every point group 
consists of a certain number of irreducible representations. 


The characters of matrices in the different irreducible 
representations of a point group can be listed in a table known as 
character table. The character table for a group can be 
constructed with the knowledge of properties of irreducible 
representation. The properties of irreducible representation can 
be obtained from the great orthogonality theorem. 


PROPERTIES OF IRREDUCIBLE REPRESENTATIONS 


A knowledge of the properties of irreducible representations is 
essential to construct the character table of a point group. 


Let us consider a point group consisting of h symmetry 
operations. These operations are divided into K classes. The 


13: lk are dimensions of these representations. The 
orthogonalilty theorem is used to obtain the following rules for 
the irreducible representations. 


l. The number of irreducible representations in a group is 
equal to the number of classes in the group. 
2: The sum of the squares of the dimensions of the 


irreducible representations of a group is equal to the order of the 
group. 


K 
NET 
i-l 
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In a matrix of order | there are square elements thus each 

irreducible representation I; will provide 1! 2 = h 

dimensional vector. The basic theorem requires this set of 

lebt sb vectors to be mutually orthogonal. NOTES 

Since there can be more than h dimensional, h orthogonal 

vectors the sum of l1? +p? & 13 ......... may not exceed 

h. 


3. The sum of the squares of the characters of identity operation in 
the irreducible representation is equal to the order of the group h. 
k 


2 [rao] -h 


4. The vectors whose components are the characters of two different 
irreducible representation are orthogonal. 


S Oo; Q(Q =0 ; Where ij 


5. In a given representation the characters of all matrices belonging 
to operation in the same class are identical. 


Symbols used for representing irreducible representation: 


1. Bethe's method is used to label irreducible representation such as 
I1,I?,I35. 

2. Mullikan s method is based on symmetry of irreducible 
representation. 

3. Symbol A or B represents one dimensional representation. E 
represents two dimensional representations and I represents three 
dimensional representation. 

4. If I'2 +1 with respect to principle axis, the symbol A is used. If I 
= -] the symbol B is used. If there are several irreducible 
representation subscript 1, 2, 3 etc are used. 

Example: 
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If there is only one representation the subscript one 
is not necessary. 


5. If the operation is inversion and I' = +1 the 
subscript g is used (g=gerade). If the operation is inversion 
and T = -1 the subscript v is used (ungerade). 

6. If I = 1 for all the symmetric operation of a point group 
Ai or Ag is used which is called totally symmetric 
irreducible representation. 


T. Single prime denotes that the irreducible representation 
is symmetric with respect to on. 
8. Double prime denotes, that the irreducible representation 


is antisymmetric with respect to on. 


10.3 Character table of C2v and C3v point groups 


CHARACTER TABLE 


The characters of matrices in the different irreducible 
representations of a point group can be listed in a table known as 
character table. 


CONSTRUCTION OF CHARACTER TABLE FOR C?V 
POINT GROUP 


C2V point group consists of 4 elements such as E, C2, ovx; and 
Ovyz, each is in a separate class. Hence, as per rule 


(1) ie. the number of irreducible representations of a group is 
equal to the number of classes of the group. 


CoV E C Ov(xz) Ov(yz) 
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There are 4 irreducible representations of this group but it is also required 
as per rule. 


(2) *The sum of the squares of the dimensions of the irreducible 
representation of a group is equal to the order of the group h". Thus we 
are looking for a set of 4 ‘+ ve’, integers , li, l2, 13 and 14 which satisfy 
the relation, 117+ 13+ 134+ L2. = h. clearly the only solution is l= l= l = 
l4 =1. Thus the group C2V have 4 one dimensional irreducible 
representations. 


On the basis of the vector properties of the representations and the rules 
derived above, one suitable vector in four space which has a component 
of one corresponding to I; will obviously be 1,1,1,1 for 


Des =h 


1? +1? +1? +1? =4 
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| CoV E C2 Ov(xz) Ov(yz) 


Unit-10 
Character table Ii 1 
I? 1 
NOTE I5 1 
I4 1 


Thus satisfying rule 


(3) i.e. the sum of the squares of their characters in any 
irreducible representation equals h. 


C2V Co Ov(xz) Ov(yz) 


I^ 1 
I5 1 
I4 1 


Now all other representations will have to be such that 
X[X(R)] -h-4 which can be true only if each Xi(R)- +1. 
Moreover in order for each of the other representations to be 
orthogonal to I'1. According to the rule 


(4) i.e. the vectors whose components are the characters of two 
different irreducible representation are orthogonal. 


ie. X Xi(R) Xi (R)=O when iz. 
R 
There will have to be 2 “+1’S” and 2”-1’S”. 
(1) Q (0) CD (0) CD) «(0 (0) 20 
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Material 
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Therefore we can write the other irreducible representation as follows. 


Unit-10 
Character table 
CV |E C» owxz) Ovoz) 
Ii 1 1 1 1 
I^ 1 -l -1 1 
T3 1 -l 1 E NOTES 
T4 1 1 -1 -1 


CONSTRUCTION OF CHARACTER TABLE FOR C3V POINT GROUP 


C3V point group consists of 6 symmetry elements. E, C3', Cs’, ov, ov' 
and ov". It can be listed by classes as E 2C3, 3ov. Rule (1) is the number of 
irreducible representation = number of classes. 


CGV JE 2C3 3ov 
I 
I? 
I5 


We therefore known that there are 3 irreducible representations, by rule (1), if 
we denote the dimensions by li, l2 and 13 we have li?» 124 13? = h= (6). The 
only values of the l; which will satisfy this requirement are 1,1and 2. 


CV IE 2C. 3ov 


Thus as by rule (2), the sum of the squares of the dimensions of the irreducible 
representation of a group is equal to h. 
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1? 41? +27 = 14144 2 6 


Always in any group there will be a one dimensional 
representation whose character is equal to ve one. Therefore 
one dimensional irreducible representation in C3V point group is. 


C3V E 2C3 30v 


Thus by rule (3) 


The sum of the squares of the characters of identity operation in 
the irreducible representation is equal to the order of the group h. 


LEXi(R)]? zh = -101 + 200 + 3(1)? = 6 


As per rule (4), the vectors whose components are the characters 
of two different irreducible representations are orthogonal. 


We now look for a 2" vector in space whose components are 
equal to + or -lwhich is orthogonal to Tı . The components of 
such a vector must consists of 3’+1’S’ and 3’-1’S’. Since X(E) 
must always be ‘+ve’ and since all elements in the same class 
must have representations with the same character, the only 
possibility here is 
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CGV| E 2C3 3ov 
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Character table 


I3 2 NOTES 


Now our 3"! representation will be of two dimension i.e.X3(E) = 2. In 
order to find out the values of X3(C3) and X3(ov), we make use of orthogonality 
relationship. 


2 00250 = 0 
DAR) ZB) = (DD) + DOC) + G)X00059y ) = 0 
=2+273C,+ 3730) =0 
DAR) ZAR) = (D) + QD %BC3) + B(-DGov ) = 0 


—2--2045C;,-3250y-0 


245C, -3550y =-2 >(1) 
2:5C, -3% =-2 >(2) 
473C, + =-4 
gc-75-a 

3 3 —4 


Substitute /5C,- -1 in equation (1) 
2(-1) +3 50, =-2 
-2+3 730, =-2 


350,- -2+2 
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3450,20 


A50 ,-— 0 


Thus the complete set of characters of the irreducible 
representation is 


C3M E 2C3 3o0v 


10.4 The great orthogonality theorem and character 


This theorem is concerned with the elements of matrices 
constituting the irreducible representation of a point group. Let us 
consider two irreducible representations i and j of a point group. 
Let l; and 1; be the dimensions of these representations, h is the 
order of the point group. R denotes a particular symmetry 
operation in the group. (T; ( R))mn is an element in the m" row 
and n" column of a matrix in the i" irreducible representation. 
The complex conjugate of the element in the m™ row and n™ 
column of a matrix in the j" irreducible representation is denoted 
by TIR) mae 


The elements (Ti ( R))u and (T3 R)) mn are related to h, li and lj 
by the orthogonality theorem as follows. 


$0», TRY =F [GF mie 


dij, Omm’, Onn’ denote the kronecker delta symbols. The kronecker 
delta symbol óij has the meaning 6; = 0 for iZj and 6; = 1 fori =j. 
This theorem covers three cases. 


1) If (Mi (R))mnand (Tj( R))mn represent two real elements 
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in the m" row and n" column of the matrix for the operation R in 
the representations i and j, then 
Unit-10 
Character table 


hl Jt, 8, = 09 (GT, (D), = h / AU j= 0 > (a) 


Th ti lied to the i ibl tation of 
e equation (a) can be applied to the irreducible representation of a NOTES 
group 
containing the operations E and A. The matrix for the operations E and 
A in the two representations are given below: 


g=] E" En a|^" An | 
A 


la al Ps » | 


Equation (a) is applied to these representations. 


Equation E11Exx + A11Axx= 0 is obtained as a result. 


2) If (Ti( R))mnis the element in the m" row and n" column of a matrix for 
operation R in the P representation and (Ij( R))mn is the element in the 
m’" row and n" column of a matrix for operation R in the same 


representation, then 


X TRY mn TRY py = (0/1), Ó,, 20 


R 


If (Ti( R))mnis the element in the m" row and n" column of a matrix for 
operation 


R in the i” irreducible representation, then 


aa», Y = All; 
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| 10.5 Direct product representation 


Suppose that R is an operation in the symmetry group of a 
molecule and X, X2, Xs ....... Xm and Y:, Y2 ....Ya .are the two 
set of functions which are bases for representations of the 


group. 


It is also true that 
RX;Y, = 25 XXX Y, 


j-ll-l 


— b» Z aX nf 
j d 


Thus the set of functions XiYx, called the direct product of 
Xi and Yk also forms a basis for a representation of the 
group. The zi, are the elements of a matrix of order (mn) x 
(mn). 


We now have a very important theorem about the 
characters of the L matrices for the various operations in the 
group. 


"The characters of the representations of a direct product 
are equal to the products of the characters of the 
representations based on the individual sets of functions". 


The above theorem can be explained by taking the 
example CV point group. The direct product 
representations of the irreducible representations of a 
group can be obtained using the character table for the 
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C3V |E 2C3 30v Unit-10 
Character table 


elii d NOTES 


The direct product representation A1A» is given by 
E 2C3 30v 
1x1 1x1  1x(-1) 
AiA2 1 1 -1 


Similarly, the direct product representations A,E,A2E and E? are obtained 
as 


E 2C3 30v 
AiE 2 -1 0 


A2E 2 -1 0 


The product representations Aı, A2, A,E and AE are irreducible 
representations. The direct product representation E? is reducible. 


The representation of a direct product, l'as, will contain the totally 
symmetric representation only if the 


irreducible Ta = irreducible l'a 
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Let la = l'ezA2 in Cav group, then we get 
Tas = A2A27 A1 (a) 


Thus in equation (a), the product representation [apis 
totally symmetric. 


10.6 Check Your Progress 


1. What is meant by representation? 

2. Define character table. 

3. What is the concern of great orthogonality theorem? 
4. How much elements represent in C2V point group? 


10.7 Answers To Check Your Progress Questions 


l. Representation is a set of matrices which represent the 
operations of a point group. It can be classified in to two types, 
which are reducible representations and irreducible 
representations. 


2. A character table is a 2 dimensional chart associated with a 
point group that contains the irreducible representations of each 
point group along with their corresponding matrix characters. It 
also contains the Mulliken symbols used to describe the 
dimensions of the irreducible representations, and the functions 
for symmetry symbols for the Cartesian coordinates as well as 
rotations about the Cartesian coordinates. 


3 This theorem is concerned with the elements of matrices 
constituting the irreducible representation of a point group. 


4 CoV point group consists of 4 elements such as E, C2, ovz and 
Ovyz, each is in a separate class. 


10.8 Summary 


* [n this unit describes the character table for a group can be 
constructed with the properties of irreducible representation. It can 
be found from the great orthogonality theorem. 

e This unit is concerned with the elements of matrices constituting 
the irreducible representation of a point group. It can expressed in 


170 


the best way for our understanding. 


10.9 Keywords 


e Types of representation - reducible representations and irreducible 
representations. 

e Character table - It is a 2 dimensional chart associated with a point 
group that contains the irreducible representations of each point group 
along with their corresponding matrix characters. 

e CV point group - It consists of 4 elements such as E, C2, ovxz and ovy, 
each is in a separate class. 


10.10 Self-Assessment Questions and Exercises 


1. Give a brief review about reducible and irreducible representations. 

2. Give a brief review about Character table of C2v and C3v point groups. 
3. Explain about the great orthogonality theorem. 

4. Explain direct product representation. 


10.11 Further Readings 


1. Chemical Application of Group Theory, F.A. Cotton, John Wiley 
and Sons Inc. New York,1971. 

2. Group theory and its applications to Chemistry, K.V. Raman, 
Tata McGraw-Hill Publishing Company,1990. 

3. Irreversible Thermodynamics, J. Rajaram and J.C. Kuriacose, Lal 
Nagin Chand, New Delhi, 1989. 
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11.0 Introduction 


Group theory is an important component for understanding the 
fundamentals of vibrational spectroscopy. The molecular or solid 
state symmetry of a material in conjunction with group theory form 
the basis of the selection rules for infrared absorption and Raman 
scattering. In this unit we learn about the Application of group theory 
to IR and Raman spectra - H20 and NH3 molecules and application of 
group theory to electronic spectra (HCHO and C2H4) 


11.1 Objectives 


After going through this unit, you will be able to; 


To know the Application of group theory to IR and Raman spectra - 
H20 and NH3 molecules 

Understand the Application of group theory to electronic spectra 
(HCHO and C2H4) 


11.2 Application of group theory to IR and Raman 
spectra - H2O and NH3 molecules 


It is not necessary for a molecule to have a permanent electric 
dipolemoment for IR absorption. But during the absorption, there 
must be a change in the dipolemoment of the molecule. This change 
in the dipole moment can be predicted in terms of symmetry using 
group theory. For this purpose direct product representation of the 
functions involved in the transition moment integral must be 
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determined. If the value of the integral is non-zero, the transition is 
allowed and if it is equal to zero, it is forbidden. 


The transition moment integral is given by 
P= [d P^ġ di 


D, and ® ‘are the wavefunctions of the ground and i" excited state of the 
vibrational level. 


P" is the electrical dipolemoment operator. 


Since dipolemoment is a vector quantity and also it is the product of charge 
and distance between the charges, it can be written as, 


P=P, + Py +P, 
P = eføixø,dt 
P =e ó*yd'dt 
" P =e ó'zó'dc 


where, ‘e’ is the charge and x,y and z are the distances in the respective 
directions. 


The value of the above integrals in general can be obtained by considering 
the symmetries, of the functions involved in the integral. Taking Pn as an 
example, the direct product representation of My’, X and V must be totally 
symmetric in order that Px £0. 


Since ®,£ is always totally symmetric, y! must have the same symmetry 
that of x. In general the symmetry corresponding to excited mode of vibration 
(®,') must match with the symmetry of any one of the Cartesian coordinates x,y 
and z for a molecule to be IR active. The above concept can be applied in the 
case of H20 and NH3 molecules. 


WATER MOLECULE 


There are three vibrational modes in the water molecule. Two of them have 
Ai symmetry and one has B» symmetry. Using group theory, it is possible to 
predict which of them are IR active. 
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A1mode: 
Px: Direct product of Tg [x Ti = AiBiAi = Bi .«.Px = 0 
Py: Direct product of Ig Ty Ti = AiB2Ai = B2 Py = 0 


Pz: Direct product of II-l; = A1 Ai A1= A1 P220 


Thus A; mode in water is IR active since PzZ 0 


B» mode: 

Px: Direct product of l's Ix [i= A1B1B? = AiA2 = A» «Px =0 
Py: Direct product of l's I yli = AiB2B2= A1A1 = Ai -Pyz0 
Pz: Direct product of Tg Iz l'i2 A1 A1 B2 = B2 P220 
Therefore B2 mode is also IR active since PyZO 

AMMONIA MOLECULE 


Ammonia molecule has four modes of vibrations, two of them 
have Ai symmetry and the other two have E symmetry. 


A1mode: 


In the case of NH3, the coordinates X and Y together 
transforms into E representations. 


PxPy: Direct product of lg Ixy Ii = A1EA1 = E ..Px Py=0 
Pz: Direct product of '5D;T;2 A1 A1 A1= A1 ..P24 0 
Since P74 0, A1 mode is IR active. 

E mode: 

P,P; : Direct product of [glx I;= AiEE = E? (41 0) 


(4 1 0) corresponds to the characters of the reducible 
representations. This can be split into various irreducible 
representations using standard reduction formula. Such an attempt 
shows that E? = A1 + A»4 E. 


Since A, forms a part of the Px, must be totally symmetric. 


<. P,P;Ż 0 
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P;:Direct product of ';D;T;i2 A1 A1E= E 
P= 0 
Since we found that Px Py Z 0, E mode is also IR active. 


It is therefore observed that all the four modes of NH3 molecule are IR 
active and all the three modes of H2O molecule are IR active. 


APPLICATIONS OF GROUP THEORY TO VIBRATIONAL RAMAN 
SPECTRA 


A charge in the dipolemoment causes the molecule to be IR-active. The 
change in the induced dipolemoment causes the molecule to be Raman active. 


The induced dipolemoment is defined as 
ui = QE 
where, & = polarisibility, E = Electric field strength 


Since a change in the induced dipolemoment will result in the change in the 
polarisibility, a molecule shows Raman spectra. If there is a change in the 
polarisibility due to the absorption of radiation. In order to apply group theory to 
Raman spectra, the value of the transition moment integral has to be calculated. 
The integral is given as, 


5,7 [a^ ejas 
where o^ = polarisibility operator 


The polarisibility operator is a measure of the quadratic functions of the 
cartesian coordinates (x^, y, z, Xy,yz,zx). For a molecule to be Raman active, 
the normal vibration should have the same irreducible representation as any one 
of the quadratic functions of the Cartesian coordinate. In other words, the direct 
product representation of any one of the following should be totally symmetric. 
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Ii Ix I2? Ii Io I? 
TiIls I1Iyl» 
Tı Fyz T2 I'1 lI»? 


[I1 is always Aj] 


This can be verified by considering water and ammonia molecule. 


WATER MOLECULE 


There are three vibrational modes for water molecule. Their 
representations are: 


Tyvib = 2A1 + B2 
A1 mode 


In order to find out the vibrational mode due to ‘Aj’ is Raman 
active or not, the direct product representations for the following 
are determined. 


DIilylI22 A1A2A1 =A2 TM Peto 
Iı Ixz T2= AiBiAi= Bi Ti Ty2T2= Ai Ai A= A1 
Ty Tyz P2= A1B2A1 = B2 Pil? 


It is noted that, among the 6 representations, three are totally 
symmetric. Therefore, Aj mode is Raman active. 
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B2 mode 


In a similar way the direct product representations involving ‘B?’ mode are Unit -11 


determined as follows: Symmetry oriented 
spectral methods 


Ii Ixy r2 = A1A?Bo- Bi [1 bo T2 

NOTES 
D; Dy P2= AiB2Bo= Ai rı Ty2P2= Ai Ai Bo= Bo 
IilDzl»-AiBiB? = A2 Dill» 


It is observed that only one among the six leads to totally symmetric 
representation. Therefore B2 mode is also Raman active. 


AMMONIA MOLECULE 


There are four vibrational modes in ammonia molecule. In terms of group 
theory, their representations have been known to be 2A; + 2E; Group theory i, 
applied to verify whether theses modes are active or not. 


A1mode 


The direct product representation corresponding to ‘Ai’ mode is determined 
in the case of the following, 


rir 
ner 


TiTyzT2 =Ai EAi=E rı r2 I2 = A Ai A12 A 


Among these, only one is totally symmetric and there is no way of finding 
out the representation for three of them (xy, x’, y?) 


E mode 


The direct product representation corresponding to ‘E’ mode is determined in 
the case of the following. 


D; rxT2} = AFE = E? ri r2r2= AiAiE=E 
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The result shows that one representation is not totally 
symmetric and the other two gives a reducible representation (E = 
4 1 0). This reducible representation can be split into number of 
representations and among them one is totally symmetric. 
Therefore ‘E’ mode is called Raman active. 


In the case of water and ammonia molecule, the corresponding 
vibrational modes are both IR and Raman active. But in the case 
of centro symmetric molecules (molecules having centre of 
symmetry) like CO», acetylene, N2F», if the vibrational modes are 
Raman active then they will be IR inactive and vice versa. This 
statement is known as mutual exclusion group. 


62 Application of group theory to 
electronic spectra (HCHO and C2H4) 


UV, Visible or electronic spectroscopy is mainly 
applicable to organic molecules. In such cases, transitions occurs 
between electronic energy levels. In the case of molecules, the 
electronic energy levels, correspond to molecular orbital energy 
levels. The various molecular orbitals are o and x (bonding 
molecular orbitals), c and x (antibonding molecular orbitals), n 
(non-bonding). The energy level diagram for these orbitals are 
shown in the figure, 


ENERGY 


Group theory can be applied to electronic spectroscopy by 
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considering HCHO (formaldehyde) molecule as an example. 


Before we apply this theory, the molecular orbital energy diagram . 
for the carbonyl group -C=O of formaldehyde should be taken in Unit -11 
t j j Symmetry oriented 
o consideration. 
spectral methods 


The diagram is shown in the figure. 


ic atom) AO (C +0) MO (O atom ) 
* 
n - 
E. d 
pA 
1 E 
E 
| 
n L as- 
Li A ea 
sp 
cies 
Too- bonds -*7—$- ^. n Pos VPE — 
ith Hydrogen atoms A. P 
[n ^AA- 


. TERI x * 
The possible transition are n-n and n-m . 


By applying group theory, it is possible to find out which type of transition is 
allowed or not. For the electronic transitions to be allowed, the following 


integrals must take non-zero value. 


[6,.xóac 
fo ygat 
føzġar 


where, b, and ®; are the wave functions for the ground state and excited states 
respectively. 
For the integral to take non-zero value, the direct product representation of 


Pand di 
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must have the same symmetry as that of any one of the Cartesian 
coordinates namely x, y, z. 


* eye 
n-z transition: 


The arrangement of electrons in the ground and excited state 
during this type of transition is shown below. 


0, 7 GP ym) 
b= (TNT) 


An electron in a molecular orbital has a particular irreducible 
representation. If the two electrons remain in the same orbital, the 
product of these representation is always totally symmetric. 
Therefore it will corresponds to ‘Ar’ representation. In such a case, 
the ground state representation is given as 


T= AIA = Ai 


In the excited state wave function, we find that there is one 
electron each in the non-bonding and m- antibonding orbitals. In 
order to find out the representations corresponding to the electrons 
in n and x orbitals, the shapes of these orbitals should be 
considered. 


The figure shows the shapes of x, x and n orbitals. 


Formaldehyde molecule belong, to C2V point group. The 
various operations present in this group are E, C2, ox; and oy;. 
These operations should be performed on m, x and n orbitals in 
order to find out the characters corresponding to each operation. 
Such an attempt leads to the following table. The last column of the 
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table gives the Mulliken symbol for each orbital. 


C2 Oxz Oyz Mulliken 
symbol 

-1 -1 1 B2 

-1 1 -1 Bı 

-1 1 -1 Bı 


Based on this table, the representation for the excited state is given as 


Ti = A1B2B = A2 


-. The direct product representation for the ground and excited state is A1.A2 
= A2. 


But this representation (A2) does not correspond to the representation of any 
one of the Cartesian coordinates. This shows that all the integrals take the value 
equal to zero, and hence n-x transition is not electronically allowed. 


But n-« has been found to be magnetically allowed. The integrals 
corresponding to the components of the magnetic dipolemoment operator are 
given below, 


fø ugar 
fouar 


[éucoóac 


The representation for the uœ), ly), Uœ are Rx,Ry and Rz which 
corresponds to B», Bı and A2 respectively. Since the direct product 
representation of d, and (4; is A» as seen earlier which 
corresponds to one of the representations here. .. one of the 
integrals take non-zero value and hence n-x transition is 
magnetically allowed. But the signal corresponding to this 

181 


Unit -11 
Symmetry oriented 
spectral methods 


NOTES 


Self-Instructional Material 


Unit 11 
Symmetry 
oriented 
spectral 
methods 


NOTE 


Self-Instructional 
Material 


transition will be very weak because the interaction of the electric 
component of electromagnetic radiation with formaldehyde will 
be very much greater than that of the magnetic component. 


n-T transition: 
The arrangement of electrons in the excited state is 
2. ut 
@j= (n) (n) (1 ) 


By a similar argument the representation corresponding to 
excited state is arrived at 


Tic Bı Aı Biz Ai 


-. The direct product representation of the ground and excited 
state is 


Ai Ai =Ai 


This direct product representation corresponds to the 
representation of the z-axis. 


-. One of the integrals take non-zero value. 


Hence a-1 transition is electronically allowed. 


11.4 Check Your Progress 


1. Define n-z transition 

2. How to prove water IR active? 

3. What are transition present in UV, Visible or electronic 
spectroscopy? 


11.5 Answer To Check Your Progress 


1. The arrangement of electrons in the ground and excited state 
during this type of transition is shown below. 


ó,2 G^) =T NAT) 


2. There are three vibrational modes in the water molecule. Two 
of them have A; symmetry and one has B2 symmetry. Using group 
theory, it is possible to predict which of them are IR active. 


3. UV, Visible or electronic spectroscopy is mainly applicable to 
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organic molecules. In such cases, transitions occurs between 
electronic energy levels. In the case of molecules, the electronic 
energy levels, correspond to molecular orbital energy levels. The 
various molecular orbitals are o and x (bonding molecular 
orbitals), o anda (antibonding molecular orbitals), n (non- 
bonding 


Unit -11 
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11.6 SUMMARY NOTES 


e Spectroscopy studies like IR and Raman the interaction of light and 
matter. The symmetry operations of a molecule form a mathematical 
group. Matrices that multiply the same way as the members of a group 
form a representation of the group. 

e  Acorrding to this reducible representation can be split into number of 
representations and among them one is totally symmetric . where the 
characters are the traces of the matrices of the representations. 


11.7 Keywords 


IR absorption - It is not necessary for a molecule to have a permanent electric 
dipolemoment for IR absorption. 

Electronic spectroscopy - UV, Visible or electronic spectroscopy is mainly 
applicable to organic molecules. In such cases, transitions occurs between 
electronic energy levels 


11.8 Self Assessment Questions And Exercises 


1. Define application of group theory to IR and Raman spectra? 
2. Explain application of group theory to electronic spectra (HCHO and C2H4 


11.9 Further Reading 


1. Chemical Application of Group Theory, F.A. Cotton, John Wiley 
and Sons Inc. New York, 1971. 

2. Group theory and its applications to Chemistry, K.V. Raman, Tata 
McGraw-Hill Publishing Company, 1990. 
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12.0 Introduction 

12.1 Objectives 

12:2 Reactions in solution 

12.3 Factors which influence the reaction rates in solution 
12.4 Application of ARRT to solution kinetics. 

12.5 Bronsted — Bjerrum equation Chelate effect 

12.6 Check your progress questions 

12.7 Answers to check your progress questions 


12.8 | Summary 
12.9 Keywords 
12.10 Self-assessment questions and exercises 


12.11 Further readings 


12.0 Introduction 

Chemistry, by its very nature, is concerned with change. 
Substances with well-defined properties are converted by 
chemical reactions into other substances with different properties. 
For any chemical reaction, chemists try to find out 


a) The feasibility of a chemical reaction which can be 
predicted by thermodynamics (as you know that a reaction with 
AG < 0, at constant temperature and pressure is feasible); 


b) (b) extent to which a reaction will proceed can be 
determined from chemical equilibrium; 
o (c) speed of a reaction i.e. time taken by a reaction to 


reach equilibrium. 

Along with feasibility and extent, it is equally important to know 
the rate and the factors controlling the rate of a chemical reaction 
for its complete understanding. For example, which parameters 
determine as to how rapidly food gets spoiled? How to design a 
rapidly setting material for dental filling? Or what controls the rate 
at which fuel burns in an auto engine? All these questions can be 
answered by the branch of chemistry, which deals with the study of 
reaction rates and their mechanisms, called chemical kinetics. The 
word kinetics is derived from the Greek word ‘kinesis’? meaning 
movement. Thermodynamics tells only about the feasibility of a 
reaction whereas chemical kinetics tells about the rate of a 
reaction. For example, thermodynamic data indicate that diamond 
shall convert to graphite but in reality the conversion rate is so 
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slow that the change is not perceptible at all Kinetic studies not 
only help us to determine the speed or rate of a chemical reaction 
but also describe the conditions by which the reaction rates can be 
altered. The factors such as concentration, temperature, pressure 
and catalyst affect the rate of a reaction. At the macroscopic level, 
we are interested in amounts reacted or formed and the rates of 
their consumption or formation. At the molecular level, the 
reaction mechanisms involving orientation and energy of 
molecules undergoing collisions, are discussed. In this Unit, we 
shall be dealing with average and instantaneous rate of reaction 
and the factors affecting these. Some elementary ideas about the 
collision theory of reaction rates are also given. However, in order 
to understand all these, let us first learn about the reaction rate. 


SOLUTION KINETICS 
When a reaction takes place in solution 


The solvent is in larger than the reactant 

Solvent concentration is constant during the course of reaction. 

Theories of solution kinetics are more complicated than gas phase. 

In solutions the collision between the reactants are termed as encounter. 
Some of the reaction has same rate in solution phase as well as gas phase Eg: 
N2Os, O3 decomposition etc. 

The rate of reaction in gas phase and solution phase are equal only in ideal 
solution 


12.1 Objectives 


After going through this unit, you will be able to: 
Understand about the kinetics and rate of the reaction 
Understand the reaction rate in solution 

Explain the concept of Bronsted and Bjerrum equation 


12.2 REACTIONS IN SOLUTION 


There are many reactions which take place more or less at the same rate i 
gaseous and solution phases (e.g. decomposition of CH2I-CH2I and N2Os i 
CCl; solvent; isomerization of d-pinene). This shows that the reactions procee 
by the same mechanism in both the phases. There are however many exceptions. 


Example: the reaction between (C2Hs)3N and CoHsBr to give a quaternar 
ammonium salt is considerably faster in polar solvents like nitrobenzene than i 
non-polar solvents or in the gaseous phase. The activated complex in th 
reaction is more polar than the reactants. 
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12.5 FACTORS WHICH INFLUENCE THE 
REACTION RATE IN SOLUTION 


(i) For ionic reaction 


A. Reaction between an ion and molecule the 
Arrhenius frequency factor has normal value. 
B. Reaction between same sign, the value of “A” is much 
lower. 


. Reaction between opposite sign, higher the value of *A". Hence 


the frequency factor is the direct measure of rate of the reaction. 
(ii) Approaching of species to each other 


Rate of approaching of reacting species to each other depend upon 
the rate of diffusion of the two species through solvent diffusion- 
controlled process that are very fast. The diffusion-controlled 
process depends upon the viscosity of the medium. In low viscous 
medium, higher the diffusion takes place. 


(iii) Cage effect 


The reactant species after the collision are held together by solvent 
species for about 10? to 10 5s. During this time lag, the reacting 
molecules collide with each other or make a suitable direction for 
the favorable reaction to take place. The reacting molecule held 
together for a period of time by solvent cage is said to be cage 
effect or frank-Rabinowitch effect. 


(iv) Energy and orientation of reacting species 


This is controlled by the nature of species. In general, the reaction 
would be fast in a solution in which the activated complex is more 
stable. Activated complex is more stable in polar solvent when it is 
polar and vice versa. 


(v) Effect of Solvation 


Solvation is based on polarity, higher the dipole will more 
solvated by polar solvent. 


Here, products are two separate ions and activated complex has 
partial charge as above. In polar solvent such as nitrobenzene there 
is more solvation of activated complex than the reactant. The 
effect of solvation lowers the activity co-efficient of activated 
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complex (y* a^) and thus rate become high. 
On the other hand, 
Y"+RX* e (v? "e Pus x? )>YR+X 


Here, there is decrease in polarity as the activated complex is formed. Activated 


complex is less solvated than reactant in polar solvent. And thus, decreases the 
rate. 


12.4 APPLICATION OF ARRT TO SOLUTION KINETICS 


For a general bimolecular reaction, 


A+B & X * > Product 
k* =Í (-.a= activity) 
a, dg 
+ 
= YO 5) 
AA Qs 


C*-k* PP. CC (2) 
y AB 


Rate of the reaction is rate of decomposition of activated complex i.e.,X* — P 


X* must have one of its vibrational degrees of freedom which would be highly 
unstable. This is responsible for the decomposition of X* in product. 


Rate = 7 [c i | 
y = frequency of vibration 


-«* H5. cc o (3) 


Any simple bimolecular reaction in solution 
Rate = K, CA Cg (4) 
Equation (3) = (4) 
KCC= yK * Fa sec 
S A B y AB 
K - y ts _, (s) 
: Y 


According to Boltzmann equation 
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in solution K,T 
hh 
NOTE Substitute in equation (5) 
K = K,T Tare K* > (6) 
EF 
For ideal solution ya, ys .y’ =1 
K,T 
K= 5 K* 
i n hence 


For unimolecular reaction, K,= K r u 
y 
have the solvation of reactant and activated complex are almost 
same both are having almost similar structure, then 
Ks = Kg 


Example: chemisorption of N2Os, Os. 


12.5 BRONSTED - BJERRUM EQUATION 


Bronsted Bjerrum explain the relationship between ionic 
strength and rate of a chemical reaction in solution. And this effect 
is involved in non catalytic reaction. 


For a general reaction, 


K n 
A? + B  ABV^*78* VE s Product 
Ky 


+ 
ao Fu — (1) 
C,C, VaVe 
c*-x*"M oco (2) 
AB * A B 
7 AB 

d kC*=kk*CC YAYB ----- (3, k3k*=ko 
dt 3 AB 3 A B z 

7 AB 
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( yAy a) Unit -12 


C Reactions in solution 
AB 


a YA yg | NOTE 
| rcs kinetic activity factor 
J AB 


For a reaction, A+ BP 


Y= K Ca Cp -------- (5), Equating (4) & (5) 
RCC =K VaYs c C , k =k Va Y a (6) 
rAB 0.2 A B r 07x 
) AB 7 AB 


Taking log on both sides, 


log k, = log ko + log Ya + log ys — log VakO A) 


According to Debye Huckel theory of strong electrolyte the activity 
coefficient 


related with ionic strength as log y= -AZ afu 


A- constant = 0.51dm 7? mol ^ for aqueous solution at 25°c 


Zi= charge of ‘i’ ion 
1 
H= bye "a Ci = concentration, Zi = valencies of species 
2 P3 
Thus log 7, = -AZ? Ju "S (7) 
log Ys =-AZp U -—---- (8) 
log ^ n- AGE Zi Aut ----- (9) 


Substitute equations (7), (8) and (9) in equation 6(A) 


k, 
log = -AZi Ju — AZ; ut A(Z, +Z)’ Ju 
ky 
--Aulz«z-G,zyl 
--Agt|Z -Z^-Z-z-ozz 
JF [z B A B A B | 
k, 
log" =+4.2Z,Z, Ju 

ko 

For an aqueous solution at 25°C this equation becomes 


189 


Self-Instructional Material 


Unit -12 Reactions 
in solution 


NOTE 


Self-Instructional Material 


k, 
log  -100Z,Z, Ju 


0 


when ionic strength is zero, kr= ko .This equation is known as 
Bronsted Bjerrum 


Ju 


equation. This equation shows that variation of ‘k’ with 


depends on ZaZp . 
The plot of log k/ko vs Ju will be linear. 
Case (i) 


When ZaZp=0 i.e. one of the reactant is zero charge (i.e. non 
electrolyte) 


Log k: = log ko, hence increase of Ju has no effect on 


Eg: [Cr(urea),|* +6H,O > Cri), +3urea, Za=3, Zs = 
0 


CHJCOOH + CNS — CH (CNSCNH+1 


Case (ii) 


k,. 
log "increases 


When ZaZp = +ve, ZaZp both has same size. 
ky 


with increasing 
TE SO% «21 —2I*280^,Z,Z-22 
H 2 8 2 4 A B 
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2[Co(NH ) Br |+ Hg ** + 2H O 5 2[Co(NH ) H O]+ HeBr 
35 2 35 2 2 

Z,,2,-4 

Case (iii) 


k, eon ; 
When ZaZp=-ve l0 decreases with increasing Ju 
ko 


(i) HO +2Br7+2H*—>2H0+Br Z,Z =(-1)(1)=-1 
2; 2 2 2 A B 


[Co(NH,)Br}* +OH~ —[Co(NH,)OH |+ Br- Z,,2,--2 


12.6 Check Your Progress 


1. What is mean by rate of the reaction? 
2. How does the rate of reaction depend on volume? 
3.Explain Enzyme catalysts? 


12.7 Answers to Check Your Progress Questions 


1. The rate of a reaction is the speed at which a chemical reaction happens. 
If a reaction has a low rate, that means the molecules combine at a slower 
speed than a reaction with a high rate. Some reactions take hundreds, maybe 
even thousands, of years while others can happen in less than one second. 


2.The rate of reaction r for the reaction AB is given by : r=k[A],where k is 
the rate constant and [A] is the concentration of the reactant.(I've assumed 
the reaction to follow first order kinetics.)So r is directly proportional to the 
concentration of the reactant. rate increases on increasing reactant 
concentration. However, the concentration of A is actually in moles/Liter, 
which is moles per unit volume. Hence, r is inversely proportional to the 
volume and the rate decreases on increasing volume of the reaction vessel. 


3. Enzyme catalysis is the increase in the rate of a process by a biological 
molecule, an "enzyme". Most enzymes are proteins, and most such 
processes are chemical reactions. Withinthe enzyme, 
generally catalysis occurs at a localized site, called the active site. 


12.8 Summary 


Y Chemical kinetics is the study of chemical reactions with respect to 
reaction rates, effect of various variables, rearrangement of atoms 
and formation of intermediates. 

v The rate of a reaction is concerned with decrease in concentration of 
reactants or increase in the concentration of products per unit time. It 
can be expressed Application of ARRT to solution kinetics, 

Y Mathematical representation of Bronsted Bjerrum equation.It has to 
be determined experimentally and cannot be predicted. 


12.9 Keywords 
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Cage Effect: The cage effect in chemistry describes how the 
properties of a molecule are affected by its surroundings. First 
introduced by Franck and Rabinowitch in 1934, the cage 
effect suggests that instead of acting as an individual particle, 

molecules in solvent are more accurately described as an 
NOTE encapsulated particle. 
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Solvation effect: If the polarity of the product is different from 
that of the starting material, solvation changes the thermodynamic 
properties of the reaction. If the transition state experiences the 
change in polarity (usually charge buildup), solvation changes 
the kinetics properties of the reaction 


12.10 Self-assessment questions and exercises 


1. Explain in detail about the rate constants. 
2. Discuss in detail about reaction rates in solution 
3. Distinguish the Bronsted and Bjerrum equation 


12.11 Further readings 


L. C.M. Guldberg and P. Waage,"Studies Concerning 
Affinity" Forhandlinger i Videnskabs-Selskabet i 
Christiania (1864), 35 


2, P. Waage, "Experiments for Determining the Affinity 
Law" ,Forhandlinger i Videnskabs-Selskabet i Christiania, (1864) 
92. 

3. C.M. Guldberg, "Concerning the Laws of Chemical 


Affinity", Forhandlinger i Videnskabs-Selskabet i 
Christiania (1864) 111 
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13.8 Answers to check your progress questions 
139 Summary 
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13.11  Self-assessment questions and exercises 
13.12 Further readings 


13.0 Introduction 


In our last two modules, we took up solution kinetics in detail. You 
now know that solution kinetics is the kinetic study of chemical 
reactions taking place in solution phase and the nature of solvent has a 
predominant influence on the kinetics of such reactions. 

In general, the solution phase reactions can be 
classified as Ionic reactions and Non Ionic reactions. And here, we are 
concerned with the kinetic study of ionic reactions. 

A-B-—P ...(1) 
Rate = k[A][B] ...(2) 
where k is the rate constant for the reaction. 


We found that only Classical Thermodynamic 
approach of transition state theory can be utilized to determine the rate 
constant of ionic reactions or solution phase reactions and studied the 
effect of solvent on rate constant. In this module, focus will be laid on 
primary salt effect and secondary salt effect. 


13.1 Objectives 
After going through this unit, you will be able to: 
e Understand the concept behind the salt effect in the solution. 
e Understand the mechanisms behind effect of pressure and 
volume of activation in solution. 
e [earn about the Effect of substituents on reaction rates 
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13.2 PRIMARY SALT EFFECT 


By doing a neutral salt to the reaction mixture having 
either type of ions, the rate of the reaction is considerably 
affected according to the Bronsted Bjerrum equation is said to 
be primary salt effect or neutral salt effect. 


It is experimentally found that the rates of second-order 
reactions between charged species are affected strongly by the 
ionic strength of the solution, this quantity being defined as 


]ed2x6zZ. 


where Cj is the concentration of each ion present in the solution 
and Zj is the charge number of the ion. When the reacting ions 
are of same sign, increase in ionic strength increases the rate 
whereas when the reacting ions are of opposite signs, increase in 
ionic strength decreases the rate. 


13.3 SECONDARY SALT EFFECT 


This arises if one of the reactants is a weak electrolyte that can 
be stimulated by catalytic reaction. The rate of the reaction 
depends upon the amount of salt added in a catalytic reaction. 
Actual change in concentration of reacting ion by adding the 
salt externally. 


This effect is not having a direct influence. But it will influence 
the concentration of reactant. The equilibrium constant for the 
dissociation of weak acid HA can be written as 


leder 
[HA] Yn 


Addition of salt will influence the activity co-efficient they also 
influence the concentrations of H+, A- and HA and this way the 
catalytic activity has been lowered/increased, will influence rate 
is said to be secondary salt effect. 


13.4 INFLUENCE OF INTERNAL PRESSURE 


We know the rate constant in solution has the relationship with 
gas phase rate equation. 


K DET. AVE: > (1) 
Uh y 
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K,T 


K= * K*—(2) 
E h Unit — 13 
Salt effect 
Ko = rate constant under ideal conditions. 
For an ideal solution yy yu yt #1 NOTE 


activity co-efficient according to langmuir is given by 


RTIn y -V(P'^ - p? 
1 1 1 2 


y=ViA 
Vi— Molar volume 
Pı — Pressure of solvent 
P2— Pressure of solute 
RTInyi=ViA > (3) 


RT In ya= Va Aa; RT In ys — Vs Ap; RT Iny = V^ A* 


RT In ya + RT An y- RT In y” = Va Aa+ Vn AB- V^ A 
rTin' ^4£2! =v, Ant VgAn- V* A* 
7) 


K, 
RTIn — — v, Aa + Va An- V^ A^ 


The molar volume of reactant and activated complex are equal. 
Case (i) 
If Aa = Ag = A? 

Then the solvent will not affect the rate of reaction in solution. 
Case (ii) 

P:olvent = Preactan and A* is large 
Then the reaction rate in solution would be low than ideal solution rate. 
Case (iii) 
When Poolvent & Pactivated complex , has much different from that of reactant 


than the RT Ln K;/Ko will be positive value and make the reaction rate in 
solution high, that is RT Ln Ksis large. 


Thus the kinetic study of solution is very complicated but in general, 
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(i) If the product of the reaction is polar, then the reaction is favoured in 
polar solvent. 


Unit — 13 (ii) If the product of the reaction is nonpolar, then the reaction 
Salt effect is favoured in non-polar solvent. 


(iii) If the product of the reaction having high internal pressure 
NOTE then the reaction is favoured in solvent having high internal 
pressure. 


13.5 EFFECT OF SOLVENT 


It can be classified in the following topics 


13.5.1 REACTION BETWEEN IONS 


The electrostatic forces between ions are much stronger than 
non electrostatic, forces. The pre exponential factors of ionic 
reactions depend in simple electrostatic principle. If ions are 
oppositely charged the pre exponential factors are abnormally 
high, where as if the charges are the same they are abnormally 
low. 


The reaction between ions in solution has markedly affected 
by solvent and its dielectric constant. In order to explain this 
effect there are two models into account. 


(i) Single sphere model 
(ii) Double sphere model 
Single sphere model 


Here slightly different relationship has been obtained from that 
of double sphere model. 


a) The reacting ions are regarded as becoming merged into one 
single sphere which has charge equal to the algebric sum of the 
both the ions. 


b) The rate equation for this model is derived by Born. 


c) Consider the process of charging a conducting sphere of 


6699 


radius “r” from an initial charge if zero to a final charge Ze. 


(d) This process is carried out by transporting from infinite 
distance to small increment of charge equal to edA. 
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À = Parameter which varies from zero to ‘z’. 


At any time the charge of the sphere has written as Ae, and if at a given 
instant the increment of charge is at a distance x from the ion the force acting 
on it is 
eA 
je a 
4I, £x 
AeedA | AedA 
df = 2 = ; — (1) 
AMEX — ATE, EX 


The work of moving the increment from x to dx is 


Ae dAdx 
dw = 
ATE, EX 


The total work of charging is obtained by carried out by double integration 


A 


0>Z? Xar 

=ë? oz d 1 
jus d [aa 2t a k 
o a x , "n pem E xl 

2 

n 1YA ee ez 3) 

E — , We g ree 

: FAT SREE 
Ga m "e t = (4d 4 ye 
87E Er | AEE 
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The | charge in electrostatic AG, is 
2 2 xe Z>? 
wes | amer e T sa) T Ia 
A 


8 EE 5 Y 
Salt effect \ Y í s) 


Free energy per molecule is 


NOTE AG” „= Ne (Z RE 7 Z E Z^ 


8zEE( Y ye Fal 
AG* es — AG* es + AG nes . 


G°”., = Electrostatic contribution to the Gibbs energy of 
activation 


G” „s = Non-Electrostatic contribution to the Gibbs energy of 
activation 


Substitute equation (5) in ARRT equation 


K KT ag" irrr 
Ne) (Za Za) Z1, Zig) 


ORT sman Gnon]3 9 (6) 
K, m Nh* .€ 


Taking natural logarithms we obtain 


kom RD, H e (Zit Ze) zu zu) 
2h SZEERT( Y Ya Ya 


This may be written as 
2 (72 2 
e Z°+Z°+2Z Z 2 2 
In k,= 1n kọ- |a na a 2-2/4 (8) 


SMEERT V y Ya Ye) 


, 


when Y" — y, = Ys 
The equation becomes double sphere model 


Double sphere model 


According to ARRT, general reaction is as follows 


Self-Instructional Material Ions + A.C product 
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RT | AG* 
k -k"* ^ , aG*t=-RTInk*,  k*=e™ 


It is a thermodynamic statement of ARRT 


The charged ions in solution are considered to be 
conducting sphere and solvent regarded as fixed dielectric 
constant and hence the spheres are hard and rigid 


Y 


} 


ya and yg = Radius of two ions, Zae and Zge = the charges on the 
ions 


e = 1.602x10'? C. 


Initially the ions are infinite distance apart. In this particular model 
known as double sphere model (hence no force between them) 


When the ions are separated by a distance ‘x’ the force acting 
between them is 


according to coulombs law 


Work done in moving two ions together to a distance dx is given 
by 
Zhe, 
== dx , ‘-ve’ sign implies the decrease of ‘x’. 
Ame, Ec 


dw 


The work done on the system in moving the ions from x = @ to x = das 


2 2 dAB 
B e ZZe 1 
Jaw =-J zi 7 , w=— AB _,dx 
a X 42&,€ 4 NEE” yX 
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Ze 2n zs 
w= ee | 1j w=4 ^d I —— 
4 NEE p xla 4€, €d sz 


(3) 


This work is +ve when the ions are same sign. This is — 
ve if opposite charge. Hence the work is electrostatic 
contribution to the Gibbs free energy of activation. 


AG ZZ AG. 
= £ g (electrostatic) + Gis (nonelectrostatic) 
N 4 TNE Edag N 
ZZe " 
Ag Ed aA LLL (4) 


o AB 


Substitute in equation (1) 


_AG* | ZaZge^N NY 
k,2——e m, k- RT e Tdi 4 RT 
Nh : Nh 
RT Ge o ZAZge 2N \ 
Ink, =ln|—e rre imtme 
T -59* ZZeN 
In k,=In—e R^ ——4—B8.—— -------- (5) 
: 4 é,€dy,RT 
k ZZeN 
ln. = —4 2 


ky  4€,€d,,RT 


Ko = rate constant of reaction in solvent of infinite dielectric 
constant. Plot of kr vs l/e will give straight line. It is 
experimentally verified but there is a deviation in low dielectric 
constant. From the slope the dag can be calculated, it has the 
the value of few hundred picometer (10 ^). 
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Ion-dipole 


In reality, many ions have a distribution of charges, and this makes a 
difference to the electrostatic interactions. Much interest is between ion and 
dipole molecules or between two dipole molecules. 


The electrostatic free energy contribution is from single sphere model 
Sphere of radius — r 
Net charge - Ze 


Dipole moment —u 


A's pp — 36 () 
G,, ——F LÁ 
8 7E,Er loneer 


For a bimolecular reaction 


A +B > X* > product 


3 
9 E j M 3 5f ) 
AG,,(A) = ——— 
877€ &r^ o A 
ze " 34^ > Q) 
o B 
AG,, (B) = TOMET” 
i 8 meer? gu 


A s E y e 34 9 Q) 
G..() —8z2EEF— | TOME £r^ 


So the electrostatic contribution for activation 


" poss r PAULI aga Por 
o # A B) o N * A T, 
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The non electrostatic contribution say AG nes 
AG* = AG" es + AG" nes 
From ARRT 


Ink = Ink,— 


Assume radii do not vary as much as the charge and dipole 
moments, it is a useful approximation to treat them all the 
same and equal to dap. 


Ya= Yr = Y«7 da 
Then equation (6) becomes 
z 
ge www +() 
A B 7 


Ink nK, "NUM 
4nEEd KT  loZE&l KT 
o AB o AB 


If both the species are charged the second term in equation (6) 
or equation (7) is usually much greater than the third. 


If either reactant has no charge, however, the final term 
becomes predominant. 


If the reacting species are uncharged, converted into activated 
complex has more polar than reactant (yzis large compared to 
YA and Yg, the rate constant increases with increasing dielectric 
constant. 
A P 
Au^ 2 w-uw-m () 
(1.602x10? C) 


1.602x10 ? Cm is the dipole moment of two elementary 
charges (1.602x10 ^C) separated by a distance 10 7m. 


13.6 EFFECT OF SUBSTITUENTS ON 
REACTION RATES 


How the rate of the reaction or equilibrium constant will 
affect by substitution in either m or p was explained by 
Hammet. It is a linear relationship between Gibbs free energy 
of any two set of reactions. 
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13.6.1 HAMMETT'S EQUATION 


According to Hammetts relationship the rate constant for the reaction 
of one compound is related to that for the unsubstituted parent compound in 
terms of two parameters p and o 


For rate constants the relationship is 
k 
log— = po 
ko 


© - Substituent constant is a number (ve of —ve) 
+ve = for electron withdrawing substituent. 
-ve = for electron releasing substituent. 
p - reaction constant is n number due to external condition, solvent etc. 


p is +ve if solvent will assist the electron withdrawing nature or releasing 
nature. 


p is —ve if solvent will against the electron withdrawing or releasing nature. 
Co6HsCOOH +> C;HsCOO' +H* p-l (constant) 
NO» -C6H5 COOH + NO2 -C6H5 COO’ + H* om=0.710 
log = po,- 0.710 xl 


0 


log = 0.710 


0 
log k= log kọ + 0.710 
uu log ™ 0.710 
ko 
un 5.31 
ko 
k=k,5.31 


ie. Rate of ionization of substituted benzoic acid is 5.31 times greater than 
unsubstituted benzoic acid. 


Using the same value of p and o for hydrolysis of benzamide in 60% 
ethanol at 80°. We can calculate the m-nitrobenzamide which can be 
hydrolysed 0.615 times as fast as benzamide. 
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| log k = log k,+ po, 


When we plot log k vs osis a linear relation, from that we can 


Unit — 13 calculate the value of p. 


Salt effect 
NOTE 
slope= P 


log K 


log Ko 


E 


Linear free energy relation of Hammet equation: 


The rate constant of a reaction related with free energy of 
activation by. 


K,T 5 
k= È e RT ec (1) 
h 
AG; 
wc p (2) 


Taking log for equation (1) 
k T, D 
log k = log ^. ae AP 
h | 2.303RT 
Hammett equation 


log k =logk)+ po (4) 


Substitute (3) in (4) 
T AG T AG" 
__—-———— =logk,— -. *po 
log ke- 7 303RT © ^ 2303RT 
AG* AG 7 


———— = = pe x 2.303RT 
2.308RT  Z.303R 


AG* = AG, —2.303RTpo (5) 


"" —-— 
- un 
Soininen Malesia AG* = free energy of activation for substituted compound 
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AG?? = free energy of activation for parent compound 


Equation (5) with a particular value of p, applies to any reaction involving a 
reactant having series of substituents. 


For another series has different p value and by assuming © is constant. 


AG* =AG,* —2.303RTp'o 


+ 
AG -AG I3O3RTO uror (6) 
p 
Equation (5) may be written as 
AG AG” 
+  -2——-230NTO ~~~ (7) 
p 
Subracting (7) — (6) p 
AG* * 
o 0 
P P 
AG" — p AG* =constant. 
' 0 
"P 


Then there is a linear relationship between Gibbs energy of activation with in 
any two set of homogeneous reaction. 


13.6.2 TAFT EQUATION 


When Hammet plot of log K/Ko against o are extended to the reactions of 
aliphatic compounds and to those of o-substituted benzene derivatives the 
straight line no longer result. Thus Taft derived another equation, including 
polar factor, conjugation factor and steric effect hence used for aliphatic 
compound also. 


log K/Ko ZP-C«S 
P — polar factor 
C - Conjugation factor 
S — Steric effect. 


He studied acid hydrolysis of ester as well as base hydrolysis of ester, single 
equation is not hold good for both the case. The polar effect has neglected in 
acid hydrolysis, 
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log (K/Ko)acid = C + S 
log (K/Ko)base = P + log (K/Ko)acia 


log (K/Ko)base - log (K/Ko)acia= P = P* 6^ 


This equation is known as Taft Equation. 


13.7 Check Your Progress 


1. What is the difference between the primary and secondary salt ? 
2. What is the effect of ionic strength on the rate of reaction? 


13.8 Answers to check your progress questions 


1. Primary salt effect. It has been observed that the rate of a reaction 
can be altered by the presence of non-reacting or inert ionic species 
in the solution. This effect is profound when the reaction takes 
place between ions, even at low concentrations. 

The term “secondary kinetic salt effect" is used to indicate a 
kinetic. salt effect due to a change in concentration of the reacting 
molecules on. account of a change in the inter-ionic forces.1,2. 
Such effects are very commonly met in the study of kinetic 
reactions. 


2. This question might be simple for those who might have 
understood the concept of CHEMICAL KINETICS in which there's 
a beautiful equation given to us with which we can find out rate of 
reaction in addition to this entropy also changesi.e with increase in 
IONIC strength of a compound the chemical reaction 

rate decreases.some factors affecting the reaction rates might also 
be included such a s temperature,pressure and concentration mainly 


13.9 Summary 


v In this module, focus was laid on primary salt effect and 
secondary salt effect. 
v It has been observed that the rate of a reaction can be 


altered by the presence of non-reacting or inert ionic species in the 
solution. This effect is profound when the reaction takes place 
between ions, even at low concentrations. This influence of charged 
species on the rate of the reaction is referred to as salt effect. 


v The salt effect is classified as Primary salt effect and 
Secondary salt effect. 
Y The primary salt effect takes into account the influence of 


electrolyte concentration on the activity coefficient and hence the 
rate of the reaction. 
206 


d The secondary salt effect is the actual change in th 
concentration of the reacting ions resulting from the addition of Unit — 13 


electrolytes Salt effect 
13.10 Keywords 
Primary salt effect- It has been observed that the rate of a reaction can be NOTE 


altered by the presence of non-reacting or inert ionic species in the 
solution. This effect is profound when the reaction takes place between 
ions, even at low concentrations. 

The secondary salt effect- It is the actual change in the concentration of 
the reacting ions. resulting from the addition of electrolytes. 


13.11 Self-assessment questions and exercises 
1. Explain in detail about effect of solvent with example. 


2. Explain the dielectric constant. 
3. What is meant by Hammett and Taft equations? 


13.12 Further readings 


1. Mendham, J.; Denney, R. C.; Barnes, J. D.; Thomas, M. J. K. 
(2000), Vogel's Quantitative Chemical Analysis (6th ed.), New York: 
Prentice Hall, p. 28, ISBN 0-582-22628-7 

2. Physical chemistry, Peter Atkins, Julio De Paula, gh edition, 
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Unit-14: ACID BASE CATALYSIS 
Structure 


14.0 Introduction 

14.1 Objectives 

14.2 Acid base catalysis- 

14.3 Acidity functions 

14.4 Bronsted relations 

1.4.5 | Zucker Hammett hypothesis 

1.4.0 Enzyme catalysis 

1.4.7 Effect of pH 

14.8 Influence the Temperature on enzyme catalyzed reactions 
14.9 Check your progress questions 

14.10 Answers to check your progress questions 
14.11 Summary 

14.12 Keywords 

14.13  Self-assessment questions and exercises 
14.14 Further readings 


14.0 Introduction 


The acid catalysis and base catalysis, a chemical 
reaction is catalyzed by an acid or a base. By Brgnsted—Lowry 
acid-base theory, the acid is the proton (hydrogen ion, H^) donor 
and the base is the proton acceptor. Typical reactions catalyzed by 
proton transfer are esterification’s and aldol reactions. In these 
reactions, the conjugate acid of the carbonyl group is a 
better electrophile than the neutral carbonyl group itself. Depending 
on the chemical species that act as the acid or base, catalytic 
mechanisms can be classified as either specific catalysis and general 
catalysis. Many enzymes operate by specific catalysis. 


14.1 Objectives 


After going through this unit, you will be able to: 


e Understand about the Acid and base catalysts. 

e Understand the methods of determining the Enzyme 
catalysis 

° Explain the concept of effect of pH and temperature on 


enzyme catalyzed reactions 


14.2 ACID BASE 


A reaction is catalysed by H* or H3O' ion is said to be 
specific acid catalysis reaction. Some reaction is catalysed by 
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any Bronsted acid is said to be general acid catalysis reaction. 
If a reaction is catalysed by OH- ion is specific, base catalysis 
or by any base are said to be general base catalysis reaction 
Some reaction catalysed by both acid and base are said to be 
acid-base catalysis. 


If a reaction is carried in Ens acid the r 


r=K,[S ] - K, is lo) 


If OH also included 
r-K [s]+K bd ]e x son ]2 (2) 
Divide through out by [S] 
K-K +K |H* |+ Kalon |> (8) 

Note: 

r- K [s ] 

r 

[3 

Skrabal diagram: 
The plot of log K vs pH of the solution is known as Skrabal diagram. 


For mutarotation of Glucose 


CASEI 


For acid catalyzed reaction: 


K-k,|m] 


LogK = LogK ,.+ Log|H" | 


LogK = LogK „+ — pH 
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CASE II K=K_ [on | 


LogK = LogK ,,,.  Log| OH | 


K 
LogK = LogK y + Log." 


LogK = LogK ,,,-K,,+ pH 


OH 
CASE III 


The horizontal line with zero slope corresponds to uncatalysed 
reaction. Zero order with respect to acid/base 


14.3 ACIDITY FUNCTIONS 


The first acidity function proposed, other than the pH, is due to 
Hammett and Doyrup, is based on equilibrium of the type 


B+ H* <> BH* 
Eg: C,H,NH, + H ' 6 CH ,NH, 
Equilibrium constant for such reactions is 
[sa] Y. 
x 
(fae 


Where the Y’s are the activity co- efficient. Where B, [BH*] can 
be distinguished 
H *] 


Spectrophotometrically, it is possible to measure , and K 


also can be measured in dilute solutions. 
Taking log in Equation (1) 


[BH *] a 
lo BH 


log k = log ———— + log____ 
[BIH] YY, 
[B + ] + P 
log k - log —-log[H ]-logY . + log 
[B] 7 Y, 


Y, 
log k — log [BH "] -= -logga  -log * 
[B] e 


BH* 
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log k — log [BH *] =-loga , ES 
[5] = m 


Fi 
Y dH Sas, 


Left hand side of equation can be measured by experimental way. The 
function on the right hand side can be measured in acidic solution by 
introducing a suitable indicator and measuring the concentrations of the two 
species. The quantity on the right-hand side. 


H =- log | Gy Vs | Hammett acidity function 


BH* 
S-H* SH+ (rapid) >[SH* 7 
SH * — product(slow) 


If the second step is slow and rate controlling step, the overall rate is 
proportional to the ion concentration of activated complex [(SH*)]* 


y-k"[SH^Y 
. [SH 'T Y* 
Rss cete 
[SH+] Y, 
K/[SH "ly. 
K*- A en. [SH ^l = S SET Mtt. (2) 
` Coy Y 


a’ and Y, are activity coefficient. 


T a. [SH']Y. 
For the pre equilibrium, j su = s /—  —  -—— (3) 
EP 7 [S] Y;a 


From equation (2) y-k^[SH T 


k*k*[SH* ]Y 
rate - * ” 
y* 
| [SH * ] Y. 
From equation (3) k =s 
° [S ] Ysa,. 


14.4 BRONSTED RELATIONS 


Since catalysis by acids and bases usually involves the transfer 
of proton from or to the catalyst, it is natural to seek a 
correlation between the effectiveness of catalyst and its 
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strength as an acid or base. The most satisfactory relationship 
was given by Bronsted termed as Bronsted relations, 


(i) The catalytic constant is related with dissociation constant Ka 
by, 


a a a 


Where, Ga and a are constan,  kacatalytic constant, Ka 
Dissociation constant, a always less than unity. 


(ii) Similarly, for a base catalysis 
k =GK* 
b b b 
| (1Y E 
Gell : E 
Da 


where k» is the dissociation constant for the base; ka that for its 
conjugate acid; p is a constant that is again less than unity. 


Modification of above equation: 


If we applied this relation to an acid catalysis reaction of 
dibasic nature i.e. acid has more than one ionizable proton or base 
more than one centre to accept proton, the modification should be 
made. 


CH3(CH2)COOH - Long chain fatty acids 
HOOC-CH2(CH2)} COOH - Dibasic acid 


(i) In the dibasic acid there is negligible interaction between the 
two carboxyl groups. The dissociation constant of the dicarboxylic 
acid is twice that of the monocarboxylic acid. Since the ion can be 
formed by loss of either of two proton. For the same reason the 
catalytic activity or catalytic constant for dibasic acid is twice that 
of monobasic acid. 


Since the ratio of 2 in the acid strengths leads to the ratio of 2" in 
the catalytic constant and this is generally less than 2. This anomaly 
is avoided if both the acid strengths and catalytic constant are 
divided by the number of proton involved in the dissociation. 
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(i) In the case of the two acids HOOC(CH2)nCOOH and 
HOOC(CH2)COO , both of which have one dissociable proton, the 
catalytic strength are the same. On the other hand, the acid strength of 
second is one half of first since the ion OOC(CH2).COO in which the 
second dissociates has two points at which a proton may be added, where 
as the ion of the first acid has only one site. To remove the inconvenience it 
is necessary to multiply the dissociation constant of HOOC(CH2)COO by2 
before inserting it into the equation 


k (qK Y 
pop, 
p 
Ed 
q NK j 


p — number of dissociable proton bound equally strongly in the acid, while 
q is the number of equivalent positions in the conjugate base to which a 
proton may be attached. 


The Bronsted relationships are special cases of the linear Gibbs energy 
relationship. 


log k = log K,+ óp 
log k = log K, -óp 
5 is same in both the cases but the reaction constant p and p` are different 
in equilibrium and rate constant equations. 
1 1 
_logk=_logk+ 6 
p p e 
1 1 
and —logk- _logk+ ô 
p p 
1 1 
Subtraction leads to —108k - —log k = constant 
4 
k? 
And therefore log —r = cons tan t 
K^ 
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This may written as log E — constant 
Unit — 14 K^ 
Acid Base catalysis 
or k=GKo 


NOTE 14.5 ZUCKER- HAMMETT HYPOTHESIS 


The acidity function is used to check the existence of 
correlation between Ho and rate constant, a pre-equilibrium of the 
type as give below can be established. 


K 
X+H* <> XH* (fast) XH‘ — P (slow) 


The second step is slow and rate determining step. 
The activated complex has equilibrium with further activated 
complex. 


K! 


XH* &y[XH'Y' 5p 


K'XH*]y . 
From (2) [XH Ý »—^7 EE dits (3) 
y* 
kk XH* Vy, 
aa rcov (4) 


[XH ^]y. 
[XH ly, Yx 


From the pre equilibrium k,, = 


[xe = T Pan 


Meran (5) 
kk k y[X]a y 
Substitute (5) in (4) y= 2 eq s H*XH* 


XH* 


kk k ya 
y= e xati X] 


y* 
Self-Instructional Materiųl 
For a first order rate, k= r /[s] 
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kk kya 


k = ened XH 


Y Unit - 14 
[S] 2 IX] Acid Base catalysis 
Taking natural logarithm for the above 
a. 

logk=logkk k +log VG sa (6) NOTE 

2 eq eq y 
logk=logkk k —H--------- (7) 

2 eq eq 0 


Ho is acidity function. The equation shows that a correlation will exist 
between log k and Ho, reaction mechanism involving pre-equilibrium.This 
same mechanism tentatively modified by ‘Zucker’ by including water 
(solvent) in the slow process as follows 


K K 


X+H* < XH+ [XH+] > product Slow/H20 


follow the derivation equation (6) becomes, 


' da a 
logk-logkk k pig” Rem 
2 eq eq E NE 


The suggestion was that it might be possible to decide between the two 
mechanism by seeing whether the rate constant showed better correlation 
with p" or with Ho, Unfortunately this Hammett-Zucker hypothesis did not 
prove reliable. 


14.6 ENZYME CATALYSIS 


Enzymes are proteins with high molecular mass of order 10000 or even 
more they are derived from living organism. 


Enzyme catalysed reactions are specific in nature. 
CH O V4"352E(0H +2C0 
6 12 6 2 
urea LAE" 2NH, + CO, 


Mechanism and kinetics of enzyme catalysed reaction 
Ky 
E+S SES 
K-1 
[ES]|LY—P+s 


E= enzyme, s =substrate 
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Apply steady state approximation for [ES ]. 
Unit- 14 K [E][S]- K-,LES] - K,[ES]=0---------- (1) 
Acid Base catalysis 
[ES]K + K,= K,[E][S] 
NOTE [ES] = K,[E]LS ] > (2) 
rd 4 
[E] = [Eo]-[ES] 
[E]=Equilibrium constant 


Substitute in equation, Ki[S][(Eo)-(ES)]=K-.i[ES]-K2[ES] 


[ES] = K[S][Eo] Rate of the reaction r = k [ES] 
K_,+K,+K,[S] ? 


- KKISIE] | 
K,[S] + K , - K, 


Divide the equation (3) by ki 


r-2  Kj[E,IIS ] 


K, 
[S]+ K- + K, 


r= KEJIS] -------- (5) 
[S] +K, 


where km = k-1+k2/kı --- Michaelis constant. This rate equation is 
known as Michaelis -Menton equation. 


At very high concentration of substrate the rate will be 
maximum. 


From equation (5) km««[S] 


„kE ILS] „ ELEME] 
[S]+K,, db [S] 


Substitute (6) in (5) 
eS -Ja e 
[S]+K,, [S]+K,, 


Limiting cases: 


Case (i) 


Self-Instructional Materiql 
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When km >> [s] 


r= Y maz [S] — k [5] 
K, 


m 


Reaction is first order with respect to [s]. 


Case (ii) 


[s]>>km, p = y, = constan t 
s 
Case (iii) 
If s =km, Vmax [5] _ 1 


r= == 


2i] .2** 


When we plot a graph Kmis defined as concentration of substrate at which 
the rate of formation of product is half the maximum rate at high concentration 


of substrate. 


Why the rate of enzyme catalyst reaction changes from first order to 
zero order as the concentration of [s] increases? 


Each enzyme molecule has one or more active sites at which the substrate 
must be adsorbed in order that catalytic action must occur. 


At low [s] 


As the [s] increases the number of sites which are occupied increases and the 
rate also increases. 


At high [s] 
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All the active sites are already occupied and the further increase in 
[s] can not further increase in rate and become zero order. It is 
very difficult to determine the Ymax and Km from the plot y vs [S]. 


(i) Line Weaver-Burk Method: 
NOTE 5 Vmax lS] 


Unit — 14 
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In F Ps Int 
MM ———À, dBIK, 
r Ymax [5] 
1 1 +_k 
~ TA Jas. Iis] Plot between 1/r vs 1/[s] 


14.7 EFFECT OF pH AND TEMPERATURE ON ENZYME 
CATALYZED REACTIONS 


Influence of pH 


The pH of the solution usually has a very marked effect on the 
rate of an enzyme reaction. In most cases the rates of enzyme 
reaction passes through a maximum as the pH is varied. The pH 
corresponding to the maximum rate is known as optimum pH. Its 
value varies with the nature of substrate and with the substrate 
concentration. 


Effects of pH are irreversible if the acidity or basicity become 
too high, since the ternary substance of the protein is destroyed. 
Reversible pH changes occur when the pH is not taken too far from 
the pH optimum. 


With in a certain pH range the pH can be changed back and 
forth without any permanent affects ensuring. 


This behaviour was first explained by Michaelis to postulate at 
least two ionizing groups as playing an important role at the active 


Self-Instructional Materitil E . 
j centre, these groups are -NH3* and -COOH, the ionization at 
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active centre may be represented as 


COOH NH; COO NH; COO^ 
\ / X / \ 
—— L a tos af | 
| Enzyme | —= | Enzyme | ———— O Enzyme 
I —— | L-.-- Li 
(EH) (EH) (E) 
— Low pH High pH — 


The pH behaviour can be explained by postulating that the intermediate, 
Zwitter ion, form is enzymatically active, but that the species to left and 
right are inactive. The concentration of intermediate goes through a 
maximum as the pH is varied so that the rate passes through a maximum. 


Mechanism: 
(a) 
S 
K, * Eo 
So BH gi 
k IE 
K RI ..K 
EH;S = -> EHS ES 
kzl 
"+ 
EH +Z 
(b) 


Applying steady state approximation for the above mechanism . 


At low substrate concentration 
[Eo] KS[S] | 
| 147 + H n 
MS 
N H K, ) | 


"E 
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[H *] 


If we study the reaction at low pH is very high. 


Unit - 14 K, 
Acid Base catalysis 
Equation (1) becomes, r= kalbo MS) c (2) 
NOTE K,| | 
© ul 
kk, [Eo ILS] 
ee 3 
r K [H t ] =} ( ) 


Taking log for equation (3) 
log r = constant — log[H*] 
log r = Const +pH is at low pH . 


(ii) At intermediate PP the km value cannot be neglected from 
equation (1) 


k,[E, ILS 
r= AEST — (4) pH has no effect on rate 


m 


(iii) At high pH value the term k,/H" is large in equation (1) 


_ K;[E,][S] 
K, 11-5 QUE] 
(8E KJ 
r= KEJIS] ------ (4a) 
mae 


log r = Constant + logH* 


log r = Constant — pH — (5) 


Self-Instructional Materiql - e mE ur 
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In the graph AB has +ve slope satisfy the eqn (3a) 


Unit - 14 
BC is by equation (4) Acid Base catalysis 
CD is by equation (5) 
At point B the eqn (3a) and (4) are equal NOTE 
kK [Eo JIS ] _ LES IS ] 
KH K 
[H*] = Ko 
pH = PK» 


At point ‘C’ the equation (4) and (4a) are equal 
k [E ][S] k[E ]ESTIH *] 
2 0 — 2 0 
K K,k 


m m'a 


[H']2Ka, PP -pKa 


14.8 INFLUENCE OF TEMPERATURE 


At temperatures of 35°C or higher the enzyme may undergo 


Rapid deactivation during the course of a kinetic experiment and then low rate 
of transformation of a substrate is observed. 


Generally the rates of enzyme catalyzed reaction frequently pass 
through a maximum as the temperature is raised. The temperature at which the 
rate is maximum often referred to as the maximum temperature (optimum 
temperature). 


At low temperature no appreciable inactivation occurs or by making 
correction for the inactivation, it is possible to determine the effect of 
temperature on enzyme catalyzed reaction itself. 

„— BIES] 
[s] + K,, 


(1) At high concentration of [s] 


r = k [Eo] 


K,=Ae T, In r 2 In k, 4 In E, 
Inr=In AE -5 
° RT 
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High the energy of activation lower the rate 


(ii) At low concentration of substrate 


r= PEIS] 
0 


kn 
kik, 
r=— 2 [E JES 
menu o ILS] 
If k1 » » Ko 
r= E qs] 
E 7 


The energy of activation is now equal to E; + E2 -E.i. 


14.9 Check Your Progress 


1. What is meant by Acid base catalyst? 
2. Explain PH effect. 


14.10 Answers To Check Your Progress Questions 


1. In acid-base catalysis, the chemical reaction is accelerated by the 
addition of an acid or a base, and the acid or base itself is not 
consumed in the reaction. ... Proton donors and acceptors, 

i.e. acids and base may donate and accept protons in order to 
stabilize developing charges in the transition state 


2. Enzymes are affected by changes in pH. The most 
favorable pH value - the point where the enzyme is most active - is 
known as the optimum pH. This is graphically illustrated in Figure 
14. Extremely high or low pH values generally result in complete 
loss of activity for most enzymes 


zm pH 


k 
P? ÁN 
— 


14.11 Summary 


General catalysis is indicated when 
m rate of the reaction depends upon the concentration of the 
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buffer as well as the pH of the solution, it is most likely to be observed when 

the pH of the solution is near the pKa of the buffer, and the reaction 

mechanism involves proton transfer in the rate-determining step. Unit- 14 

The PH effect causes polydentate complexes to be thermodynamically more Acid Base catalysis 
stable than their monodentate counterpart 


14.12 Keywords NOTE 
Acid base catalysts: In acid-base catalysis, the chemical reaction is 

accelerated by the addition of an acid or a base, and the acid or base itself is 

not consumed in the reaction..... Proton donors and acceptors, 

i.e. acids and base may donate and accept protons in order to stabilize 

developing charges in the transition state 


PH pH is a measure of hydrogen ion concentration, a measure of the acidity 
or alkalinity of a solution. The pH scale usually ranges from 0 to 14. 
Aqueous solutions at 25°C with a pH less than 7 are acidic, while those with 
a pH greater than 7 are basic or alkaline PH 


14.13 — Self-assessment questions and exercises 

1. Discuss the influence of pH and temperature on enzyme catalyzed 
reactions. 

2. Discuss the effect of subsistent on reaction rates in terms of Hammett 
equation and Taft equation 

3. Discuss the effect of ionic strength on reaction rates. 


14.14 Further readings 
Shriver and Atkins., Physical Chemistry, 5th ed : W. H. Freeman and 
Company New York. 
C.M. Guldberg and P. Waage," Studies Concerning Affinity" Forhandlinger i 
Videnskabs-Selskabet i Christiania (1864), 35 
P. Waage, "Experiments for Determining the Affinity Law" ,Forhandlinger i 
Videnskabs-Selskabet i Christiania, (1864) 92. 
. C.M. Guldberg, "Concerning the Laws of Chemical 

Affinity", Forhandlinger i Videnskabs-Selskabet i Christiania (1864) 111 
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(init 14 M.Sc. Degree Examination 
Acid Base catalysis Advanced Physical chemistry (CBCS 2018 — 19 Academic year onwards) 
Time: Three hours Maximum: 75 marks 
NOTE SECTION A Answer All Questions 
(10 X2 = 20) 
1) What is simple definition of entropy? 
2) How is vibrational partition function calculated? 
3) What is difference between Fermi Dirac and Bose Eistein statistics? 
4) Write the wave equation formula. 
5) Explain the paulis exclusion principle. 
6) Give the two application of HMO method. 
7) Define point group. 
8) Define the reducible and irreducible representations. 
9) What is difference between primary and secondary salt effect. 
10) Mention the acid base catalysis. 
SECTION B Answer ALL questions, choosing either (a) or (b) 
(5 x5 =25) 
i, (a)Write the short notes on negative absolute temperature. 
(or) 
(b) Mention the equilibrium constant from partition function. 
12. (a) Define the approximation and variation methods of quantum 
mechanical treatments. 
(or) 


(b) How will you find out the Bose Eistein Distribution law? 
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13. (a) Write the short notes on HMO methods. 


(or) Unit — 14 
(b) What are character table of C2V and C3 V point groups. Acid Base catalysis 


14. (a) Write the applications of group theory of electronic spectra. 
(or) 


(b) What are the factors affected in reaction rate of solutions? NOTE 


15. (a) Mention the Hammett and Taft equations. 
(or) 


(b) What is bronsted relation of acid-base catalysis reactions? 


SECTION C Answer any THREE questions (3 x 10 = 30) 
16) How is Boltzmann distribution calculated? 
17) How do solve the one dimensional wave equation? 
18) (a) predict the IR and NMR spectra of H2O and NH3 molecules. 


(b) Define point groups, and mention the classification of point 
groups. 


19) Explain the Bronsted Bjerrum equation. 
20) Explain the following 

a. Zuker Hammet equation 

b. Michaelis Menton Equation 


c. Lineweaver-Burke Equation. 
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